12 PDF
12 PDF
m 1 of 1
function x = Gauss(A,b)
n = size(A,2);
%FORWARD ELIMINATION
for k = 1:n-1
for i = k+1:n
factor = A(i,k)/A(k,k);
for j = k+1:n
A(i,k)=0;
A(i,j) = A(i,j) - factor*A(k,j);
end
b(i,1) = b(i,1)- factor*b(k,1);
end
end
%BACK SUBSTITUTION
x(n,1) = b(n,1)/A(n,n);
for j = n-1:-1:1
sum = 0;
for i = n:-1:j+1
sum = sum + A(j,i)*x(i,1);
end
x(j,1) = (b(j,1)- sum)/A(j,j);
end
end