Algorithme factorisation LU

Bonjour,
j ai ce algothime je veux le faire en pascal ou en c
poir j de 1 a n faire
pour i de 1 a j faire // calcul des rij matrice U
rij =aij

pour k de 1 a i-1 faire
rij=rij-l jk rkj
fin pour
fin pour
pour i de j+1 a n faire // calcule lij matrice L
lij = aij
pour k de 1 a j-1 faire
lij = lij-lik rkj
fin pour
lij= lij/rij
finpour
fin pour

1 0 0 r r r a a a
r 1 0 0 r r a a a
l r r 1 * U 0 0 r = a a a a





svl aide mpi