函數s =牛頓(x，y，x0，nn)
%牛頓插值，其中x是插值點，y是其函數值。
%x0是待插值點的x坐標值，nn是牛頓插值多項式的階。
nx =長度(x)；
ny =長度(y)；
如果nx~=ny
警告(“向量x和y應該有相同的長度”)
返回
結束
m =長度(x0)；
%根據公式計算所需插值點x0的每個元素。
因為i=1:m
t = 0.0
j = 1；
YY = y；
kk = j；
而(kk & lt=nn)
kk = kk+1；
對於k=kk:nx
YY(k)=(YY(k)-YY(kk-1))/(x(k)-x(kk-1))；
結束
結束
t = YY(1)；
對於k=2:nn
u = 1.0；
jj = 1；
while(jj & lt；k)
u = u *(x0(I)-x(jj))；
jj = jj+1；
結束
t = t+YY(k)* u；
結束
s(I)= t；
結束