(abs(sin(-2*atan(sqrt(b/a))))/sqrt((b^2-a^2)*(cos(-2*atan(sqrt(b/a)))+1)/2+a^2))