static double GAMMADIST(double x, double alpha, double beta, bool cumulative)
{
GammaDistribution dist = new GammaDistribution(alpha, beta);
return(cumulative
? dist.CumulativeDistribution(x)
: dist.ProbabilityDensity(x));
}