Un método de factorización de matrices ( descomposición LU), es el método de doolitle, en el cual la matriz triangular inferior L tiene todos los elementos de su diagonal principal iguales a 1.
PSEUDOCODIGO
Leer a, b
n=longitud (a)
L=matriz identidad
Para k=1 hasta n
suma=0
Para p=1 hasta k-1
suma=suma+l(k,p)*u(p,k)
fin para
u(k,k)=(a(k,k)-suma)
Para i= k+1 hasta n
Suma=0
Para r=1 hasta k-1
suma=suma+l(i,r)*u(r,k)
Fin Para
l(i,k)=(a(i,k)-suma)/u(k,k)
Fin Para
Para j = k+1 hasta n
suma = 0
Para s = 1 hasta k-1
suma=suma+l(k,s)*u(s,j)
Fin Para
u(k,j)=(a(k,j)-suma)/l(k,k)
Fin Para
Fin Para
Z = sustitución Progresiva(L, b, n)
X = sustitución Regresiva(U, Z, n)
Imprima X
Fin Programa