%% 求解微分方程
%龙格库塔法
% 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