| | Used this? Share your thoughts... Please login to comment. |
|---|
10 Mar 2008 | 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); |
akela_p501usa
