查看单个帖子
旧 2009-02-25, 16:37   #3
carptin
初级会员
 
注册日期: 2008-12-14
年龄: 50
帖子: 6
声望力: 0
carptin 正向着好的方向发展
默认 回复: 【求助】关于ode45的调用

今天心情好,给你个例子好好看看 三个参数的
方程D2y-u*(1-y^2)*Dy+y=0
FUNCTION dy=myfun(x,y,u);
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=u*(1-y(1)^2)*y(2)-y(1);

如果求u=10,如下:
ode45(@myfun,[0 20],[2;0],[],10)
carptin 当前离线   回复时引用此帖