function [xdot] = sect_3_prob_3_J_ode_fn(t,x, omega, alpha_0, alpha_1, alpha_2)
%   
% Written by:
% -- 
% John L. Weatherwax                2006-08-28
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

n = length(x); 
xdot = zeros(n,1); % x = [ v(t), gamma(t) ] 

v     = x(1); 
gamma = x(2); 

xdot(1) = ( 10000 * cos( alpha_0 + alpha_1 * sin( omega * t ) + alpha_2 * cos( omega * t ) ) - 32*sin( gamma ) )/20.;
xdot(2) = ( 10000 * sin( alpha_0 + alpha_1 * sin( omega * t ) + alpha_2 * cos( omega * t ) ) - 32*cos( gamma ) )/(20.*v);