function XVdot = Force(Time,XV)

global F Omega Lambda Alpha Gamma 

X  = XV(1); V  = XV(2);   

Xdot = V; Vdot = -2*Lambda*V-Omega^2*X-Alpha*X^3 + F*cos(Gamma*Time);

XVdot = [ Xdot; Vdot ];