1 PROGRAM CFC 



2 C *• THIS PROGRAM COMPUTES A SET OF SERIES COEFFICIENTS AND THEN 



3 C »* COMPUTES THE CORRESPONDING CONTINUED FRACTION COEFFICIENTS. 



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



5 DIMENSION C0EF(21 ,23,3), CHECK(20), C(82), S(10), A(20). B(20) 



6 C(1) = 1. 



7 BOTTOM = 1 . 



8 TOP = 1 . 



9 DO 2 I = 1 ,45 



10 X = 48 * I 



11 Y=9*(I+I-1 )**2 - 4 



12 C(I+1) = C(I) * Y / X 



13 2 CONTINUE 



14 PRINT 20, (C( I ) , I • 1 ,40) 



15 20 FORMAT (5G20.9) 



16 11 FORMAT (/) 



17 DO 100 I = 1 , 11 



18 COEF(I ,1,3) = 0. 



19 COEFd , 1 + 1 ,3) =0. 



20 COEFd , 1+2,3) = 0. 



21 100 CONTINUE 



22 A(1 ) = C(2) 



23 C0EF(2,2,3) = 1. 



24 DO 140 I = 3,21 



25 DO 1 10 J = 2, I 



26 COEFd, J, 1) = C0EF(I-1 ,J,3) 



27 COEFd, J, 2) = C0EF(I-2,J,3) 



28 COEFd, J, 3) = C0EF(I-1 .J-1 ,3) 



29 110 CONTINUE 



30 IF (I .EO. 3) GO TO 150 



31 CON = 0. 

 32 

 33 

 34 

 35 

 36 

 37 

 38 

 39 



40 1 

 41 

 42 

 43 



44 1 

 45 

 46 

 47 

 48 

 49 

 50 

 51 

 52 



53 BT = C(K) * COEFd ,0-1 ,1 ) 



54 130 CONTINUE 



55 PRINT 160, CON, AT, BT 



56 PRINT 11 



AT = 0. 











BT = 0. 











K = I - 3 











DO 120 J = 



3,1 









K = K + 1 











CON = C(K) 



* COEFd 



,J-' 



1,3) 



+ CON 



AT = C(K) 



♦ COEFd , 



J-1 



.2) ■ 



f AT 



BT = C(K) 



* COEFd , 



J-1 . 



,1) ■ 



f BT 



CONTINUE 











PRINT 160, 



CON, AT. 



BT 







CHECK(I-2) 



= BT 









A(I-2) = - 



(CON + C(K+1)) / 



AT 



CONTINUE 











CON = 0. 











AT = 0. 











BT = 0. 











K = I - 2 











DO 130 J = 



3, I 









K = K + 1 











CON = C(K) 



* COEF( I 



,J-' 



1,3) 



+ CON 



AT = C(K) 



* COEFd , 



J-1 



.2) 



+ AT 



91 



