Kozák Miklós: A szabadfelszínű nempermanens vízmozgások számítása digitális számítógépek felhasználásával (Akadémiai Kiadó, Budapest, 1977)
Mellékletek. Programok, szubrutinok, nyomtatási sémák 313 - 4. fejezethez
04.03. A NÉGYISMERETLENES LINEÁRIS EGYENLETRENDSZER MEGOLDÁSÁNAK SZUBRUTINJA (EQSOL4) LN oooi SU3R0UTÍNE EQSOLL (A, N> LN 0002 C AZ EGYENLETRENOSZER MEGOLDÁSA LN 0003 DIMENSION A <N,5> LN 00 0L 1 = 0 LN 0005 K=2 LN 0006 1 1=1+1 LN 00 07 P=1./A(I,2> LN 00 08 R=A (Kt 2) LN 0009 DO 2 M=3,5 LN 0010 T=A <ItN) +P LN Odll A (I t M) =T LN 0012 2 A (K t M) =A (K, M) -T*R LN 0013 I=< LN 001L K=K + 1 LN 0015 P=1 t/A (I t 3) LN 0016 A (I, L ) =A (ItL) *P LN 0017 A(I t 5 )=A (I t 5)*P LN 0018 IFd.EO.NÍG OTO 3 LN 0019 R=A (K, 1) LN 0020 A (K, 2) =A (K, 20-A(I,L) «R LN 0021 A(K,5)=A(K,5)-A(I,5)*R LN 0022 K=K + 1 LN 0023 R=A (K11) LN 002L A(K,2)=A(Kt2>-A(I,L),R LN 00 25 A(K,5! = A(K,5)-A(I,5)»R LN 0026 GOTO 1 LN 0027 3 I=N-1 LN 0028 AU, 5 MACI, 5)-A(I,3)*A(N,5) LN 00 29 K=N LN 0030 L M=I-1 LN 0031 ACM, 5 MACH, 5)-A (1,50 «A(H,L> LN 0032 K=M LN 0033 H=M -1 LN 003L A(M,5)=A (Nt 50-A (1,5) «ACH,L> -A <K,5> »A (M,3> LN 00 35 I*M LN 0036 IF( M,NE, 10 GOTO h LN 0037 RETURN LN 0038 END USASI FORTRAN DIAGNOSTIC RESULTS FOR EOSOLL NO ERRORS 330
