%Método de Factorización Triangular "LU"
%Felipe Malavet Ardila, Juan Felipe Rodriguez Navarro, Natalia Romero
Ortiz
%Grupo 2
clear,clc
disp ('Bienvenido al Matriz-Solver LU.')
%Se ingresa el tamaño del sistema de ecuaciones a solucionar
n = input ('Ingrese el tamaño de la matriz de coeficientes "A": ');
k=2;
q=0;
p=0;
w=0;
e=0;
%Se genera una matriz auxiliar para evitar indeterminaciones
AU=transpose(zeros(n,1));
%Se genera una matriz cualquiera como "plantilla"
A=magic(n);
%Se genera un vector de ceros del tamaño del vector solución del
sistema de ecuaciones
B=zeros(n,1);
%Se genera una matriz identidad del mismo tamaño de "A"
L=eye(size(A));
%Se generan dos vectores "plantilla" para "X" y "Y" de ceros
Y=zeros(n,1);
X=zeros(n,1);
%Dentro del siguiente código se pide al usuario que ingrese, dato por
dato,
%los componentes de cada una de las filas de la matriz "A" de
coeficientes
%del sistema
for (i=1:n)
for (j=1:n)
A(i,j)=input ('Ingrese el elemento i,j de la matriz de
coeficientes "A" (Filas:) ');
end
end
%
for (j=1:n)
for (i=k:n)
if(A(j,j)==0)
for(o=j+1:n)
while(A(o,j)~=0)
AU(1,:)=A(j,:);
A(j,:)=A(o,:);
A(o,:)=AU(1,:);
m= A(i,j)/A(j,j);
A(i,:)=-A(j,:)*m+A(i,:);
L(i,j)=m;
end
end
end
m= A(i,j)/A(j,j);
