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