MATLAB求解微分方程

%% 求解微分方程
%龙格库塔法
% t0=0;tf=10;          %范围0-10
% y0=2;                   %初始值y0
% [t,y]=ode23('funt',[t0,tf],y0);  %2、3阶龙格库塔法
% y1=sqrt(t+1)+1;    %原方程
% plot(t,y,'b.',t,y1,'r-')
%% 求解高阶微分方程
% t0=0;tf=20;                                  %范围0-20
% [t,y]=ode23('sys',[t0,tf],[1,0])        %3、4阶龙格库塔法(核心将高阶转化为一阶，即令y1=y',y2=y,有y2'=y1)
% subplot(1,2,1);plot(t,y(:,2));
% subplot(1,2,2);plot(y(:,2),y(:,1));

作者：TianSL