if b**2-4*a*c<0: r=(-b)/(2*a) u=((4*a*c-b*b)**(1/2))/(2*a) x1=complex(r,u) x2=complex(r,-u) print(x1,x2)