10楼: Originally posted by 月只蓝 at 2013-11-27 22:29:46
function solve_ode
clear all;clc
y0=0;
tspan=0:0.01:5;
=ode45(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


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



untitled.jpg
...

10楼: Originally posted by 月只蓝 at 2013-11-27 22:29:46
function solve_ode
clear all;clc
y0=0;
tspan=0:0.01:5;
=ode45(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


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



untitled.jpg
...