fdiff4.m

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

   function [t,yd]=fdiff3(S,gam)
   tt=S.breaks; M=[ones(3,1); 0]; 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,:)*[[(t(i)-tt(k))*ones(3,1);t(i)].^aa];
       %       yd(i)=(S.coefs(k,:).*aa1)*[(t(i)-tt(k))*M.^(aa-gam)];
       yd(i)=(S.coefs(k,:).*aa1)*[[(t(i)-tt(k))*ones(3,1); t(i)].^[aa-gam]];
   end