34
DR=FX/F
X=X+DR
IF(ITER.GT.10) GO TO 4
ITER=ITER+1
IF(DABS(DR/X).GT.1.D-6) GO TO 1
4 CALL COMPL1(X)
RETURN
END
SUBROUTINE PCONT1(P,T)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/Z/Z
Z=-1D0
IF(T.LT.250D0.OR.T.GT.340D0)Z=0D0
IF(P.LE.0D0.OR.P.GT.12D0) Z=0D0
RETURN
END
SUBROUTINE COMPL1(RO)
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 KM
INTEGER BN,CN
COMMON/KM/KM/COEF2/B,C(53)/BCKN/BN(53),CN(53),KN(53)/AI1/AO,A1
ROR=KM*RO
S1=0D0
S2=0D0
S3=0D0
DO 1 N=8,53
EXP=DEXP(-CN(N)*ROR**KN(N))
IF(N.LE.13) S1=S1+C(N)
S2=S2+C(N)*(BN(N)-CN(N)*KN(N)*ROR**KN(N))*ROR**BN(N)*EXP 1
S3=S3+C(N)*(-CN(N)*KN(N)**2*KM*ROR**(KN(N)-1)*ROR**BN(N)*
*EXP+(BN(N)-CN(N)*KN(N)*ROR**KN(N))*BN(N)*KM*ROR**(BN(N)-1)*
*EXP-(BN(N)-CN(N)*KN(N)*ROR**KN(N))*ROR**BN(N)*EXP*CN(N)*KN(N)*
*KM*ROR**(KN(N)-1))AO1=B*RO-ROR*S1
AO=AO1+S2
A1=AO+AO1+RO*S3
RETURN
END
BLOCK DATA DCAGA8
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 KD
INTEGER BN,CN,GN,QN,FN
COMMON/EFD/ED(19),KD(19),GD(19),QD(19),FD(19)
*/BCKN/BN(53),CN(53),KN(53)/UN/UN(53)
*/AN/AN(53)/GQFN/GN(53),QN(53),FN(53)
DATA ED/151.3183D0,244.1667D0,298.1183D0,337.6389D0,324.0689D0,
*99.73778D0,241.9606D0,296.355D0,370.6823D0,365.5999D0,
*402.8429D0,427.5391D0,450.6472D0,472.1194D0:488.7633D0,
*2.610111D0,26.95794D0,105.5348D0,122.7667D0/
DATA KD/.4619255D0,.5279209D0,.583749D0,.6341423D0,.6406937D0,
*.4479153D0,.4557489D0,.4618263D0,.6798307D0,.6738577D0,
*.7139987D0,.7503628D0,.7851933D0,.8157596D0,.8389542D0,
*.3589888D0,.3514916D0,.4533894D0,.4186954D0/