A convenient approximated analytic solution is proposed for the problem of the motion of a body under a resistive force, acting in the magnitude of the squared velocity of the body. This solution is an explicit function of time, that keeps a good behavior both near the initial state and far from the initial state. To obtain a general analytic solution, we firstly used a reduction principle to be able to manipulate scalar objects, and we analyzed limit behaviors, both near the initial state and far from the initial state. Secondly, we proposed an approximated analytic solution with heuristics based on the built knowledge. Finally, a robust and stable integration scheme is proposed, based on the obtained analytic solution. We compared the scheme with other standard integration schemes.