PI=3.1415926 T=.094 R=.007 Z=0.0 TOL=.1E-5 DO 5 I=1,5 D=.0005*FLOAT(I) II=0 K=1 GM=PI/2.1 3 IF(GM.GT..5*PI) GM=GM-.5*PI SN=SIN(GM) CS=COS(GM) IF(ABS(CS).LT..1E-5) GO TO 4 SIG=(1.0-SN-D*SN/(2.*R))*2.*R/CS PS=(180./PI)*ATAN(SIG/T) CALL ITER8(I,K,GM,XM1,SIG,YM1,Z,DX,II,TOL) IF(II) 4,4,1 4 GM=GM+DX GO TO 3 1 GAM=GM*180./PI 5 WRITE(5,2) I,D,GAM,SIG,PS 2 FORMAT(/,I2,3X,'D ',F10.4,3X,'GAM ',F10.4,3X,'SIG ' 2,F10.4,3X,'PSI ',F10.4) END