function dxdts = sec_ord_forced_F_t(t,x,ft,F_t,params)
% dxdts is the column vector of the first order differential equations
% We pass in the parameters for z and wn in the vector params
z = params(1);
wn = params(2);
F_t = interp1(ft,F_t,t);
dxdts = zeros(2,1);
dxdts(1) = x(2);
dxdts(2) = -(2*z*wn*x(2) + wn^2*x(1)) + F_t;