function dvdt = ode2a_global(t,v);
% ODE2A_GLOBAL: computes dvdt. 
global m_value g_value k_value
dvdt = g_value-(k_value)/(m_value)*sqrt(v);