R语言二分法编程
发布网友
发布时间：2022-05-21 14:44
共1个回答
热心网友
时间：2023-08-27 23:19
bisect(function(x) x^3-x-1,1,2)
结果：1.324718追问你这个肯定对，不过我就是做个练习，自己编，能告诉我我编的程序为什么运行不出正确结果吗？谢谢
追答你的程序输入了汉语状态下的“（”
fzero0)
list(fail="finding root is fail!")
else {
repeat{
if(abs(b-a)<eps)
break
x<-(a+b)/2
if(f(a)*f(x)<0)
b<-x
else
a<-x
}
list(root=(a+b)/2,fun=f(x))
}
}
f<-function(x) x^3-x-1
fzero(f,1,2,1e-6)
结果：1.324718和bisect函数得到的一致。
