function x = gauss_elim_vector(A,b) % A should be square % size(A) = [length(x) length(x)] n = size(A,1); for k = 1:(n-1) for i = (k+1):n mult = A(i,k) / A(k,k); A(i,k:n) = A(i,k:n) - mult * A(k,k:n); b(i) = b(i) - mult * b(k); end end x = zeros(n,1); for k = n:-1:1 sum = b(k) - A(k,(k+1):n) * x((k+1):n); x(k) = sum / A(k,k); end