请教一个解方程组的问题

prince12th

用solve语句解这样一个方程组
* R" e0 ?5 E; E5 U! D(x1 - a1)^2 == b^2*(3*x1^6 + 6*x1^2*x2^4)
(x2 - a2)^2 == b^2*(3*x2^6 + 6*x2^2*x1^4)
% e) a. J. P# h/ @) f( v. ra1,a2已知，b是0.00171053
8 m: u! w6 a! E5 Z1 b; J* w
当设a1,a2都为1时，结果出来的是
$ U6 e- J1 g0 E6 U{{{x1 -> 1373.46, x2 -> 13.4299},
0 p. J- J. i' X5 k{x1 -> 10987.7,    x2 -> 15.3663},
  l. l& Y& p) Y/ a6 O* c) z+ g{x1 -> 1.05482*10^6,   x2 -> 17.8499},
- Z/ Q5 T9 |+ j/ Z/ c- X{x1 -> 2.81284*10^6, x2 -> 17.8426}}}
9 S% e% Y! T. B" q0 a
% B( L- y0 Q, g% ix1,x2不是应该对称的么，而且把结果代回原来的式子里也不正确，为什么呢？

sycxyc

有小数点的计算精度是有限的，并且高次方程没有精确解，如果一开始b的值就是以小数的形式，则中间的计算过程误差会越来越大，你需要将0.00171053化成分数的形式
1 A5 [* g5 z& ?/ V; x% F0 T' P. N1 YRationalize[0.00171053, 10^-16]