标准正态分布函数的近似计算

之前想写个程序自动分析数据的分布，但卡在无法求正态分布的分布函数了，无意中复习概率论课程，发现在附录中居然有近似的计算公式！太高兴了记录下来

#define pi (3.1415926535898)
#define a0  (0.33267)
#define a1  (0.4361836)
#define a2 (-0.1201676)
#define a3 (0.9372986)
double test(double x)
{
double t = (1.0/(1.0+a0*fabs((x))));
return exp(-x*x/2.0) * t * (a1 + a2*t + a3*t*t) / sqrt(2.0*pi);
}
double testfun(double x)
{
if(x == 0)
{
return 0.5;
}
else if(x > 0)
{
return 1.0 - test(x)/sqrt(2*pi);
}
else
{
return test(x)/sqrt(2*pi);
}
}