matlab Newton向前插值

**********************本文章由计算机专业的孙旭编写，复制请谨慎

function Fun()
clc;
X=[0 1 2 3];
Y=[1 2 17 64];
x0=0.5;
h=1;
New_Front(X,Y,x0);
end


%子函数
function New_Front(X,Y,x0,h)
len_x=length(X);
len_y=length(Y);
if (len_x~=len_y) 
    disp('你输入的维度不相等');
    return;
else
    syms s;
    p=Y(1);
    q=1;
    f=1;
    for i=1:1:len_x-1
        Y(i+1:len_x)=Y(i+1:len_x)-Y(i:len_x-1);
        A(i)=Y(i+1);
        q=q*(s-i+1);
        f=f*i;
        p=p+(q/f)*A(i);
    end
    
    %化简多项式
    simplify(p);
    
    if nargin==4
        disp('所求的插值函数的值为');
        p0=subs(p,'s',(x0-X(1)/h));
        vpa(p0,6),
    else 
        disp('所求牛顿的向前插值多项式为: ');
        p=collect(p);
        p=vpa(p,6),
    end
end
end