diff options
Diffstat (limited to 'CarpetExtra/IDHydroToy/src/gammp.f')
-rw-r--r-- | CarpetExtra/IDHydroToy/src/gammp.f | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/CarpetExtra/IDHydroToy/src/gammp.f b/CarpetExtra/IDHydroToy/src/gammp.f new file mode 100644 index 000000000..99a938857 --- /dev/null +++ b/CarpetExtra/IDHydroToy/src/gammp.f @@ -0,0 +1,16 @@ + FUNCTION gammp(a,x) + implicit none + REAL*8 a,gammp,x +CU USES gcf,gser + REAL*8 gammcf,gamser,gln + if(x.lt.0.or.a.le.0)pause 'bad arguments in gammp' + if(x.lt.a+1)then + call gser(gamser,a,x,gln) + gammp=gamser + else + call gcf(gammcf,a,x,gln) + gammp=1-gammcf + endif + return + END +C (C) Copr. 1986-92 Numerical Recipes Software t4. |