138 static double P[] = {
139 2.46196981473530512524E-10,
140 5.64189564831068821977E-1,
141 7.46321056442269912687E0,
142 4.86371970985681366614E1,
143 1.96520832956077098242E2,
144 5.26445194995477358631E2,
145 9.34528527171957607540E2,
146 1.02755188689515710272E3,
147 5.57535335369399327526E2
150 static double Q[] = {
152 1.32281951154744992508E1,
153 8.67072140885989742329E1,
154 3.54937778887819891062E2,
155 9.75708501743205489753E2,
156 1.82390916687909736289E3,
157 2.24633760818710981792E3,
158 1.65666309194161350182E3,
159 5.57535340817727675546E2
162 static double R[] = {
163 5.64189583547755073984E-1,
164 1.27536670759978104416E0,
165 5.01905042251180477414E0,
166 6.16021097993053585195E0,
167 7.40974269950448939160E0,
168 2.97886665372100240670E0
171 static double S[] = {
173 2.26052863220117276590E0,
174 9.39603524938001434673E0,
175 1.20489539808096656605E1,
176 1.70814450747565897222E1,
177 9.60896809063285878198E0,
178 3.36907645100081516050E0
181 static double T[] = {
182 9.60497373987051638749E0,
183 9.00260197203842689217E1,
184 2.23200534594684319226E3,
185 7.00332514112805075473E3,
186 5.55923013010394962768E4
189 static double U[] = {
191 3.35617141647503099647E1,
192 5.21357949780152679795E2,
193 4.59432382970980127987E3,
194 2.26290000613890934246E4,
195 4.92673942608635921086E4
198 #define UTHRESH 37.519379347
214 y = 0.5 + 0.5 *
erf(
x);
229 double p,
q,
x,
y,
z;
299 return (1.0 -
erfc(
x));