RUNT VERSION FEB Ti B 17:12 Oit/Zi/lb 



FUNCTION LN6QF(«,N,Nl,A,9»0TRMNT,Z) 

 C SILVES SIMULTANEQLIS LINEAR EQUATIONS Rr r^USSTAN REDUCTION. 

 C. FORTRAN IV EQUIVALENT OF LNEQS. 



12 REAL A(M»M),8(1»1) . Z ( M ) • DTRMNTt RM A¥f RMC XT, W. OHV 



12 NMi«N-l 



1* DO 40 J«1,NM1 



15 J1»J+1 



C. FIND ELEMENT OF COL J, ROWS J-N, WHtr'H HAS -^AX ABSOLUTE VALUF, 



17 LMAX.J 



20 RKAX-ABS( »( J> J) ) 



34 DO 6 K'Jl.N 



35 RNEXT»ABS(A(K,J) ) 



52 IF (RMAX .GE. RNEXT) GO TO 6 



55 RMAX-RNEXT 



57 L^'AX«K 



60 6 CONTINUE 



63 IF (LMAX ,N6. J) GO TO 10 



C MAX ELEMENT IN COLUMN IS ON DIAGONAL 



65 IF (A(J>J)> 20. '44. 20 



C. MAX ELtMENT IS NOT JN DIAGONAL. EXCHtNGF ROUS J AND LMAX. 



73 10 DO 12 L'J.N 



75 W'A(J,L) 



102 A( J>L>«A(LMAX>L) 



113 12 A(LMAX,L)«W 



124 DO 14 L=1»N1 



125 W=6(J,L) 



132 B( J,L)'B(LMAX.l ) 



143 14 B(LMAX,L)=W 



154 DTRMNT • -OTRINT 



C. ZERO COLUMN J BcLOV Trl£ DIAGONAL. 



155 20 Z(J>»1./A(J»J) 

 165 DO 30 K=Jl,N 



167 IF (A(K,J)) 22.30,22 



175 22 W — Z( J)*A(K. J) 



205 DO 24 L'Jt.N 



2U7 24 A(K.L)«W*A( J.L)*A(K»L) 



230 00 26 L'l.Nl 



231 26 B(K.L)«W*9( J.L)*B(K.L) 

 252 30 CONTINUE 



255 40 CONTINUE 



257 IF (A(N,N)) 42,^4,42 



3b5 42 Z(N)«i./A(N.N) 



C OBTAIN SOLUTION SY iACK SUBSTITUTION. 



275 DO 50 L»1,N1 



277 50 6(N,L)»Z(N)*8(N,L) 



315 DO 60 K-1,NM1 



316 J*N-K 



317 Jl'J+1 



321 DO 58 L'l.Nl 



322 W-O. 



323 DO 56 I'Jl.N 



325 56 V>«A( J. I)*B( I.L)4-W 



342 58 B( J,L)«(B{J.L)-»J)*Z( J) 



362 60 CONTINUE 



Table D-2. Continued 



129 



