function dv_dt = eqn_forced_osc(t,v)

% v(1) = x
% v(2) = dx/dt

dv_dt = [v(2); 
    0.3.*cos(t) - 0.15 .* v(2) + v(1) - v(1).^3];