Abstract
This work is devoted to find the reliable numerical solution of an epidemic model with constant vaccination strategy. For this purpose, a structure preserving numerical scheme called the Grunwald-Letnikov nonstandard finite difference scheme is designed. The proposed technique retains all the important properties of the continuous epidemic model like boundedness, positivity, and stability. This behavior of the proposed numerical scheme is validated mathematically and graphically. The role of the vaccination in controlling the disease dynamics in the population is verified through numerical simulations. The stability of the system under discussion is also examined at the disease free equilibrium point and the endemic equilibrium point. Finally, the outcome of this study is furnished with concluding remarks and future directions of research.