fsolve accuracy with complex variables [SOLVED]

[JMB]

Thank you very much for your answers.

My Epsilon is set to 1e-12. Changing it to different values made no difference.

The solution provided by Parisse works just fine. After the definitions:

z:=x+i*y
eq:=|z|^2-200*z+500+500*i

Both solve and fsolve return the correct solution:

solve([re(eq),im(eq)],[x,y],[190,1]) -> [197.43587635 2.5]

fsolve([re(eq),im(eq)],[x,y],[190,1]) -> [197.43587635 2.5]