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