Abstract
The harvested logistic model with a slow variation in coefficients has been considered. Two cases, which depend on the harvest rate, were identified. The first one is when the harvest is subcritical, where the population evolves to an equilibrium. The other is supercritical harvesting, where the population decreases to zero at finite times. The single analytic approximate expression, which is capable of describing both harvesting cases, is readily and explicitly obtained using the multi-time scaling method together with the perturbation approach. This solution fits for a wide range of coefficient values. In addition, such an expression is validated by utilizing numerical computations, which are obtained by using the fourth-order Runge-Kutta technique. Finally, the comparison shows a very good agreement between the two methods.