fn rgamma(z: f64) -> f64
1/Γ(z) for moderate real z: DLMF 5.5.1 shifts into [1, 2), then a Taylor polynomial. The shifts multiply by z + k exactly where Γ has poles, so those zeros are exact.
1/Γ(z)
z
[1, 2)
z + k
Γ