TM No. 377 -2 

 3?00 FORTRAN (2,1) / 



St^m«miT I Np CORWFCT (JJ) — — — 



COMMOM U(?400)» W(?400) 

 RPAL LOGF ' — '"■ "~ 



WQTTP(61.inO0) 

 " — Dn 999 J-l.JJ " 



le- (W(J)) A»99.4 

 A Y«!l (JV/¥(,i^ 



YY=A8SF(Y) 



rp (YY-i.tsft49) lo.iitii - - 



10 Xa( ( (-. 71 (S?4037»YY-3. 5271 192 )»YY* 1.1 61 6579 )*YY*<f 7.905764) •YY 

 XaX-,00nie,4998S 



GO TO 9n 



11 lF(YY-?9. 743803) lll»lll»llll 

 im XaATANJF(YYl 



X^X««^7,?QB788 ^ - 



GO TO 90 

 ill X=(LOr,F(YY-1.205B0 34) ♦9, 7816200 >/. 14598791 

 90 TmFTA=X 



IP (Y> 3»99,2 

 ? IP (THFTA-45.) lft.18»19 



- ""l-R TrVTH=T^TDnv+.0 00(^'5*TTHETA')»(THFTA)^I.Tr(J ~ 



A=90,-THFTA 



iP(A-7n.r i22ti^2.i2i 



12? RxrHa-.0 0O9«(A*A) ♦ , 1 238*A* . 758 



60 TO 30 

 121 RXTH = -,0 0'^«(A»A) ♦,07l*A-l,34 

 - - -Grr-TO-^n - - - .- - . — - . 



19 If (ThFTA-70,) nw,118.119 

 11^ HvTMa-,C10no9*(THETA»TWTA) ♦ . 01238»TH€TA* . 758 

 Aago.-THFTA 



GO TO 1?0 — 



ng RyTHa(-,ooo5«(THETA»THETA) ♦.071*THETA-1 .34) 



i;jO RyTHa(-,001 ♦.0 00**5*A)»(A) *1,00 

 ?0 IP (U(J) ) 7«98.«< 

 7 TmfTA»270.-THFTA 



GO TO 9 

 ^ T4FTAr90.-THETA 



- g ttf,r)-=n(>t>/»xTii- - 



W ( ))«W(J)/PYTH 



6n TO 75 



9M IF (W(J)) 87t88.8fl 

 'arr V^TUa?7r), 



W (J) a-W(J) 

 UU )« n .- 



yo TO 25 

 8« ThFTA=90, - - - . 



W(J)aw(J) 



crrDTo. 



GO TO 25 



!r-tP — (T>1gT A-45. ) 28. ? 8t? 9 " 



?fl RYTH»(-.001f.0 00 085«THETA)»THETA*1.00 



A390,-THFTA- -■ - 



lP(A-70.) 222*222,221 

 TTf RyTffr(-. oooq 9 »C»*> »^ ,>llJ «» »». r M> 



Go TO 30 



