newtint.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 18 行

clear
clc
%function newtint(x,y,xhat)
x=[0.4 0.55 0.65 0.8 0.9];
y=[0.41075 0.57815 0.69675 0.88811 1.02653];
xhat=0.895;%0.596;
n=length(y);
c=y(:);
for j=2:n
    for i=n:-1:j
        c(i)=(c(i)-c(i-1))/(x(i)-x(i-j+1));
    end
end
yhat=c(n);
for i=n-1:-1:1
    yhat=yhat*(xhat-x(i))+c(i);
end
fprintf('N(%0.3f)=%f\n',xhat,yhat);