% bms925, process modeling 6



tstart = 0;
tend = 100;

x_init = -1;
xdot_init = 1;


[t,v] = ode45(@eqn_forced_osc,[tstart tend],[x_init; xdot_init]);

figure(1)
plot(t,v(:,1))
figure(2)
plot(t,v(:,2));