function dy=beadf(t,y);
global Omega alpha

dy=[y(2);Omega^2*(cos(y(1))-alpha)*sin(y(1))];