erfc.cpp

代码

double erfc(double x)
{
	double temp;
    if (x < 0.0)
	{
		temp= 1.0 + gammp(0.5, x*x);
	}
    else
	{
        temp = gammq(0.5, x *x);
	}
	return temp;
}