def sqrt_complex(c): r=sqrt(abs(c)) phi=phase(c)/2 x=r*cos(phi) y=r*sin(phi) return complex(-x,-y),complex(x,y)