m*dv/dt = m*g -k*(v)^2, v(0)=0