fn rgamma_core(h: f64) -> f64
1/Γ(3/2 + h) for |h| <= 1/2 by Estrin’s scheme (short dependency chains).
1/Γ(3/2 + h)
|h| <= 1/2