Trailing-Edge
-
PDP-10 Archives
-
decuslib10-01
-
43,50144/errlt.src
There are no other files named errlt.src in the archive.
00010 C ***PROGRAM TO EVALUATE THE ERROR FUNCTION AND ITS COMPLEME
00020 + NT
00030 FUNCTION ERR(W)
00040 DIMENSION A(25),B(30)
00050 M=24
00060 A(1)=16443152242714E-13
00070 A(2)=-9049760497548E-13
00080 A(3)= 0643570883797E-13
00090 A(4)= 0196418177368E-13
00100 A(5)=-0001244215694E-13
00110 A(6)=-0009101941905E-13
00120 A(7)=-0001796219835E-13
00130 A(8)= 0000139836786E-13
00140 A(9)= 164789417E-13
00150 A(10)= 39009267E-13
00160 A(11)=- 00893145E-13
00170 A(12)=- 03747896E-13
00180 A(13)=- 01298818E-13
00190 A(14)=-00136773E-13
00200 A(15)= 00077107E-13
00210 A(16)= 00046810E-13
00220 A(17)= 00011844E-13
00230 A(18)=- 0005E-13
00240 A(19)=- 1384E-13
00250 A(20)=- 0652E-13
00260 A(21)= 0145E-13
00270 A(22)= 0010E-13
00280 A(23)= 0024E-13
00290 A(24)= 0011E-13
00300 A(25)= 0002E-13
00310 X=ABS (W)
00320 IF (ABS (X)-.01) 1, 2, 2
00330 1 XERR=2.0/(3.0*1.77245385)*X*(3.0-X**2)
00340 GO TO 6
00350 2 Z=(X-1.0)/(X+1.0)
00360 DO 3 I=1,30
00370 B(I)=0.
00380 3 CONTINUE
00390 DO 4 I=1,M
00400 M1=(M+1)-I
00410 B(M1)=2.0*Z*B(M1+1)-B(M1+2)+A(M1+1)
00420 4 CONTINUE
00430 F=-B(2)+Z*B(1)+.5*A(1)
00440 XERR=1.0-(1.0/1.77245385)*(EXP (-(X**2)))*F
00450 IF (ABS (X)-.01) 6, 7, 7
00460 6 CERR=1.0-XERR
00470 GO TO 5
00480 7 CERR=(1.0/1.77245385)*(EXP (-(X**2)))*F
00490 5 IF (W) 9, 8, 8
00500 8 ERR=XERR
00510 GO TO 13
00520 9 ERR=CERR
00530 13 RETURN
00540 END