202LUFactorisation

LU Factorisation

%LU

clear all;

clc;

a = input("enter matrix: ")

n = length(a);

l = eye(n);%identity matrix

for i=1:n-1

for j = i+1:n

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

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

for k=1:n

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

end

end

end

l

a