import sympy as sp
from IPython.display import display
z1,z2 = sp.symbols("z1 z2",complex=True)
x ,y  = sp.symbols("x y",real=True)
z1 = 3*x+(2*x-3*y)*1j
z2 = 6+1j
eq1 = sp.Eq(sp.re(z1),sp.re(z2))
eq2 = sp.Eq(sp.im(z1),sp.im(z2))
Soll=sp.solve([eq1,eq2],[x,y])
display(Soll)