pub fn intkambe(n: i32, z: Complex, eta: Complex) -> Result<Complex>Expand description
The Kambe integral I_n(z, eta) = integral_eta^infinity t^n exp(-z^2 t^2/2 + 1/(2t^2)) dt.
Upstream: treams.special.intkambe. Differences: |n| above 260 and non-finite
arguments give an error. The odd orders -3 and -1 come from a series over scaled
incomplete gamma functions, and lower orders from that series or the downward
recurrence of treams, whichever bounds the rounding error more tightly.