fn approx_parabola_integral<S: Scalar>(x: S) -> S
Compute an approximation to integral (1 + 4x^2) ^ -0.25 dx used in the flattening code.