ode45、ode15s、ode15i怎么选？
%主程序-求解器：ode45clear allclose allclcy0 = 0;tspan = [0:0.01:10];[t,y] = ode45(@subfun,tspan,y0);figure(1)plot(t,y);% 画出每次积分结果与时间的图形%主程序-求解器：ode15sclear allclose allclcy0 = 0;yp0 = 1;tspan = [0:0.01:10];options = odeset('RelTol',1e-12);[t
