“…The data used for the computational simulations are extracted from the literature or assumed. The parameter values for the computational simulations are (α * ) α = 2 (assumed), (β * ) α = 0.2 (assumed), M α S = 667.685 [28], M α D = 4.1 (assumed), β α 1 = 0.25 (assumed), ε = 0.007 [10], η α = 0.1 [10], γ α = 0.0015 [10], δ α = 0.002 [10], µ α = 1 70.5 [26], d α = 0.07 (assumed), t H = 1.02 (assumed), t D = 1.03 (assumed), µ α d = 0.013 (assumed), α α 1 = 0.1 (assumed), α α 2 = 0.35 [45,46], α α 3 = 0.4 [45,46], and α = 0.3, 0.5, 0.7, 0.9, 1. The initial conditions are S = 874.1400, O w = 1.2000, O b = 1.5000 and D = 100.0000 on the scale of 10000 individuals.…”