R语言实现牛顿前插与后插

Question

#这是只编了前插的一部分

N.Diff<-function(x,f,X){
  k<-length(x)           #定义节点个数
  answer<-0              #定义结果初始值
  CFt<-matrix(NA,k,k)  #定义差分表
  if(abs(X-x[1])<=abs(X-x[k])){      #向前差分
    t<-(X-x[1])/(x[2]-x[1])      #定义间距
    CFt[i,1]<-f[i]    #差分表第一列为f
    for(j in 2:k){
      for(i in 1:k){
       if (i+j>(k+1)) break
        CFt[i,j]<-CFt[i+1,j-1]-CFt[i,j-1]    #差分表的元素
        }
      }
    lc<-rep(1,k);cf<-rep(0,k)
    for(j in 2:k) {lc[j]<-prod(t-0:(j-2))/prod(1:(j-1))}    #定义连乘式
    for(j in 1:k) {cf[j]<-CFt[1,j]} 
    Nn<-cf%*%lc
    list(CFTable=CFt,Nnx=Nn)}}

 

x<-seq(0,0.5,by=0.1)
f<-c(1,0.995,0.98007,0.95534,0.92106,0.87758)
X<-0.048
N.Diff(x,f,X)

 

但是显示找不到对象，想问问是为什么，具体哪里出了问题

Answer 1

不知道你这个问题是否已经解决, 如果还没有解决的话:

• 可以查看手册：r语言 获得R 中的内容

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 以帮助更多的人 ^-^