END IF END Decompose big SUB Pivot (a, o, s, n. k) p=k ABS(2o(k),k/So(k)) DOFOR ii = k + 1, n dummy = ABS(ac(ii).k/So(ii)) IF dummy > big THEN dummy p =ii END IF END DO dummy Op Op = OK Ok = dummy END Pivot big SUB Ludecomp (a, b, n, tol, x, er) DIM On, Sn er = 0 CALL Decompose(a, n, tol, o, s, er) IF er -1 THEN CALL Substitute(a, o, n, b, x) END IF END Ludecomp SUB Decompose (a, n, tol, o, s, er) DOFOR 1 - 1. n 0; = i Si = ABS(1,1) DOFOR j = 2, n IF ABS(ai,j)>S; THEN S; = ABS(aj) END DO END DO DOFOR k = 1 n - 1 CALL Pivot(a, o, s, n. k) IF ABS(Bork).)) big THEN dummy p =ii END IF END DO dummy Op Op = OK Ok = dummy END Pivot big SUB Ludecomp (a, b, n, tol, x, er) DIM On, Sn er = 0 CALL Decompose(a, n, tol, o, s, er) IF er -1 THEN CALL Substitute(a, o, n, b, x) END IF END Ludecomp SUB Decompose (a, n, tol, o, s, er) DOFOR 1 - 1. n 0; = i Si = ABS(1,1) DOFOR j = 2, n IF ABS(ai,j)>S; THEN S; = ABS(aj) END DO END DO DOFOR k = 1 n - 1 CALL Pivot(a, o, s, n. k) IF ABS(Bork).))