功能: 计算拉格朗日函数在一点的值。
优点: 逐点进行运算,后面步骤用带前面的计算结果。
输入: 给定的n+1个点 x0,x1,... xn 及y0,y1,... yn 和待求点xt
输出: 过给定n+1个点的拉格朗日函数L(x)在xt点的函数值
算法描述:
设I(n1,n2,n3,...,nk)表示过x(n1),x(n2),...x(nk)的插值结果
从两个点开始每怎加一个点计算一行,依次作线性差值,得到下表
叠代公式:
代码:
/* --- Neville ---
* input: n+1个节点 x0, x