主函数比较
% 根据文献自己编写
q = 0.5; %分数阶阶数
fdefun = @(t,y) [2/gamma(3-q)*t^(2-q)-1/gamma(2-q)*t^(1-q)-y+t^2-t]; %一元微分方程
y0 =0;%初值y0为列向量
h = 2^(-6);%步长
tspan = [0,2];
[t,y] = fdewfc(q,fdefun,tspan,y0,h);
figure(1)
plot(t,y(1,1:end)) ;
xlabel('t'); ylabel('y(t)');
hold on
plot(t,t.^2-t+0.5,'r-.')%解析解
%plot(t,-2./(t+1),'r')
%%%%%%与意大利 Roberto Garrappa, University of Bari, Italy 结果比较
t0=tspan(1); tfinal = tspan(2);
[t, y_fde12] = fde12(q,fdefun,t0,tfinal,y0,h);
plot(t,y_fde12+1,'k.')
Jfdefun = @(t,y)[-1];%Jacobi阵 此处退化为导数
[t, y_flmm2] =