function XVdot = ForS(Time,XV)

global F Omega0 Lambda Omega

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

Xdot = V; Vdot = -2*Lambda*V-Omega0^2*sin(X) + F*cos(Omega*Time);

XVdot = [ Xdot; Vdot ];