Quadratic substitution for a square-root endpoint singularity (source code)

= Quadratic substitution for a square-root endpoint singularity

If $f(x)=x^{-1/2}g(x)$ with $g$ smooth near zero, the substitution $x=s^2$ gives
$$
\int_0^1f(x)\,dx=\int_0^12g(s^2)\,ds.
$$
The transformed integrand is smooth, so an ordinary second-order midpoint rule recovers its usual convergence rate.