-85- 



c 

 c 



SOBBOOTINE flDNR 



DATA 



DATA 



DATA 

 1 



DATA 

 1 



DATA 



DATA 

 1 



DATA 



DATA 

 1 



DATA 



GO TO 40 

 10 IF (X.LT. . 15) 



IP {X.LT. 1.85) 



SIGHA = -1. 



A = 2.-X 



B = X-1. 



GO TO 50 

 20 Z = U-X 



SIGHA = SIGN(1. 



Z = ABS{Z) 



GO TO 80 

 30 SIGHA = 1. 



A = X 



B = 1<.-X 



GO TO 50 



no INT = 3 

 X = 2- *P 

 GO TO 10 



50 W = SQET(-ALOG (A+A*B) ) 

 IF (W. LT. 2.5) GO TO 70 

 IF (W.LT- n. ) GO TO 60 



IS(P,Y) 



SQBT2/1. I»1«2136/ 



A 1,A2,A3/-. 575 17029, -1.89651 33, -.0 5496 2605/ 



BO, B1,B2,B3/-. 11 377303, -3. 2934740, -2. 3749959, 



-1.1875145/ 

 CO, C1,C2,C3/-. 11 466659, -.13147744, -.236820 10, 



-050739749/ 

 DO, D1,D2/-44. 279769, 21.985462,-7-5861027/ 

 E0,E1,E2,E3/-. 0566842208,. 39370209, -.3 16650 10, 



.062089629/ 

 F0,F1,F2/-6. 2667859,4. 6662627,-2-96 28832/ 

 GO, G1,G2,G3/. 000 1851 1591, -.0020281 520, 



-. 14983844,. 010786386/ 

 HO, HI, H2/. 099529751, .521 17329,-. 068883 009/ 



GO TO 30 

 GO TO 20 



-Z) 



WI = 1./H 

 SS = ((G3*iri+G2 

 SD = ((HI-t-H2)*H 

 F = » ■<■ H«(GO+S 

 GO TO 90 



)*WI+G1)«BI 

 I+H1) ♦HI+HO 

 N/SD) 



60 SB = ((E3«W+E2) *H*E1) *H 

 SD = ((H + F2)*H + F1)*ll*F0 

 F = H + H»(E0+SN/SD) 

 GO TO 90 



70 SH = {(C3»8*C2) *B*C1) ♦« 

 SD = ((B+D2)*H+D1) *W+DO 

 F = H ♦ B»(CO+SH/SD) 

 GO TO 90 



INVEBSE GAOSSIAN ENTBY 



REDOCED ABGOHENT IS IN (.85,1.), 

 OBTAIN THE TBANSFOEHED VABIABLE 



W GBEATER THAN 4., APPEOX. F BY A 

 BATIONAL FUNCTION IN 1./B 



B BETBEEN 2.5 AND 4., APPEOX. F 

 BY A RATIONAL FONCTION IN B 



B BETBEEN 1-13222 AND 2.5, APPEOX. 

 F BY A BATIONAL FaNCTION IN W 



Z BETBEEN 0. AND -85, APPBOX. F 

 BY A EATION&L FONCTION IN Z 



Z2/(B1*Z2+A2/(B2+Z2-fA3/(B3 + Z2) ) ) ) 

 FOEH THE SOLUTION BY 

 THE PROPEE SIGN 



GO TO 100 



MOLT. F BY 



HDNE0010 

 HDNB0020 

 HDNB0030 

 HDHB0040 

 HDHR0050 

 HDNE0060 

 HDNE0070 

 HDNEOOBO 

 HDNE0090 

 HDNR0100 

 HDNE0110 

 nDNE0120 

 HDHBOnO 

 HDNE0140 

 HDHB0150 

 HDNB0160 

 SDNR0170 

 HDNB01B0 

 HDNE0190 

 nDNB0200 

 HDNB0210 

 HDNB0220 

 HDNE0230 

 tlDNR0240 

 HDNE0250 

 HDNB0260 

 HDNE0270 

 HDNE02B0 

 HDNE0290 

 HDNB0300 

 aDNB0310 

 HDNB0320 

 MDNB0330 

 HDNB0340 

 HDNR0350 

 aDNB0360 

 HDNS0370 

 HDNB0330 

 HDNE0390 

 HDNB0400 

 HDNB0410 

 nDNB0420 

 HDNE0430 

 {1DNR0440 

 flDNR0450 

 HDNR0460 

 BDNE0470 

 aDNB0480 

 HDNE0490 

 HDNR0530 

 {1DNB0510 

 ODNE0520 

 HDNE0530 

 HDKB0543 

 HDNB0550 

 HDNB0560 

 HDNE0570 

 HDMB0580 

 HDNR0590 

 aDNB0600 

 nDNE0610 

 HDNE0620 

 flDNB0630 

 HDNE0640 

 HDNS0650 

 HDNE0660 

 HDNE0670 

 {IDNE0680 



