runge kutta method for differantial equations systems

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)

Answer 1

[Jan]

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