5 * Revision 1.1.1.1 1996/04/01 15:01:57 mclareni
10 #if defined(CERNLIB_SINGLE)
11 FUNCTION C309R5(X,ETA,ZL,EPS,FCL,TPK1,ETANE0,LIMIT,ERR,NFP,
14 C *** Evaluate CF1 = F'(ZL,ETA,X)/F(ZL,ETA,X) (COMPLEX)
16 IMPLICIT COMPLEX(A-H,O-Z)
18 REAL EPS,ERR,FPMIN,FPMAX,ABSC,SMALL,PX
21 ABSC(W)=ABS(REAL(W))+ABS(AIMAG(W))
29 IF(ABSC(F) .LT. FPMIN) F=FPMIN
35 C *** begin CF1 loop on PK = k = lambda + 1
48 IF(ABSC(C) .LT. FPMIN) C=FPMIN
49 IF(ABSC(D) .LT. FPMIN) D=FPMIN
54 IF(ABSC(FCL) .LT. SMALL) FCL=FCL/SMALL
55 IF(ABSC(FCL) .GT. FPMAX) FCL=FCL*FPMIN
57 IF(REAL(PK) .LE. PX) THEN
58 IF(ABSC(DF-1) .GE. EPS) GO TO 10
60 ERR=EPS*SQRT(REAL(NFP))
63 IF(LPR) WRITE (6,1000) LIMIT,ABS(X)
67 1000 FORMAT(1X,'***** CERN C309 CCLBES ... CF1 (COMPLEX) HAS FAILED ',
68 1'TO CONVERGE AFTER',I10,' ITERATIONS AS ABS(X) =',F15.0)