4楼: Originally posted by watertxf at 2013-11-27 19:07:38
t的取值范围为1~5...

function solve_ode
clear all;clc
y0=0;
tspan=1:0.01:5;
[ts ys]=ode15s(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


function dydt=odefun(t,y)
dydt=1-(1+16*(1-exp(-t)-y))*y;

2楼: Originally posted by 月只蓝 at 2013-11-27 18:34:09
没有解析解。可以求数值解。

2楼: Originally posted by 月只蓝 at 2013-11-27 18:34:09
没有解析解。可以求数值解。