Abstract
This work attempts to study the numerical solution of nonlinear delayed Immunized Susceptible Latent Infected and Recovered (MSLIR) epidemic model of HBV disease. Reproduction number, equilibria and stability are discussed. Three different numerical techniques, Euler, RK-4 and the non-standard finite difference (NSFD) techniques are used for the numerical solution of the model. The proposed technique is independent of the size of the time step, while forward Euler and RK-4 depend on the size of a time step and retains all essential characteristics of the continuous MSLIR epidemic model like positivity and stability of equilibrium, while well-known forward Euler and RK-4 cannot sustain these characteristics. Therefore, the proposed (NSFD) technique becomes a more efficient and reliable numerical technique than the forward Euler and RK-4 scheme. Numerical simulations are presented for the validation of the obtained results.