34 SAVE erfcs, sqrtpi, nterf, xbig, sqeps, first
35 DATA erfcs( 1) / -.0490461212 34691808e0 /
36 DATA erfcs( 2) / -.1422612051 0371364e0 /
37 DATA erfcs( 3) / .0100355821 87599796e0 /
38 DATA erfcs( 4) / -.0005768764 69976748e0 /
39 DATA erfcs( 5) / .0000274199 31252196e0 /
40 DATA erfcs( 6) / -.0000011043 17550734e0 /
41 DATA erfcs( 7) / .0000000384 88755420e0 /
42 DATA erfcs( 8) / -.0000000011 80858253e0 /
43 DATA erfcs( 9) / .0000000000 32334215e0 /
44 DATA erfcs(10) / -.0000000000 00799101e0 /
45 DATA erfcs(11) / .0000000000 00017990e0 /
46 DATA erfcs(12) / -.0000000000 00000371e0 /
47 DATA erfcs(13) / .0000000000 00000007e0 /
48 DATA sqrtpi /1.772453850 9055160e0/
52 nterf =
inits(erfcs, 13, 0.1*r1mach(3))
53 xbig =
sqrt(-
log(sqrtpi*r1mach(3)))
54 sqeps =
sqrt(2.0*r1mach(3))
63 IF (y.LE.sqeps)
erf = 2.0*x/sqrtpi
64 IF (y.GT.sqeps)
erf = x*(1.0 +
csevl(2.*x**2-1., erfcs, nterf))
69 20
IF (y.LE.xbig)
erf = sign(1.0-
erfc(y), x)
70 IF (y.GT.xbig)
erf = sign(1.0, x)
function inits(OS, NOS, ETA)
std::string dimension(void) const
octave_value log(void) const
octave_value sqrt(void) const