Differential equations using ODE45

please help me in solving this question..
Q. solve:
dy1/dt = y2
dy2/dt = 1000*y2*(1-y1^2)-y1
given initial conditions : y1(0) = 0 and y2(0) = 1.

 采纳的回答

更多回答(1 个)

tspan = [0 20];
Y_int=[0 1];
[t,x] = ode45(@odefcn,tspan,Y_int,[],a, b, c, u);
figure(1)
plot(t,y(:,1),'r',t,y(:,2),'b')
function dydt = odefcn(t,y)
dydt = zeros(2,1);
dydt(1) = y(2);
dydt(2) = 1000*y(2)*(1-y(1)^2)-y(1);
end

2 个评论

Please don't post the full answers to questions that sound like homework assignments, particularly if the poster has not shown what they've tried first.

请先登录,再进行评论。

类别

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by