na-lab-lu

code

clc;

A=input("enter matrix");

B=input("enter rhs");

n=length(B);

m=eye(n);

x=zeros(n,1);

% i=column j=row

for i=1:n-1

for j=i+1:n

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

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

end

end

L=m;

U=A;

y=inv(L)*B;

x=inv(U)*y;

disp(x);