matlab code

function x=GEsolver(a,b)

%Gaussian elimination code based our Textbook p.78-80

n=size(a,1);

for j = 1 : n-1

if abs(a(j,j))<eps; error('zero pivot encountered'); end

for i = j+1 : n

mult = a(i,j)/a(j,j);

for k = j+1:n

a(i,k) = a(i,k) - mult*a(j,k);

end

b(i) = b(i) - mult*b(j);

end

end

%back-ward substitution

x=zeros(n,1);

for i = n : -1 : 1

for j = i+1 : n

b(i) = b(i) - a(i,j)*x(j);

end

x(i) = b(i)/a(i,i);

end