from sympy import *
from IPython.display import display
z,z1,z2,z3=symbols('z z1:4',complex=True)
f=z1**2-(1+sqrt(3)*I)
z2=solve(f,z1)
for z3 in z2:
    z=re(z3)+im(z3)*I
    display(z)