Matlab Nevilles插值

@Matlab

Neville’s Iterated Interpolation

In mathematics, Neville’s algorithm is a method of computing interpolating polynomials, developed by mathematician Eric Harold Neville.

Algorithm step

To evaluate the interpolating polynomial P on the n + 1 distinct numbers x₀,…,x_n at the number x for the function f:

Input
numbers x,x₀,x₁,…x_n;
values f(x₀),f(x₁),…,f(x_n) as the first column Q_0,0,Q_1,0,…,Q_n,0 of Q.
Output
the table Q with P(x)=Q_n,n.

Step prompt
Step 1
for i = 1,2,…,i
for j = 1,2,…,j
set Q_i,j=((x-x_i-j)Q_i,j-1-(x-x_i)Q_i-1,j-1)/(x_i-x_i-j).

Step 2
Output(Q)
Stop

The code is as follows

function feval = NevilleInterpolation(xs, fs, xeval)
n