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 - 5. fejezethez

05.03.

LN 0001 LN 000? LN 0003 LN 00 01* LN 00 0? LN 0006 LN 00 07 LN 00 08 LN 0009 LN 0010 LN 0011 LN 0012 LN 0013 LN 0011* LN 0015 LN 0016 LN 0017 LN 0018 LN 0019 LN 0020 LN 00 21 LN 00 22 LN 0023 LN 0021« LN 0025 LN 0026 LN 00 27 LN 0028 LN 00 29 LN 0030 LN 0031 LN 0032 LN 0033 LN 0031« LN 0035 LN 0036 LN 00 37 LN 0038 LN 0039 LN 00<*0 LN 001*1 LN 00<«2 LN 001*3 LN 00<*1|

A HÉTISMERETLEN ES LINEÁRIS EGYENLETRENDSZER MEGOLDÁSÁNAK SZUBRUTINJA (EQSOL 7)

SUBROUTINE EQUS0L7 (A ,SZ, L2) 05.03

C A LINEARIS EGYENLETRENDSZER MEGOLOASA 7 C SOLUTION OF LINEAR EGUTICNS 7 INTEGER SZ,SZM1,S.P.T DIMENSION A (L2* 8)

K=0

SZH1 = SZ-1 DO 135 1 = 1, SZM1 DO 135 M=3,6 K=K +1

C1=1./A(K*M)

N 2= M *1

00 105 S=N2,8 105 A(K,S)=A(K,S)*C1

IF(M.LT.5> GO TO 107 GO TO 11<*

107 N3= 6-M

DO 112 S = i, N3 Ml= K*S C2= A (Ml, M)

K9=H*1

00 111 P=K9,8

111 A(M1,P) = A(M1,P>-C2*A(K*P)

112 COMTINUE 11A CONTINUE

IFTM.E0.5) GO TO 115 GO TO 126 115 C2= A ( K+l * H)

00 117 S = 6, 8

117 A(Kfl,S> =A( K«-i,S)-C2«A(K,S>

IF (I.LT.SZM1) GO TO 119 GO TO 126 119 DO 12A S=2, 5 M 2= K+S C2=A (M2, 1)

OO 123 P=2, 3 T =P *A

123 A(M2,P) = A(M2,P)-C2*A(K,T)

P=8

T=8

A (M2, P> = A(M2,P) -C2*A (K,T)

12A CONTINUE

126 IFÍM.LE.5) GO TO 135

127 IFd.LT.SZMDGO TO 128

338