thp.m

来自「牛顿插值法」· M 代码 · 共 19 行

function r=thp(x,y)
 m=length(x);
 h=0.1;
  for k=1:m
      s1=spline(x,y,x(k)+2*h);
      s2=spline(x,y,x(k)+h);
      s3=spline(x,y,x(k)-h);
      s4=spline(x,y,x(k)-2*h);
      if k==1
           r(k)=(-3*y(1)+4*s2-s1)/(2*h);
      else if  k==m
           r(k)=(s4-4*s3+3*y(m))/(2*h);
              else
           r(k)=(s2-s3)/(2*h);
       end
   end
end