Evaluate the polynomial interpolant by Neville’s method in Python
''' p = neville(xData,yData,x).
Evaluates the polynomial interpolant p(x) that passes
trough the specified data points by Neville's method.
'''
def neville(xData,yData,x):
m = len(xData) # number of data points
y = yData.copy()
for k in range(1,m):
y[0:m-k] = ((x - xData[k:m])*y[0:m-k] + \
(xData[0:m-k] - x)*y[1:m-k+1])/ \
(xData[0:m-k] - xData[k:m])
return y[0]
开心洋葱 , 版权所有丨如未注明 , 均为原创丨未经授权请勿修改 , 转载请注明Evaluate the polynomial interpolant by Neville’s method in Python！
