function dvdt = ode2(t,v);
% ODE2A: computes dvdt.
g = 9.81; k = 9; m = 1;
dvdt = g-(k/m)*sqrt(v);