25 static double p[7] = { 220.2068679123761, 221.2135961699311,
26 112.079291497870, 33.91286607838300, 6.37396220353165,
27 0.7003830644436881, 0.352624965998910e-1
29 static double q[8] = { 440.4137358247522, 793.8265125199484,
30 637.3336333788311, 296.5642487796737, 86.78073220294608,
31 16.06417757920695, 1.755667163182642, 0.8838834764831844e-1
33 static double cutoff = 7.071, root2pi = 2.506628274631001;
35 double pee, queue, pdf;
53 expntl = exp(-0.5 * zabs * zabs);
54 pdf = expntl / root2pi;
57 pee = expntl * ((((((p[6] * zabs + p[5]) * zabs + p[4])
58 * zabs + p[3]) * zabs + p[2]) * zabs +
60 / (((((((q[7] * zabs + q[6]) * zabs + q[5]) * zabs + q[4])
61 * zabs + q[3]) * zabs + q[2]) * zabs + q[1]) * zabs + q[0]);
63 pee = pdf / (zabs + 1.0 / (zabs + 2.0 / (zabs + 3.0 / (zabs + 4.0 /