SUHBOUTTNF BKGStEND.0FLT.NInLfT,UINFl.O,ITAaLtiT) 

 C ROUTINE TO SOLVfc 4 SfT OF S I MUT 4l-.t OUS DIFFtRENTIAL FQUATIONS 

 C ADAPTFD F«OM ScIfcNTTFIC SUBROUTINE PaCKaGEi IBM, 1970 



Cr)NMoN/NU^M/V(5).'^ERyCb)iXiMTiIWT,ztTA,HS 



CnM'10M/NUMU/HNK(3m) 



DIMfNSinN Aiix(e.5),4{8).B(8)fC(8),PHMT(5)iAMlNHJ) 



MOlMBMIMtT+l 

 PRHT( 1) ■ 1 • 



RRMT(2)=tNn 

 PRr<T( 1 

 PPHTC' 

 IF(T.f 



(T.Lt.i 



t"LT 



,1 



hOOO 



PELTHb3600, 

 10.) DFlTH»T/9, 

 no (12^ jNil .NINUET 



U22 Uf 



y( 



J N ) » , 1 



Hv( jN)=p,noi 



Nnl 

 KYt 



)«1.0.FLDAT(MNUtT)»0.00l 



MDIM 

 f),066h*6667»riERY(I) 



xtNniPHMTt2) 



■JNTI^ 

 t£)EC 



(ill) 



(«l = r 

 C11 = J 



[i) = l 



!: ( 4 M I N I ) 

 ^Ii-X))18.i7.2 



"Sas<)32 

 707) 07 

 16666667 



tuia 

 ■} 3 

 JXf 1 



.IX (3 

 IX(6 



i.2''?8<)32 



1 ,7U71(17 



0.5 



Isl .NDiM 



• n»y(n 



nePfRY(I) 



n»o. 



H,XtN0)»H)7t6»S 



b CONTINUE 



INLET 



17b 



INLET 



176 



INLET 



177 



INLET 



176 



INLET 



179 



INLET 



180 



INLET 



181 



INLET 



162 



INLET 



163 



INLET 



leu 



INLET 



185 



INLET 



166 



INLET 



167 



INLET 



168 



INLET 



189 



INLET 



190 



INLET 



I'l 



INLET 



192 



INLET 



193 



INLET 



19U 



INLET 



195 



INLET 



196 



INLET 



197 



iNLtT 



198 



INLET 



199 



INLET 



200 



INLET 



201 



Inlet 



202 



INLET 



203 



INLET 



20U 



INLET 



205 



INLET 



206 



INLET 



207 



INLET 



206 



INLET 



209 



INLET 



210 



INLET 



211 



INLET 



212 



INLET 



213 



INLET 



21U 



INLET 



215 



INLET 



216 



INLET 



217 



INLET 



218 



INLET 



219 



INLET 



220 



INLET 



221 



INLET 



222 



INLET 



223 



INLET 



22U 



INLET 



225 



INLET 



226 



INLET 



227 



35 



