fxdot = @(t,x) [x(2); (1/5)*(sin(t)-4*x(1)-7*x(2))];
[t,x] = ode45(fxdot,[0,6],[0,9])
plot(t,x)