fdiffint.m

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

   function F=fdiffint(A,B,t,gam,a,b)
   a0=A(1); A=A(2:end); n=length(B); L=(b-a)/2;
   if gam>=0, F=0; else, F=a0*t^gam/gam; end
   for i=1:n
       an=i*pi/L; bn=gam*pi/2;
       F=F+an^gam*(A(i)*cos(an*t+bn)+B(i)*sin(an*t+bn));
   end
   if a+b, F=subs(F,t,t-L-a); end