function y=hermite(x0,y0,y1,x)
%hermite insert ★x0為輸入節(jié)點的向量；y0為y的值向量，y1為相應(yīng)節(jié)點一階倒數(shù)的函數(shù)值的向量，x為所要求的插值節(jié)點。
%只適用一階倒數(shù)。
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