网友采纳 哦谢谢啦那最后检验的第三个解为什么会出现一个False
汪柏林的回答:
网友采纳 ……好吧,我都没注意到。我运算精度没取够:eqn={x^2-y==1,x+y^2==0}sol=Solve[eqn,{x,y}];eqn/.N[sol,16]取到16位精度就可以了。其实这里最正确的做法是在最后一行使用FullSimplify[eqn/.sol]但是这个的结果根式比较复杂,这样子计算速度太慢了,我不愿等……等下,我想到个满足全部要求的方法了。那就是先避免对原方程组(本质上是个四次方程组)的显式求解(不使用根式表示解,而使用方程表示解),然后单独化出一个显式根式解,再对非显式解去用FullSimplify,这样化简会很快:eqn={x^2-y==1,x+y^2==0}sol=Solve[eqn,{x,y},Quartics->False]ToRadicals@solFullSimplify[eqn/.sol]