FUNCTION AIGAU15(XT) IMPLICIT REAL*8 (A-H,O-Z) IF(XT.LT.-27)THEN AIGAU15=0 RETURN ENDIF IF(XT.GT.7)THEN AIGAU15=1 RETURN ENDIF Z=-ABS(XT) IF(Z.EQ.0)THEN AIGAU15=0.5D0 ERT=0 RETURN ENDIF IF(Z.LT.-3)THEN ALGAU15=-Z*Z-LOG(-2*Z)-.5723649429247001D0 AIGAU15=EXP(ALGAU15)*GLAGU(Z) ELSE AIGAU15=AIGAUP(Z) ENDIF IF(XT.GT.0)THEN AIGAU15=1-AIGAU15 ENDIF RETURN END FUNCTION GLAGU(Z) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION XI2(17),AI2(17) C THE LAST THREE POINTS OF THE 20 HAVE BEEN REMOVED AS FOR C THIS FUNCTION THEY CONTRIBUTE NOTHING DATA (XI2(I),AI2(I),I=1,17)/.7053988969198875D-1, 2 .1687468018511139D0,.3721268180016114D0,.2912543620060683D0, 3 .9165821024832736D0,.2666861028670013D0,1.707306531028344D0, 4 .1660024532695068D0,2.749199255309432D0,.7482606466879237D-1, 5 4.048925313850887D0,.2496441730928322D-1,5.615174970861617D0, 6 .6202550844572237D-2,7.459017453671063D0,.1144962386476908D-2, 7 9.594392869581097D0,.1557417730278120D-3,12.03880254696432D0, 8 .1540144086522492D-4,14.81429344263074D0,.1086486366517982D-5, 9 17.94889552051938D0,.5330120909556715D-7,21.47878824028501D0, A .1757981179050582D-8,25.45170279318691D0,.3725502402512321D-10, B 29.93255463170061D0,.4767529251578191D-12,35.01343424047900D0, C .3372844243362438D-14,40.83305705672857D0,.1155014339500399D-16/ GLAGU=0 DO I=1,17 ARG=XI2(I)/(2*Z) GLAGU=GLAGU+AI2(I)*EXP(-ARG*ARG) ENDDO RETURN END FUNCTION AIGAUP(Z) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION ALGVAL(9) DATA ALGVAL/0.110452484992927D-04,0.102687868070607D-03, 2 0.731358293340571D-03,0.400497116494002D-02, 3 0.169474267623447D-01,0.558058841491461D-01, 4 0.144422183173242D+00,0.297941545282589D+00, 5 0.500000000000000D+00/ I=(Z+3)*8/3+.5D0+1 X=(I-1)*3D0/8-3 CALL GLEG7(X,Z,AI) AIGAUP=ALGVAL(I)+AI*0.5641895835477563D0 RETURN END SUBROUTINE GLEG7(ABEG,AEND,AI) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION X7(4),A7(4) DATA X7/0.9491079123427585D0,0.7415311855993944D0, 2 0.4058451513773972D0,0D0/ DATA A7/0.1294849661688697D0,0.2797053914892767D0, 2 0.3818300505051189D0,0.4179591836734694D0/ XMID=(ABEG+AEND)/2 RAT=(AEND-ABEG)/2 AI7=A7(4)*EXP(-XMID*XMID) DO I=1,3 XL=XMID-X7(I)*RAT XH=XMID+X7(I)*RAT AI7=AI7+A7(I)*(EXP(-XL*XL)+EXP(-XH*XH)) ENDDO AI=AI7*(AEND-ABEG)/2 RETURN END