(* Set the values of the various parameters *)
w0 = 1;

(* Solve the differential equations with the given set of initial conditions. *)
sol = NDSolve[
{x'[t] == v[t],
v'[t] == -w0*w0*Sin[x[t]],
x[0] == 0,
v[0] == 1.5},
{x,v},{t,0,50*2*Pi},MaxSteps -> 20000];

(* Plot the solution of the differential equations *)
ParametricPlot[Evaluate[{x[t],v[t]} /. sol],{t, 0, 5*2*Pi},
AxesLabel -> {"x (rad)","v (rad/s)"}]