The aim of this paper is to construct a non-standard finite difference (NSFD) scheme that can be used to calculate numerical solutions for an epidemic model with vaccination. Here, the NSFD method is employed to derive a set of difference equations for the epidemic model with vaccination. We show that when, the discretized model preserves positivity and global dynamics of the continuous model. Numerical simulation illustrates the effectiveness of our results. [ABSTRACT FROM PUBLISHER]