绘制内维尔的方法

Question

我在R中写了这个Neville的方法，但我也熟悉Matlab。所以你在R或Matlab中的任何评论对我来说真的很棒。

基于这些数据点：

x= 1.5
xi = c(1,1.3,1.6,1.9,2.2)
fi = c(0.7651977, 0.6200860, 0.4554022, 0.2818186, 0.1103623)

以下是我在R中的脚本：

neville <- function(x,xi,fi){
  n <- length(xi) - 1
  Q <- matrix(0, n+1, n+1)
  Q[,1] = fi
  for (i in 1:n) {
    for (j in 1:i){
       Q[i+1,j+1] <- ((x - xi[i-j +1])*Q[i+1,j] - (x - xi[i+1])
                 *Q[i,j])/(xi[i+1] - xi[i-j+1])
    }
  }
  return (Q)
}

显然，我的输出是矩阵形式。那么，任何人都可以告诉我如何绘制这个Neville方法的插值来表示我的xi和fi点吗？