[x1,y1,x2,y2]=solve('x1^2 + y1^2=1','x2^2-8*x2 +y2^2 +15=0','x1*x2 + y1 * y2=1','x1 + x2 =a','x1','y1','x2','y2') [x1,y1,x2,y2] = solve(... x1^2/r1^2 + y1^2/r2^2 == 1,... (x2-a)^2 + (y2-b)^2 == r3^2,... x1*x2/r1^2 + y1*y2/r2^2 == 1,... (x2-x1)^2 + (y