r语言程序题

R语言学习的人不多。第二题： td<-function(fx,a,b,r){ #(fx:是表达式，方程是fx=0,变量是x；a，b，区间[a,b],r是精度。 c<-(a+b)/2x<-cif(abs(eval(fx))

R语言学习的人不多。

第二题：

td<-function(fx,a,b,r){ #(fx:是表达式，方程是fx=0,变量是x；a，b，区间[a,b],r是精度。

c<-(a+b)/2

x<-c

if(abs(eval(fx))<r) x else

{

x<-ada<-eval(fx)

x<-bdb<-eval(fx)

x<-cdc<-eval(fx)

if(da*db<0) td(fx,a,c,r) else td(fx,c,d,r)

}

}

调用函数：

fx<-expression(x*x-x-1)

a<-1b<-2r<-1e-5

td(fx,a,b,r)

---