LU decomposition: Difference between revisions

m
(→‎{{header|Tcl}}: Added zkl)
Line 2,469:
fcn lu(A){
n:=A.len();
L:=make_arrayeye(n,n,0.0);
U:=make_array(n,n,0.0);
P:=pivotize(A);
Line 2,475:
 
foreach j in (n){
L[j][j]=1.0;
foreach i in (j+1){
U[i][j]=A[i][j] - (i).reduce('wrap(s,k){ s + U[k][j]*L[i][k] },0.0);
(i).reduce('wrap(s,k){ s + U[k][j]*L[i][k] },0.0);
}
foreach i in ([j..n-1]){
Anonymous user