newton.m

来自「基于matlab实现的实质分析实例的计算 对基本的插值的计算 可运行」· M 代码 · 共 20 行

%牛顿插值算法程序(针对电子教案中例题)
x = [0	1	2	3	4];
y = [1	4	15	40	85];
v = 3.5;
n = 4;
N = y(1);
t = 1;
for k = 2:(n+1)
    z = y;
    for i = 0:(n+1-k)
        
        y(k+i) = (z(k+i)-z(k+i-1))/(x(k+i)-x(k+i-(k-1)));
    end
    t = t*(v-x(k-1));
    N = N+y(k)*t;
end
N
N1 = v^3+v^2+v+1