部分代码:
%% 应用三级三阶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);
1