function dzdt = Damped_nonlin_sys(t,z,m,f,k1,k3)
%Damped_nonlin_sys computes the derivative of the states z
dzdt(1,1) = z(2);
dzdt(2,1) = -(f*z(2)+k1*z(1)+k3*z(1)^3)/m;
end

