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
(Folytatás) 05.02. EGYENLETÁLLANDÓK SZÁMÍTÁSÁNAK SZUBRUTINJA, FELDENAM!KUS MÓDSZERHEZ (EQCONSTP) LN 00 65 A(M-i,7)=ü. LN 0066 A (H,7)=0. L N 0067 77 CONTINUE LN 0068 M=H ♦ 1 LN 0069 A (M, 1» =0 . LN 0070 A(H,5)=0. L N 0071 A (HU, 1) =0. LN 0072 A (MH,4> =0. LN 0073 A(H+1,5> =0. LN 0074 A (M,4)=2*0X LN 0075 A(M , 6) -1• LN 0076 A (M , 2>=-l. LN 0077 A(M,3>=B2*W*B(I»16> LN 0078 A (M«7) =B 2+W +8 (I,16) LN 0079 A (M,8)=A (M,3>* (B (I,6>»B (IU ,6mB(I, 5) -8 (1+1,5) LN 0080 M=H + i LN 0081 AC= 0 LN 0082 BC=0 LN 0883 DC* 0 LN 0084 FC= 0 LN 0085 P*0 LN 0086 R=0 L N 0087 00 80 K=P,7,7 LN 0088 H1= 5 +R LN 0089 M2= 13 + K LN 0090 AC= AC + B(I,Ml)/8 (I,H2)**2+B( I+l.Ml) /B (I+1,H2) **2 LN 0091 BC= BC + 1. /8(I,M2>+l./B(I + i,MZ> LN 0092 DC=DC + B(I, Ml)+*2/B(I,HZ)**3*B(IH,Ml>**2/8<I+l,M2)*«3 LN 0093 FC= FC + B (I + 1 ,M2) -B(I , M2 ) L N 0094 8 0 R=(f LN 0095 IF( B (1,7» .GT. 0. ) AC=AC/2. LN 0096 C1=ABS(B(I,5))/B(I,15) LN 00 97 C2=ABS(BII*1,5) )/B(I*l,15) L N 0098 C3=ABS(B(I,9))/B(I,22> LN 0099 C4=ABS(B(If 1,9) )/B(I*l,22) 336
