Consider the following non-dimensional equations of projectile motion in presence of air resistance,
I understand it can be numerically solved using the mid-point method, but can I also use Verlet method? Why midpoint method is more favourable in this case? Is this because the truncation error in velocity equation is of order $O(t^2)$ and $O(t^3)$ for position when midpoint method is used, which are $O(t^2)$ and $O(t^4)$ respectively when Verlet method is used? Am I thinking in the right direction?