Metodo de Doolittle
disp('Metodo de Doolittle');
Metodo de Doolittle
A=input('Ingrese la matriz A en []: ');
b=input('Ingrese el vector b en []: ');
n=size(A);%Tamaño de A
L=eye(n);%Matriz de 1
U=eye(n);%Matriz de 1
cond(A)%verificar si esta bien condicionada la matriz
ans =
4.5228
% halla la fila 1 de U
for j = 1:n
U(1,j) = A(1,j);
end
Warning: Colon operands must be real scalars. This warning will become an error in a future release.
% halla columna 1 de L
for i = 1:n
L(i,1) = A(i,1)/U(1,1);
end
Warning: Colon operands must be real scalars. This warning will become an error in a future release.
for k=2:n % Se hallan las demas filas desde la 2
for j=k:n
Acumulador2 = 0; %Se inicializa el acumulador en cero para el
despeje de U
for p=1:k-1
Acumulador2 = Acumulador2 + L(k,p)*U(p,j); % Acumulador es como
el multiplicador
end
U(k,j) = (A(k,j) - Acumulador2);
end
for i= k:n
Acumulador = 0; %Se inicializa el acumulador en cero para el despeje
de L
for p=1:k-1
Acumulador = Acumulador + L(i,p)*U(p,k);
end
L(i,k) = (A(i,k) - Acumulador)/U(k,k);
end
end
1
