#include "incomplete-gamma.h"
#include <iostream>
using namespace std;

int main()
{
    double alpha=67.81076;
    double sigma=0.006545152;

    for (double x=0.3; x<2.0; x+=0.01)
    {    cout << x 
    		<< "\t" << GammaP(alpha,x/sigma) 
		<< "\t" << GammaQ(alpha,x/sigma) 
		<< "\n";
    }
}

// Change Log:
// 30 March 2004 Kevin Karplus
//	Modified to match current .h files.
//	Interpretation of test still not obvious to the uninitiated.