fn pi_cot_scaled(r: f64) -> f64
πr cot(πr) for |r| <= 1/2, accurate (and finite) as r → 0.
πr cot(πr)
|r| <= 1/2
r → 0