NSolve无法求解隐式微分：x=1处斜率求解的替代方案咨询

你遇到的核心问题是NSolve不适合求解超越方程——你的方程x Sin[y] - y Sin[x] == 2在x=1时转化为Sin[y] - y Sin[1] == 2，这是一个包含三角函数的超越方程，没有解析解，而NSolve主要用于处理多项式类方程，所以无法给出结果。我们可以用数值求解工具FindRoot来解决这个问题，步骤如下：

1. 回顾已完成的隐式微分 #

你已经正确通过隐式微分求出了导数表达式：

eqn[x_, y_] := x*Sin[y] - y*Sin[x] == 2
yPrime = Solve[D[eqn[x, y[x]], x], y'[x]] /. {y[x] -> y, y'[x] -> y'} // Simplify

得到的导数为：

{{y' -> (y Cos[x] - Sin[y])/(x Cos[y] - Sin[x])}}

2. 用FindRoot求解x=1时的y值 #

因为方程是超越方程，我们需要在给定区间-5<=y<=5内用数值迭代法找根。首先可以先可视化函数在y∈[-5,5]的取值，确认根的大致位置：

Plot[Evaluate[eqn[1, y] /. Equal -> Subtract], {y, -5, 5}, 
     AxesLabel -> {"y", "f(y)"}, GridLines -> Automatic, 
     PlotLabel -> "f(y) = Sin[y] - y Sin[1] - 2"]

从图中可以看到，函数在y≈-2和y≈-4附近穿过0，这就是我们要找的根。接下来用FindRoot指定初始值求解：

(* 求解第一个根，初始值设为-2 *)
sol1 = FindRoot[eqn[1, y], {y, -2}]
(* 求解第二个根，初始值设为-4 *)
sol2 = FindRoot[eqn[1, y], {y, -4}]

运行后会得到两个数值解：

• {y -> -1.97573}
• {y -> -3.94151}

3. 代入导数计算x=1处的斜率 #

把得到的y值代入之前的导数表达式，计算对应斜率：

(* 计算第一个解对应的斜率 *)
slope1 = y' /. yPrime[[1]] /. x -> 1 /. sol1
(* 计算第二个解对应的斜率 *)
slope2 = y' /. yPrime[[1]] /. x -> 1 /. sol2

最终会得到两个斜率值：

• 约-0.2765
• 约0.7045

为什么NSolve不行？ #

简单来说，NSolve的设计目标是求解多项式方程或可转化为多项式系统的方程，它无法处理包含三角函数、指数函数等的超越方程。而FindRoot是基于牛顿迭代法的数值求解工具，只要给定合适的初始值，就能找到方程在该初始值附近的根，非常适合这类超越方程的求解。

内容的提问来源于stack exchange，提问作者Derek Hassan