Newtons_secantmethod <- function(fun,x1,x2,ep=1e-5,max=100){
i<-0 ;k <- 1
while(k<=max){y1<-fun(x1,x2)$a; y2<-fun(x1,x2)$b;
x1 <- x2-y2*((x2-x1)/(y2-y1))
y1 <-fun(x1,x2)$a
x2 <- X1-y1*((x1-x2)/(y1-y2))
e <- sqrt((x1-x2)%*%(x1-x2))
if(e<ep){
i<- 1 ;break
}
k=k+1
}
list(root=x2,num=k,i=i)
}
fun<-function(x1,x2){
a<-x1^3-x1-1
b<-x2^3-x2-1
list(a=a,b=b)
}
Newtons_secantmethod(fun,1,2)
代码如上了,显示找不到对象,求大佬帮忙解决。