inversehaversine(inversebetaregularized((x-sin(x)(-2/3 cos(x/2)^2-1))/pi,1/2,5/2))-x