gammq.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 21 行

double gammq(double &a, double &x)
{
	double temp,gamser,gammcf,gln;
    if (x < 0.0 || a <= 0.0)
	{
        cout<< "pause";
        exit(1);
    }
    if (x < a + 1.0)
	{
        gser(gamser, a, x, gln);
        temp= 1.0 - gamser;
	}
    else
	{
        gcf(gammcf, a, x, gln);
        temp = gammcf;
	}
	return temp;   
}