“…Different simulations can be carried out using various values of parameters. In the present numerical approach, we use the following parameters values taken from [9]: α 1 = 0.0023, α 2 = 0.0015, α 3 = 0.0018, α 4 = 0.00013, β 1 = 0.000022, β 2 = 0.000026, β 3 = 0.000027, β 4 = 0.000024, β 5 = 0.000025, λ 1 = 0.00003, λ 2 = 0.00002, λ 3 = 0.00018, λ 4 = 0.00004, δ 1 = 0.00002, δ 2 = 0.00002, δ 3 = 0.00022, δ 4 = 0.000028, µ 1 = 0.0002, µ 2 = 0.0003, µ 3 = 0.0004, µ 4 = 0.0004, γ 1 = 0.00035, γ 2 = 0.0001, γ 3 = 0.0001, τ 1 = 2000, and τ 2 = 3290.…”