r - 在 R 中求解超越方程

Question

R 中有求解超越方程的函数吗？

例如我想解下面的方程

x = 1/tan(x)

有什么建议吗？我知道解决方案有多个根，所以我也希望能够恢复给定时间间隔内的所有答案

Answer 1

我会绘制函数曲线并查看它以查看它的外观:

R > y = function(x) { x - 1/tan(x) }
R > curve(y, xlim = c(-10, 10))
R > abline(h = 0, color = 'red')

然后我看到0和3之间有一个根，我会用uniroot来得到我想要的根:

R > uniroot(y, interval = c(0, 3))
$root
[1] 0.8603

$f.root
[1] 6.612e-06

$iter
[1] 7

$estim.prec
[1] 6.104e-05