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