fdiff3.m

高等引用数学的matlab求解的配套程序

   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);