Mathematica求解二元四次多项式根时Roots、FindRoot报错如何解决

报错核心原因 #

• 语法层面：Mathematica的表达式调用统一使用方括号[]，你写的f(x1,x2)会被系统判定为f乘以(x1,x2)的非法运算，不符合语法规范
• 命令适配层面：Roots仅支持单变量多项式求根，无法处理二元多项式；FindRoot为数值求根工具，必须为每个变量指定初始猜测值，仅填写{x1,x2}缺少必要的初始参数，所以触发报错

正确求解方案 #

首先明确：单个二元四次多项式f[x1,x2]==0的解是二维平面上的连续四次曲线，存在无穷多组解，没有有限个孤立根。如果需要得到有限个解，需要联立另一个二元方程组成方程组求解。

场景1：求解联立方程组的解析/数值解 #

如果有两个二元多项式方程f[x1,x2]==0和g[x1,x2]==0，可以用以下命令：

• 求精确解析解：
  Solve[f[x1, x2] == 0 && g[x1, x2] == 0, {x1, x2}]
• 求数值解（无需化简解析表达式时使用）：
  NSolve[f[x1, x2] == 0 && g[x1, x2] == 0, {x1, x2}]

场景2：求指定位置附近的单组数值解 #

如果只需要找某个初始猜测附近的一组解，使用FindRoot时给每个变量添加初始值即可，比如猜测根在x1=0、x2=1附近，写法为：
FindRoot[f[x1, x2] == 0, {x1, 0}, {x2, 1}]

附加：零点分布可视化 #

如果需要直观查看f[x1,x2]==0的零点曲线，可使用等高线绘图命令：
ContourPlot[f[x1, x2] == 0, {x1, -5, 5}, {x2, -5, 5}]
可根据需要调整x1、x2的取值范围。

内容的提问来源于stack exchange，提问作者Pramod Sharma