部分代码： %% 应用三级三阶Runge-Kutta 方法与复合Gregory求积公式到问题(1.2,1.3,1.4) %% 通用函数 [t,y]=main(ddefun,kernelfun,initialfun,lag,tspan,dimensional) %% 其中 ddefun为右端函数，kernelfun为积分核函数，initialfun为初始函数,lag为延迟量， %% tspan为求解区间，dimensional为问题维数 % %%%% 问题1.2 % [t,x]=main(@fun_f1,@fun_g1,@fun_varphi1,0.8,[0,8],1); %plot(t,x) %xlabel('\it{t}'); %ylabel('\it{x(t)}'); %%%% 问题1.3 %[t,x]=main(@fun_f2,@fun_g2,@fun_varphi2,0.8,[0,8],1); %plot(t,x) % xlabel('\it{t}'); % ylabel('\it{x(t)}'); %%% 问题1.4 [t,x]=main(@fun_f3,@fun_g3,@fun_varphi3,1/6,[0,5],2); x1 =x(:,1); =x(:,1);