>> [x,y]=solve('sqrt(x^2+y^2)-sqrt((x-100)^2+y^2)=a','sqrt((x-200)^2+y^2)-sqrt((x-100)^2+y^2)=b','x,y')
Warning: possibly spurious solutions [solvelib::checkVectorSolutions]
x =
(a^2*b - a*b^2 + 30000*a + 10000*b)/(200*a + 200*b)
(a^2*b - a*b^2 + 30000*a + 10000*b)/(200*a + 200*b)
y =
(- (a^4*b^2)/10000 + a^4 + (a^3*b^3)/5000 - 2*a^3*b - (a^2*b^4)/10000 + 6*a^2*b^2 - 50000*a^2 - 2*a*b^3 + 20000*a*b + b^4 - 50000*b^2 + 400000000)^(1/2)/(2*a + 2*b)
-(- (a^4*b^2)/10000 + a^4 + (a^3*b^3)/5000 - 2*a^3*b - (a^2*b^4)/10000 + 6*a^2*b^2 - 50000*a^2 - 2*a*b^3 + 20000*a*b + b^4 - 50000*b^2 + 400000000)^(1/2)/(2*a + 2*b)
Warning: possibly spurious solutions [solvelib::checkVectorSolutions]
x =
(a^2*b - a*b^2 + 30000*a + 10000*b)/(200*a + 200*b)
(a^2*b - a*b^2 + 30000*a + 10000*b)/(200*a + 200*b)
y =
(- (a^4*b^2)/10000 + a^4 + (a^3*b^3)/5000 - 2*a^3*b - (a^2*b^4)/10000 + 6*a^2*b^2 - 50000*a^2 - 2*a*b^3 + 20000*a*b + b^4 - 50000*b^2 + 400000000)^(1/2)/(2*a + 2*b)
-(- (a^4*b^2)/10000 + a^4 + (a^3*b^3)/5000 - 2*a^3*b - (a^2*b^4)/10000 + 6*a^2*b^2 - 50000*a^2 - 2*a*b^3 + 20000*a*b + b^4 - 50000*b^2 + 400000000)^(1/2)/(2*a + 2*b)