function xdot = jj2(t,x)
global  betac Ib
%gammac = 1/sqrt(betac);
%xdot = zeros(2,1);
xdot(1) =  (-sin(x(2))./(betac)) -1/betac.*x(1)+Ib/betac ;
xdot(2) = x(1);

