为什么不推荐用R语言uniroot解一元二次方程，报错“位于极点边的f()值之正负号不相反”的原因及解决办法

问题描述

如果需要同时解多个一元二次方程，如果在网上查的话99%的人会推荐使用R语言的unitroot函数，但实际上这个函数非常不适合满足这类需求，很容易报如下错误：

Error in uniroot(f2, c(-10, 0), a = a, b = b, c = c, tol = 1e-04) :
位于极点边的f()值之正负号不相反

原因

首先可以先看看uniroot是怎么解一元二次方程的：

一元二次方程：ax^2+bx+c=0，设a=1，b=5，c=6，求x？

f2 <- function (x, a, b, c) ax^2+bx+c
a<-1;b<-5;c<-6
result <- uniroot(f2,c(0,-2),a=a,b=b,c=c,tol=0.0001)
result$root
[1] -2

把参数带入方程，用uniroot()函数，我们就解出了方程的一个根，改变计算的区间，我们就可以得到另一个根。

result <- uniroot(f2,c(-4,-3),a=a,b=b,c=c,tol=0.0001)
result$root
[1] -3

方程的两个根，一个是-2，一个是-3。

重点在这&#x