
function zdot = equation(t,z)
global omega_0 omega beta F0 m
zdot = [z(2); -beta./m*z(2)-omega_0^2*z(1)+F0*cos(omega*t)]; 
end


