solve ((x1-m)^2/a^2+(y1-n)^2/b^2=1, (x2-m)^2/a^2+(y2-n)^2/b^2=1), (m, n)