fdiff3.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 9 行

   function Sd=fdiff3(S,gam)
   tt=S.breaks; t=tt(2:end-1); 
aa=[3:-1:0]'; aa1=gamma(aa'+1)./gamma(aa'+1-gam);
   for i=1:length(t)
       kk=find(tt<=t(i)); if length(kk)>0, k=kk(end); end
       yd(i)=(S.coefs(k,:).*aa1)*[[(t(i)-tt(k))*ones(3,1);t(i)].^[aa-gam]];
   end
   Sd=csapi(t,yd);