回复: 【求助】关于ode45的调用
我在使用ode45的时候遇到更多参数的时候直接用global变量,就刚才兄弟的程序:
D2y-u*(1-y^2)*Dy+y=0
function mainfun
global u;
u=10; %这个就是你说的第三个参数
[t,y]=ode45(myfun,[t0,ty],[ya,yb]);
FUNCTION dy=myfun(x,y);
global u;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=u*(1-y(1)^2)*y(2)-y(1);
|