hermite.asv

来自「一阶导数hermite插值。」· ASV 代码 · 共 19 行

function y=hermite(x0,y0,y1,x)
%hermite insert ★x0为输入节点的向量；y0为y的值向量，y1为
n=length(x0);
m=length(x);
for k=1:m
    yy=0.0;
    for i=1:n
        h=1.0;
        a=0.0;
        for j=1:n
            if j~=i
                h=h*((x(k)-x0(j))/(x0(i)-x0(j)))^2;
                a=1/(x0(i)-x0(j))+a;
            end
        end
        yy=yy+h*((x0(i)-x(k))*(2*a*y0(i)-y1(i))+y0(i));
    end
    y(k)=yy;
end