blob: 99a938857f6dcc97a169305eb713b755f42b22e3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
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.
|