(exp(-0.5*pi*b)cos(pi*a/2)-a)^2+(exp(-0.5*pi*b)sin(pi*a/2)-a)^2