18 lines
292 B
C
18 lines
292 B
C
|
|
float gammp(a,x)
|
|
float a,x;
|
|
{
|
|
void gcf(),gser();
|
|
void nrerror();
|
|
float gamser,gammcf,gln;
|
|
|
|
if (x < 0.0 || a <= 0.0) nrerror("Invalid arguments in routine gammp");
|
|
if (x < (a+1.0)) {
|
|
gser(&gamser,a,x,&gln);
|
|
return gamser;
|
|
} else {
|
|
gcf(&gammcf,a,x,&gln);
|
|
return 1.0-gammcf;
|
|
}
|
|
}
|