Solving a model for the evolution of smoking habit in Spain with homotopy analysis method

We obtain an approximated analytical solution for a dynamic model for the prevalence of the smoking habit in a constant population but with equal and different from zero birth and death rates. This model has been successfully used to explain the evolution of the smoking habit in Spain. By means of the Homotopy Analysis Method, we obtain an analytic expression in powers of time tt which reproduces the correct solution for a certain range of time. To enlarge the domain of convergence we have applied the so-called optimal convergence-control parameter technique and the homotopy-Padé technique. We present and discuss graphical results for our solutions.