Evaluate the polynomial interpolant by Neville’s method in Python,,''' p = nevi


''' 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]

评论关闭