Abstract
This article analyses two-dimensional whooping cough model with diffusion term. The resulting two coupled partial differential equations are solved numerically by employing homotopy analysis method (HAM). Convergence of the derived iterative schemes is highlighted and illustrative examples have been presented. It is noticed that the numerical and graphical results are consistent with the dynamical behaviour of the disease.