glfdiffs.m

来自「高等引用数学的matlab求解的配套程序」· M 代码 · 共 8 行

   function dy=glfdiff(y,t,gam)
   h=t(2)-t(1); dy(1)=0; y=y(:); t=t(:);
   for i=2:length(t)
      L=min([i,300]);
       w=1; for j=2:L, w(j)=w(j-1)*(1-(gam+1)/(j-1)); end
      dy(i)=w*[y(L:-1:1)]/h^gam;
   end