>> [x,y]=solve('x^2+(y-100)^2-89^2','x^2+y^2-10000');
>> x=double(x),y=double(y)
x =
79.7022
-79.7022
y =
60.3950
>> s=sqrt((x-200).^2+y.^2-134^2)
s =
12.7715
252.8337
严格来说没有解，近似解为（79.7022， 60.3950）
