>> syms x y
>> s=solve('x+y=20+15*i','x*y=10+2*i')
s =
x: [2x1 sym]
y: [2x1 sym]
>> s.x,s.y
ans =
(592*i + 135)^(1/2)/2 + 10 + (15*i)/2
10 - (592*i + 135)^(1/2)/2 + (15*i)/2
ans =
10 - (592*i + 135)^(1/2)/2 + (15*i)/2
(592*i + 135)^(1/2)/2 + 10 + (15*i)/2