function dc = fun(t,c,L)%%%%%%L对应矩阵c的行数dc=zeros(L,1);dc1=zeros(L,1);N=L/2;% % NN=zeros(n,1);% % for i=1:n% % NN(i,1)=i-N-1;% % end%%%%%%%%%%%%测试L=4,N=2,时的特例。% % % % % dc(1) = -12*pi*1i*c(3)*c(1)-64*1i*pi^3*c(1);% % % % % dc(2) = 6*1i*pi*(c(3)*c(2)-c(1)*c(4))-8*1i*pi^3*c(2);% % % % % dc(3) = 0;代做C/C++实验、C语言程序代做、代写留学生C编程、C程序代写留学生% % % % % dc(4) =6*1i*pi*c(3)*c(4)+8*pi^3*1i*c(4);for i=1:Lk=i-N-1;%%%%%注意c向量对应的系数需要重新计算dc(i)=dc(i)-8*1i*pi^3*k^3*c(i);%%%%%%计算求和式的部分for m0=1:Lm=m0-N-1;n=k-m;if (n>=(-N))(nn0=n+N+1;dc1(i)=dc1(i)+n*c(m0)*c(n0);endenddc(i)=dc(i)-6*1i*pi*dc1(i);end& 转自：http://ass.3daixie.com/2018061815675167.html

原著是一个有趣的人,若有侵权,请通知删除