programa LU

function lu()
n=input('introduzca la cantidad de ecu ');
for i=1:n
    for j=1:n
        A(i,j)=input('los elementos de la matriz A ');
    end
end
for i=1:n
    b(i)=input('los elementos del vector B ');
end
U=zeros(n);
L=zeros(n);
for j=1:n
    L(j,j)=1;
end
for j=1:n
    U(1,j)=A(1,j);
end
for i=2:n
    for j=1:n
        for k=1:i-1
            s1=0;
            if k==1
                s1=0;
            else
                for p=1:k-1
                s1=s1+L(i,p)*U(p,k);
                end
            end
            L(i,k)=(A(i,k)-s1)/U(k,k);
        end
        for k=i:n
             s2=0;
           for p=1:i-1
               s2=s2+L(i,p)*U(p,k);
           end
           U(i,k)=A(i,k)-s2;
        end
    end
end
disp(A);
disp(L);
disp(U);
c=inv(L)*b';
x=inv(U)*c;
disp(x);

y se debe hacer correr es de la siguiente forma:

lu()
introduzca la cantidad de ecu 3
los elementos de la matriz A 4
los elementos de la matriz A -2
los elementos de la matriz A -1
los elementos de la matriz A 5
los elementos de la matriz A 1
los elementos de la matriz A -1
los elementos de la matriz A 1
los elementos de la matriz A 2
los elementos de la matriz A -1
los elementos del vector B 9
los elementos del vector B 7
los elementos del vector B 12


     4    -2    -1
     5     1    -1
     1     2    -1


    1.0000         0         0
    1.2500    1.0000         0
    0.2500    0.7143    1.0000


    4.0000   -2.0000   -1.0000
         0    3.5000    0.2500
         0         0   -0.9286


   -1.3077
   -0.2308
  -13.7692