function [xdot] = sect_3_prob_5_J_ode_fn(t,x, t_f,T,xair,k1,mu)
%
% Checked 2011-11-23 (no errors found)
%
% 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.
% 
%-----

xdot = -(1/T)*( x - xair ) - (mu*(k1^2)/2) * exp( (t-t_f)/T );