%%
%定义 多项式插值的Newton形式的函数
function Lout = Newton(x,xn,yn)
n = length(xn);
f = zeros(n,n);
f(:,1) = yn;
for k = 1:n-1
for i = k: n-1
f(i+1,k+1) = (f(i+1,k)-f(i,k))/(xn(i+1)-xn(i+1-k));
end
end
Lout = yn(1);
for k= 2:n
lk = 1;
for i = 1:k-1
lk=lk*(x-xn(i));
end
Lout = Lout+f(k,k)*lk;
end
评论0
最新资源