Abstract
Whooping cough is a highly transmitted disease around the world. According to the World Health Organization (WHO), 0.15 million cases had reported globally in 2018. Most of the Asian and African states are infected regions. Through the study, we investigated the whole population into the four classes susceptible (S), exposed (E), infected (I), and vaccinated or recovered (R). The transmission dynamics of whooping cough disease are studied analytically and numerically. Analytical analyses are positivity, boundedness, reproduction number, equilibria, and local and global stabilities. In numerical analysis, we developed an implicit numerical integration scheme consistent with the biological problem's properties. The analysis of the implicit method for the said model is dynamically consistent, positive, and bounded. Furthermore, an implicit numerical integration scheme is suitable for studying a particular epidemic model such as the whooping cough disease.