34
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/
DATA GD/0D0,.0793D0,.141239D0,.281835D0,.256692D0,.027815D0,
*.189065D0,.0885D0,.366911D0,.332267D0,.432254D0,.512507D0,
*.576242D0,.648601D0,.716574D0,0D0,.034369D0,.038953D0,.021D0/
DATA QD/6*0D0,.69D0,12*0D0/,FD/16*0D0,1D0,2*0D0/
DATA AN/.1538326D0,1.341953D0,-2.998583D0,-.04831228D0,
*.3757965D0,-1.589575D0,-.05358847D0,2.29129D-9,1576724D0,
*-.4363864D0,-.04408159D0,-.003433888D0,.03205905D0,.02487355D0, -
*.07332279D0,-.001600573D0,.6424706D0,-.4162601D0,-.06689957D0,
*.2791795D0,-.6966051D0,-.002860589D0,-.008098836D0,3.150547D0,
*.007224479D0,-.7057529D0,.5349792D0,-.07931491D0-1.418465D0,
*-5.99905D-17,.1058402D0,.03431729D0,-.007022847D0,.02495587D0,
*.04296818D0,.7465453D0,-.2919613D0,7.294616D0,-9.936757D0,
*-.005399808D0,-.2432567D0,.04987016D0,.003733797D0,1.874951D0,
*.002168144D0,-.6587164D0,.000205518D0,.009776195D0,-.02048708D0,
*.01557322D0,.006862415D0,-.001226752D0,.002850906D0/