runge kutta method for differantial equations systems

4 views (last 30 days)
How can I solve following differential equations systems and display the charts for the variables y1 and y2 as a function of time for t between 0 to 0.1 s? (with runge kutta method)
dy1/dt = (10sin(314t)-5y1-y2) / 0.01
dy2/dt = 1000y1
dp:the initial conditions are written yi(0)

Accepted Answer

Jan
Jan on 9 Feb 2021
Edited: Jan on 9 Feb 2021
y0 = [1, 2];
[T, Y] = ode45(@OdeFcn, [0, 0.1], y0);
plot(T, Y);
function dy = OdeFcn(t,y)
dy = [(10 * sin(314 * t) - 5 * y(1) - y(2)) / 0.01; ...
1000 * y(1)];
end

More Answers (0)

Community Treasure Hunt

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

Start Hunting!

Translated by