Here is code written in java for IncGamma.This code uses my gamma stirling method if are you interested i can post the code.
[bash]public double IncGammaLow(double x,double a){
double sum = 0;
double result = 0;
if(x < 0.0f){
throw new IllegalArgumentException("arg x must be greater than 0 t" + x);
}else if( a < 0.0f || a == 0.0f || a < 1.0f){
throw new IllegalArgumentException("arg a must be not less than 1.0 t" + a);
}else{
double exp,gamma_arg,gamma,nom,denom,power,gamma2;
long k;
int i,j;
double[]array = new double[32];
double run_value;
k = -1;
run_value = 0;
for(i = 0;i
= run_value;
}
for(j = 0;j;
}
exp = Math.exp(-x);
power = Math.pow(x, a);
gamma2 = this.GammaStirlingApprox(a);
result = power*gamma2*exp*sum;
}
return result;
}[/bash]
