myDataBus

Summary
Comments

Posted 10 Mar 2008 by

akela_p501usa
static public double igamc( double a, double x )
throws ArithmeticException {
double big = 4.503599627370496e15;
double biginv = 2.22044604925031308085e-16;
double ans, ax, c, yc, r, t, y, z;
double pk, pkm1, pkm2, qk, qkm1, qkm2;

if( x <= 0 || a <= 0 ) return 1.0;

if( x < 1.0 || x < a ) return 1.0 - igam(a,x);

ax = a * Math.log(x) - x - lgamma(a);
if( ax < -MAXLOG ) return 0.0;

ax = Math.exp(ax);

Feel free to add your comments.