na-lab-gausselimination

code

clc;

A=input("enter coefficient matrix");

B=input("Enter RHS");

Aug=[A B];

n=length(B);

x=zeros(n,1);

% i=column j= row

for i=1:n-1

for j=i+1:n

m=Aug(j,i)/Aug(i,i);

Aug(j,:)=Aug(j,:)-m*Aug(i,:);

end

end

x(n)=Aug(n,n+1)/Aug(n,n);

for i=n-1:-1:1

x(i)=(Aug(i,n+1)-Aug(i,i+1:n)*x(i+1:n))/(Aug(i,i));

end

disp(x);