LU decomposition: Difference between revisions

Content added Content deleted
Line 2,169: Line 2,169:
<lang parigp>matlup(M) =
<lang parigp>matlup(M) =
{
{
my (P = matid(#M), L = P, U = M);
my (L = matid(#M), U = M, P = L);


for (i = 1, #M-1, \\ pivoting
for (i = 1, #M-1, \\ pivoting