ode45 second order
26 views (last 30 days)
Show older comments
please, I want to plot second order and use function ode45. So just now I can't write function.
y" = 3.6sin(t) - 60y - 5y'
and this is my question 👆👆
function dy=model(t,y)
dy = [y(2); 3.6*sin(t) - 60*y(0) - 5*y(1)]
end
and this 👆👆 is my error function
0 Comments
Answers (1)
Alan Stevens
on 18 Apr 2021
Matlab indices start at 1, so you want
dy = [y(2); 3.6*sin(t) - 60*y(1) - 5*y(2)];
3 Comments
Alan Stevens
on 18 Apr 2021
"how can i add sin input to this order"
I'm not sure what you mean. Your second order ode already has 3.6*sin(t) as a forcing function. Do you mean you want to add another term to the forcing function e.g.
y'' + 5*y' + 60*y = 3.6*sin(t) + 2*sin(2*t)
If so, then just make
dy = [y(2); 3.6*sin(t) + 2*sin(2*t) - 60*y(1) - 5*y(2)];
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!