In this paper, a NonStandard Finite Difference (NSFD) scheme is constructed, which can be used to determine 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 difference equations have the same dynamics as the original differential system, such as the positivity of the solutions and the stability of the equilibria, without being restricted by the time step. Our proof of global stability utilizes the method of Lyapunov functions. Numerical simulation illustrates the effectiveness of our results.