Numerical Investigation of a Chlamydia Epidemic Model

In this article, we analysed the numerical study of five compartmental epidemic model with Chlamydia infection. The population is divided into five classes: susceptible, exposed, infected in the asymptomatic phase, infected in the symptomatic phase and recovered (SEIAISR). This model was shown to have two equilibrium points: a disease-free equilibrium and an endemic equilibrium. The local stability was defined using the computed effective reproduction number R0, as well as the sensitivity of variables. When R0 < 1, the disease-free equilibrium ( ) is locally asymptotically stable, and if R0 > 1, the endemic equilibrium ( ) is locally asymptotically stable. This model is solved numerically using three numerical techniques: forward Euler, RK-4, and proposed non-standard finite difference (NSFD) techniques. The (NSFD) technique becomes a more efficient and reliable numerical technique than the forward Euler and RK-4 techniques. The NSFD technique retains all essential characteristics of a continuous (SEIAISR) chlamydia epidemic model, like positivity and stability of equilibria. In contrast, well-known forward Euler and RK-4 techniques cannot sustain these characteristics. Furthermore, the NSFD technique is independent of time step size, while forward Euler and RK-4 depend on the time step size. The numerical simulations with a numerical test were represented to validate all the results.