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
...

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

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