解决多项式分数微分方程的一阶隐乘积梯形法则求解 - Matlab开发

在区间[t0,T]上求解多项式分数阶微分方程的初值问题lam_Q D^(al_Q) y(t) + ... + lam_1 D^(al_1) y(t) = f(t,y(t)) y(0) = y0(1), y'(0) = y0(2), ... y^m(0) = y0(m)，其中m是大于max(al_1,...,al_Q)的最小整数。该问题通过收敛阶数为1的矩形类型的隐式积分规则解决。更多信息请参见以下论文[1] Garrappa R.：分数阶微分方程的数值解：调查和软件教程，数学2018, 6(2), 16 doi: https://doi.org/10.3390/math6020016可下载的pdf： http : //www.mdpi.com/2227-7390/6/2/16/pdf