Impact of general incidence function on three-strain SEIAR model

We investigate the behavior of a complex three-strain model with a generalized incidence rate. The incidence rate is an essential aspect of the model as it determines the number of new infections emerging. The mathematical model comprises thirteen nonlinear ordinary differential equations with susceptible, exposed, symptomatic, asymptomatic and recovered compartments. The model is well-posed and verified through existence, positivity and boundedness. Eight equilibria comprise a disease-free equilibria and seven endemic equilibrium points following the existence of three strains. The basic reproduction numbers $ \mathfrak{R}_{01} $, $ \mathfrak{R}_{02} $ and $ \mathfrak{R}_{03} $ represent the dominance of strain 1, strain 2 and strain 3 in the environment for new strain emergence. The model establishes local stability at a disease-free equilibrium point. Numerical simulations endorse the impact of general incidence rates, including bi-linear, saturated, Beddington DeAngelis, non-monotone and Crowley Martin incidence rates.