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;   
}