fn approx_parabola_integral(x: f64) -> f64
An approximation to $\int (1 + 4x^2) ^ -0.25 dx$
This is used for flattening curves.