本文实例讲述了python计算牛顿迭代多项式的方法。分享给大家供大家参考。具体实现方法如下:
''' p = evalpoly(a,xdata,x). evaluates newton's polynomial p at x. the coefficient vector 'a' can be computed by the function 'coeffts'. a = coeffts(xdata,ydata). computes the coefficients of newton's polynomial.''' def evalpoly(a,xdata,x): n = len(xdata) - 1 # degree of polynomial p = a[n] for k in range(1,n+1): p = a[n-k] + (x -xdata[n-k])*p return pdef coeffts(xdata,ydata): m = len(xdata) # number of data points a = ydata.copy() for k in range(1,m): a[k:m] = (a[k:m] - a[k-1])/(xdata[k:m] - xdata[k-1]) return a
希望本文所述对大家的python程序设计有所帮助。