glfode.m

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

   function y=glfode(u,t,a,n)
   h=t(2)-t(1); D=sum(a./[h.^n]); W=[]; y=zeros(length(t),1); w=y;
   for i=1:length(n)
      w=1; for j=2:length(t), w(j)=w(j-1)*(1-(n(i)+1)/(j-1)); end
      W=[W; w];
   end
   for i=2:length(t)
      for j=1:length(n), A(j)=W(j,2:i)*[y(i-1:-1:1)]; end
      y(i)=(u(i)-sum(A.*a./[h.^n]))/D;
   end