1 PROGRAM PWRTRN 



2 C ** THIS PROGRAM DETERMINS TRUNCATION POINTS FOR THE POWER SERIES. 



3 IMPLICIT DOUBLE PR ECISION ( A-H, 0-Z) 



4 DIMENSION D(50), ALOGD(50) 



5 D(1 ) = 1 . 



6 ALOGD( 1 ) = 0. 



7 P = 3. 



8 DO 50 I = 2, 50 



9 D(I) = D(I-1 ) / P / (P-2) 



10 P = P + 3. 



11 IF (D( I ) .LE 0.) GO TO 50 



12 ALOGD(I) = ALOG10(D(I)) 



13 50 CONTINUE 



14 PRINT 60, D, ALOGD 



15 60 FORMAT (10E12.6) 



16 DH = 18. 



17 M = 1 



18 DO 10 K = 2,50 



19 30 P = M - K 



20 Z = (ALOGD(K) - ALOGD(M) + DH) / 3. / P 



21 IF (P .GT. -1 .1 ) GO TO 20 



22 A = ALOGD(M) - ALQGD(M+1) - 3. * Z 

 23 

 24 

 25 

 26 

 27 

 28 

 29 

 30 

 31 

 32 

 33 





IF (A .GT. 0. ) GO TO 20 





M = M + 1 





GO TO 30 



20 



L = K - 1 





MM = M - 1 





AZ = EXP (Z * 2.3025851) 





AZSQ = AZ ♦ AZ 





PRINT 40, L, MM, Z, AZ , AZSQ 



40 



FORMAT (215, 4E15.8) 



10 



CONTINUE 





END 



90 



