In this paper, two numerical methods for solving the MSEIR model are presented. In constructing these methods, the nonstandard finite difference strategy is used. The new methods preserve the qualitative properties of the solution, such as positivity, conservation law, and boundedness. Numerical results are presented to express the efficiency of the new methods.
