用solve解含小量的非线性方程组出现不同结果

guohf

要解方程组

(x-6e-5)^2+(y-5e-5)^2=(8e-6)^2
x^2+(y-5e-5)^2=3.6e-9

方程有实数解，可采用两种方法，得不到想要的结果，且会报错，请高手给看下，先谢了

法一：

>> ax=6e-5;ay=5e-5;

>> syms x y;

>>[x,y]=solve((ax-x)^2+(ay-y)^2-(8e-6)^2==0,x^2+(y-5e-6)^2-3.6e-9==0)

x, y得到很怪的复数解。

将ax, ay的具体值代入，法二：

[x,y]=solve((x-6e-5)^2+(y-5e-5)^2-(8e-6)^2==0,x^2+(y-5e-6)^2-3.6e-9==0)

Error using mupadmex

Error in MuPAD command: Not a squarematrix. [(Dom::Matrix(Dom::ExpressionField()))::_power]

Error in sym/privBinaryOp (line 1693)

           Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});

Error in sym/mpower (line 189)

       B = privBinaryOp(A, p, 'symobj::mpower');

出现上面错误，不太理解这些“错误”指的是什么？在网上查了一下，添加了点：

[x,y]=solve((x-6e-5).^2+(y-5e-5).^2-(8e-6)^2==0,x.^2+(y-5e-6).^2-3.6e-9==0)

Warning: Explicitsolution could not be found.

> In solve at 179

x =  [ empty sym ]

y = []

guohf

是我写错了，把第二个方程中的5e-5写成5e-6了，真马虎。
其实我要做一个循环，把这次得到的x, y值作为新的ax, ay值，再进行求解，第一种方法较实用一点了