zgamr.f
620 Bytes
complex(8) function zgamr (z)
implicit none
c july 1977 edition. w. fullerton, c3, los alamos scientific lab.
c this version is an inaccurate preliminary one. eventually this
c routine should be a fundamental routine with no dependence on cgamma.
c
complex(8) z, zlngam
real(8) x
integer irold,ir
external zlngam
c
zgamr = (0.0, 0.0)
x = real (z)
if (x.le.0.0 .and. aint(x).eq.x .and. 0.0.eq.aimag(z)) return
c
call entsrc (irold, 1)
zgamr = zlngam(z)
call erroff
call entsrc (ir, irold)
zgamr = exp (-zgamr)
c
return
end