PROGRAM RCY DEM C This program demonstrates the use of the coefficients calculated C by the programs trnrox. The program will produce a table of values C for different values of lambda as set in the DATA statements. C DOUBLE PRECISION DLAMDA(7),RCY,RCY11,RCY12,RCY21,RCY22 INTEGER J DATA DLAMDA/1D0,.5D0,.25D0,.2D0,.125D0,.1D0,2D0/ WRITE(6,5) 5 FORMAT(1X, &'RESISTANCE FUNCTIONS RCY from Jeffrey & Onishi, J.F.M., 1982') WRITE(6,7) 7 FORMAT(1X,'vol 139, p261') WRITE(6,8) 8 FORMAT(/1X,'Table of Functions RCY for 0.1 and g<5> in the paper. C The g<5> seems to have a misprint of 4 instead of 2. G2=1D0*XL**2/(10D0*XLP1) G3=(43D0*XL-24D0*XL**2+43D0*XL**3)/(500D0*XLP1) ISTEP=1 ELSE G2=1D0/(10D0*XL*XLP1) G3=(43D0-24D0*XL+43D0*XL**2)/(500D0*XL**2*XLP1) ISTEP=-1 ENDIF RCY=2D0*(G2*DLOG(2D0)-G3) MMIN=1 MMAX=MAXS-1 ENDIF G2SUM=2D0*G2 G3SUM=4D0*G3 DO 100 M=MMIN,MMAX,2 IF (ISTEP.GT.0) THEN NM = (M*(M+1))/2 ELSE NM = (M*(M+1))/2 +M ENDIF ZL=XLP1**(-M) COEF=0D0 DO 50 IQ=0,M COEF=COEF+RYCCOF(NM)*ZL NM=NM+ISTEP ZL=ZL*XL 50 CONTINUE RCY=RCY+COEF-(G2SUM-G3SUM/DBLE(M+2))/DBLE(M) 100 CONTINUE IF(I.EQ.J)RETURN RCY=RCY*8D0/XLP1**3 IF(I.EQ.2)RCY=RCY*XL**3 RETURN END