function dv = deriv(t, v, m, cd)
g= 9.81;
dv= g - (cd / m) * v^2;