In this work, we investigate the numerical solution of the susceptible exposed infected and recovered measles epidemic model. We also evaluate the numerical stability and the bifurcation value of the transmission parameter from susceptibility to a disease of the proposed epidemic model. The proposed method is a chaos free finite difference scheme, which also preserves the positivity of the solution of the given epidemic model.