5 * Revision 1.1.1.1 1996/04/01 15:02:28 mclareni
12 DIMENSION C(0:128),CC(0:128),AC(0:64),BC(0:64)
14 PARAMETER (Z1 = 1, Z2 = 2, HF = Z1/2)
31 2 C(J)=F(ALFA*COS((J+HF)*C2)+BETA)
33 C COMPUTE B-COEFFICIENTS
47 3 BC(J)=C1*(A0*COS(F1+F2)-A1*COS(F2))
50 C COMPUTE NEW C-COEFFICIENTS
54 5 C(I1+J)=HF*(AC(I1-J)-BC(I1-J))
56 C CHECK IF CALCULATION IS FINISHED
70 IF(A0 .LE. EPS .AND. A1. LE. EPS .AND. A2 .LE. EPS) GO TO 9
76 C DOUBLE THE NUMBER OF COEFFICIENTS.
85 C REQUIRED ACCURACY NOT OBTAINED
90 13 DELTA=DELTA+ABS(C(J))
92 CALL MTLPRT(NAME,'E406.1','REQUIRED ACCURACY NOT OBTAINED')
95 C REQUIRED ACCURACY OBTAINED
96 C SUM NEGLECTED TERMS IN EXPANSION
102 C CHECK IF FURTHER REDUCTION OF COEFFICIENTS IS POSSIBLE.
107 8 IF(CC(NC) .GE. REST) GO TO 14