Advanced optimal control approaches for immune boosting and clinical treatment to enhance dengue viremia models using ABC fractional-order analysis

Introduction: This work focuses on the Dengue-viremia ABC (Atangana-Baleanu Caputo) fractional-order differential equations, accounting for both symptomatic and asymptomatic infected cases. Symptomatic cases are characterized by higher viremia levels, whereas asymptomatic cases exhibit lower viremia levels. The fractional-order model highlights memory effects and other advantages over traditional models, offering a more comprehensive representation of dengue dynamics.

Methods: The total population is divided into four compartments: susceptible, asymptomatic infected, symptomatic infected, and recovered. The model incorporates an immune-boosting factor for asymptomatic infected individuals and clinical treatment for symptomatic cases. Positivity and boundedness of the model are validated, and both local and global stability analyses are performed. The novel Adams-Bash numerical scheme is utilized for simulations to rigorously assess the impact of optimal control interventions.

Results: The results demonstrate the effectiveness of the proposed control strategies. The reproduction numbers must be reduced based on specific optimal control conditions to effectively mitigate disease outbreaks. Numerical simulations confirm that the optimal control measures can significantly reduce the spread of the disease.

Discussion: This research advances the understanding of Dengue-viremia dynamics and provides valuable insights into the application of ABC fractional-order analysis. By incorporating immune-boosting and clinical treatment into the model, the study offers practical guidelines for implementing successful disease control strategies. The findings highlight the potential of using optimal control techniques in public health interventions to manage disease outbreaks more effectively.

1 Introduction

Worldwide, thousands of dengue cases are reported every year. The world's tropical and subtropical regions are affected by dengue infection, which is a mosquito-carrying disease. A high temperature and flu-like symptoms are signs of mild illness or asymptomatic to stern disease. DHF (Dengue Hemorrhagic Fever) or DSS (Dengue Show Syndromes Syndromes) is a highly infectious form of dengue fever that causes serious bleeding, shock, and death. Generally, it was noticed that only one out of four dengue contagions is symptomatic. Dengue virus occurs in four major types (DENV types 1, 2, 3, and 4), all of which can cause serious illness. The usual signs of DENV type 1 are like a common cold and mild fever, which will not lead directly to DHF; conversely, later DENV types can lead to DHF (1–3).

To understand the dynamical behavior of dengue transmission, we formulated a mathematical model, particularly focusing on vector-borne disease transmission from mosquitoes to humans. Esteva and Vargas (4, 5) pioneered the creation of a fundamental dengue model and explored numerous fundamental mathematical concepts and their accompanying numerical simulations. Feng et al. (6) presented a two-strain dengue infection model and examined competitive exclusion. Researchers have conducted numerous studies to better understand the transmission of dengue fever (7–9).

The importance of fractional-order models lies in their ability to capture the complex dynamics and long-term dependencies within the transmission process. By incorporating fractional derivatives, these models provide a more comprehensive understanding of disease spread, which is crucial for designing effective intervention strategies. The fractional-order models can accommodate the nuanced behavior of dengue transmission, offering insights that integer-order models may overlook, thereby enhancing the accuracy and effectiveness of disease control measures.

The fractional order model has been conclusively demonstrated by a recent study to be capable of controlling the trend of complex diffusion disorder (10–14). Many have emphasized various mathematical models for Dengue transmission and prevention (15–19). All cited references explain the transmission process of Dengue infection from different perspectives, including dynamic analysis, evaluation of vaccination, and optimal control measures (20–24). The most updated studies on Dengue with real-life data are presented in (25, 26). The mathematical description of Dengue is briefly described in Deterministic and Stochastic terms. The evolution of dengue with asymptomatic carriers using optimal control measures was investigated in (27).

Therefore, motivated by the aforementioned literature, we propose a computational framework for the dissemination of dengue at a given viremia level. We investigated whether symptom-free people were markedly more susceptible to mosquitoes than clinically symptom-positive patients. The new idea of a mathematical model to analyse the immune-boosting factor for asymptomatic infected cases and the waning immunity that cases re-infect is reported. To make practical applications and simulations easier, we utilize the Adams-Bash forth numerical scheme, which is renowned for its accuracy and stability. This choice ensures that our model reflects real-world scenarios while maintaining computational efficacy. A key highlight of this study is the incorporation of optimal control strategies into the ABC fractional order Dengue viremia model. These strategies are designed to explore how interventions, such as self-prevention and vector control, can be optimized to curtail disease spread. The analysis extends to investigating disease-free and endemic stability, providing crucial insights into the long-term behavior of the system under various control scenarios.

This article is prepared as follows: In portion 2, we review the fundamental definitions for the fractional-order operator and provide a list of mathematical properties that were used throughout the work. The dengue viral mathematical model with fractional order was presented in portion 3. Portion 4 examines the local as well as global consistency of the suggested model through the Routh-Hurwitz criteria and the Lyapunov function. An optimal control solution and discussion are present in portion 5. The final section focuses on numerical simulations and a comprehensive conclusion.

2 Fundamental results

This section introduces fractional derivation and some of its properties, which will be used in the following components.

Definition 2.1.

Consider ψ ∈ ℍ′ (0, T) and η ∈ [0, ȶ], then Atangana-Baleanu fraction component in Caputo case is

$\begin{array}{ll}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }\psi \left(ȶ\right)=\frac{\overline{\mathcal{N}\left(\eta \right)}}{1-\eta }{\displaystyle {\int }_0^{ȶ}}\frac{d}{dx}\psi \left(\mathcal{S}\right){\mathcal{N}}_{\eta }\left[-\frac{\eta }{1-\eta }\left(ȶ-\mathcal{S}\right)\right]d\mathcal{S}& (1)\end{array}$

The method yields a variation operator Caputo-Fabrizio that replaces

$\begin{array}{l}{\mathcal{N}}_{\eta }\left[-\frac{\eta }{1-\eta }\left(ȶ-\mathcal{S}\right)\right]d\mathcal{S}\text{ by }{\mathcal{N}}_1=exp\left[-\frac{\eta }{1-\eta }\left(ȶ-\mathcal{S}\right)\right]\end{array}$

It's noteworthy that $\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }$[constant] = 0. Here $\overline{\mathcal{N}\left(\eta \right)}$ is the typical function and it is defined as $\overline{\mathcal{N}\left(0\right)}=1$ and $\overline{\mathcal{N}\left(1\right)}=1$. $\overline{\mathcal{N}\left(\eta \right)}$ depict the familiar Mittag–Leffler operator, it also reflects the exponential function generality.

Definition 2.2.

The fractional integral of $\mathcal{ABC}$ with order η given by

$\begin{array}{ll}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔗}_{ȶ}^{\eta }\psi \left(ȶ\right)=\frac{1-\eta }{\overline{\mathcal{N}\left(\eta \right)}}\psi \left(ȶ\right)+\frac{\eta }{\overline{\mathcal{N}\left(\eta \right)\Gamma \left(\eta \right)}}{\displaystyle {\int }_0^{ȶ}}\psi \left(\mathcal{S}\right){\left(ȶ-\mathcal{S}\right)}^{\eta -1}d\mathcal{S}& (2)\end{array}$

Lemma 2.1.

Consider a fractional-order system

$\begin{array}{ll}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }\mathcal{x}\left(ȶ\right)=\mathcal{f}\left(ȶ,\mathcal{x}\right),t>t_0& (3)\end{array}$

Where ηϵ(0, 1) in the initial case $\mathcal{x}\left({ȶ}_0\right).$

If $\mathcal{f}\left(ȶ,\mathcal{x}\right)$ fulfills the Lipschitz condition in relation to x, then system (Equation 3) exhibits a unique solution in the region [t₀, +∞) × φ and φ ⊆ ℝⁿ.

Lemma 2.2.

If $\mathcal{x}\left(ȶ\right)\in {ℝ}^{+}$ become an ongoing and attainable consequence. Then

$\begin{array}{ll}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }\left(\mathcal{x}\left(ȶ\right)-{\mathcal{x}}^{*}-{\mathcal{x}}^{*}ln\frac{\mathcal{x}\left(ȶ\right)}{{\mathcal{x}}^{*}}\right)\le \left(1-\frac{{\mathcal{x}}^{*}}{\mathcal{x}\left(ȶ\right)}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }\mathcal{x}\left(ȶ\right)& (4)\end{array}$

Here t > t₀, η ϵ (0, 1) and ${\mathcal{x}}^{*}\in {ℝ}^{+}.$

3 Evaluation of dengue dynamics

In this section, we expand upon the previously described Dengue SIR-SI model (18) by incorporating additional factors and refining the classification of both human and mosquito populations. Our model includes viremia levels, an immune-boosting factor for asymptomatic infected cases, and clinical treatment for symptomatic infected cases.

To study the mode of spread of dengue sickness, the human species (${\mathcal{N}}_{𝔥}$) is subdivided into four classes: susceptible (Ȿ_𝔥), symptomatic infectious (𝔗_Ȿ𝔥), asymptomatic infectious $\left({𝔗}_{{\mathcal{A}}_{𝔥}}\right)$ and recovered human populations (𝔎_𝔥). We classified female mosquito species $\left({\mathcal{N}}_{𝔪}\right)$ into Susceptible (Ȿ_𝔪) and infective mosquitoes (𝔗_𝔪). A Susceptible individual among as one who is not infected and immune, infected humans are both asymptomatic and symptomatic are those who have acquired Dengue viremia from an infected mosquito populations and are all capable of spreading dengue virus to susceptible mosquitoes. Let we examines π_𝔥 and π_𝔪 the acquisition rates of humans and mosquitoes. The proposed model, illustrated in the flowchart, demonstrates the dengue transmission dynamics. Based on Figure 1, we developed the following differential equation.

$\begin{array}{l}\frac{d{Ȿ}_{𝔥}}{dȶ}={\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\\ \frac{d{𝔗}_{{\text{s}}_{\text{h}}}}{dȶ}={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\\ \frac{{d𝔗}_{{\mathcal{A}}_{𝔥}}}{dȶ}={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ \frac{d{𝔎}_{𝔥}}{dȶ}=\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\\ \frac{{dⱾ}_{𝔪}}{dȶ}={\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\\ \frac{{d𝔗}_{𝔪}}{dȶ}={\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}& (5)\end{array}$

Where

α_𝔪- individual mosquito‘s biting rate

χ_Ȿ𝔥- Dissemination to human by mosquitoes, which leads to a symptomatic infectious in humans

${\chi }_{{\mathcal{A}}_{𝔥}}$- Dissemination to human by mosquitoes, which leads to a asymptomatic infectious in humans

χ_𝔪 - Viremia dissemination to mosquito by human species

δ_𝔥 - Human Fatality rate

τ - Symptomatic infected human treatment rate

γ - Recovering rate.

θ - Transition rate at which a recovered person becomes defenseless due to loss of immunity

ϱ − Rates of immunosuppression for asymptomatic victims

δ_𝔪 - Rate of mosquito natural mortality (an average mosquito life span)

From the basic cases $({Ȿ}_h,T_{{Ȿ}_h},T_{{\mathcal{A}}_h},K_h,{Ȿ}_m,T_m)\ge 0.$

Figure 1

Figure 1. The process diagram in dengue dynamics.

In this approach, the aggregate human and mosquito population ratios are provided by

$\begin{array}{l}{\mathcal{N}}_{𝔥}={Ȿ}_{𝔥}+{𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}+{𝔎}_{𝔥}\text{ and }{\mathcal{N}}_{𝔪}={Ȿ}_{𝔪}+{𝔗}_{𝔪}.\end{array}$

In addition, the area of biologically significance for the aforementioned dengue model is indicated and presented by the covered set

$\begin{array}{c}\Phi =\{{Ȿ}_{𝔥},{𝔗}_{{Ȿ}_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}},{𝔎}_{𝔥},{Ȿ}_{𝔪},{𝔗}_{𝔪}\in {ℝ}_{+}^6:{Ȿ}_{𝔥}+{𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}+{𝔎}_{𝔥}\\ \le {\mathcal{N}}_{𝔥},{Ȿ}_{𝔪}+{𝔗}_{𝔪}\le {\mathcal{N}}_{𝔪}\}\end{array}$

A fractional representation of the 𝔄𝔅ℭ model as

$\begin{array}{l}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔥}={\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{Ȿ}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{\mathcal{A}}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔎}_{𝔥}=\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔪}={\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{𝔪}={\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}& (6)\end{array}$

A fractional derivation of Atangana-Baleanu of order 0 < η < 1 is denoted $\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }$ in Caputo notation.

4 Model analysis

This section examines the validity, singularity and positive variance of the solution of the SIR-SI type model. Additionally, a reliability estimate for Model (Equation 6) has also been developed.

4.1 Existence and uniqueness

Theorem 4.1.

For each non- negative initial stage $({Ȿ}_{𝔥}\left(0\right),{𝔗}_{{Ȿ}_{𝔥}}\left(0\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(0\right)$, ${𝔎}_{𝔥}\left(0\right),{Ȿ}_{𝔪}\left(0\right),{𝔗}_{𝔪}\left(0\right))\in {ℝ}_{+}^6,$ then there survives a oneness solution of fractional order model (Equation 6).

Proof

Let $\Phi =\{{Ȿ}_{𝔥},{𝔗}_{{Ȿ}_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}},{𝔎}_{𝔥},{Ȿ}_{𝔪},{𝔗}_{𝔪}\in {ℝ}_{+}^6:max(|{Ȿ}_{𝔥}|$, $\left|{𝔗}_{{Ȿ}_{𝔥}}\right|,\left|{𝔗}_{{\mathcal{A}}_{𝔥}}\right|,\left|{𝔎}_{𝔥}\right|,\left|{Ȿ}_{𝔪}\right|,\left|{𝔗}_{𝔪}\right|)\le \epsilon \}$.

Define a mapping

$𝕄\left(\mathcal{x}\right)=\left({𝕄}_1\left(\mathcal{x}\right),{𝕄}_2\left(\mathcal{x}\right),{𝕄}_3\left(\mathcal{x}\right),{𝕄}_4\left(\mathcal{x}\right),{𝕄}_5\left(\mathcal{x}\right){,𝕄}_6\left(\mathcal{x}\right)\right)$ and

$\begin{array}{l}{𝕄}_1\left(\mathcal{x}\right)={\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\\ {𝕄}_2\left(\mathcal{x}\right)={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\\ {𝕄}_3\left(\mathcal{x}\right)={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ {𝕄}_4\left(\mathcal{x}\right)=\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\\ {𝕄}_5\left(\mathcal{x}\right)={\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\\ {𝕄}_6\left(\mathcal{x}\right)={\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}\end{array}$

Where $\mathcal{x}=\left({Ȿ}_{𝔥},{𝔗}_{{Ȿ}_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}},{𝔎}_{𝔥},{Ȿ}_{𝔪},{𝔗}_{𝔪}\right)\in \Phi$

For any $\mathcal{x},\overline{\mathcal{x}}\in \Phi ,$ we have

$\begin{array}{l}\left|\right|𝕄\left(\mathcal{x}\right)-\overline{𝕄\left(\mathcal{x}\right)}\left|\right|=|{𝕄}_1\left(\mathcal{x}\right)-\overline{{𝕄}_1\left(\mathcal{x}\right)}|+|{𝕄}_2\left(\mathcal{x}\right)-\overline{{𝕄}_2\left(\mathcal{x}\right)}|\\ +|{𝕄}_3\left(\mathcal{x}\right)-\overline{{𝕄}_3\left(\mathcal{x}\right)}|+|{𝕄}_4\left(\mathcal{x}\right)-\overline{{𝕄}_4\left(\mathcal{x}\right)}|\\ +|{𝕄}_5\left(\mathcal{x}\right)-\overline{{𝕄}_5\left(\mathcal{x}\right)}|+|{𝕄}_6\left(\mathcal{x}\right)-\overline{{𝕄}_6\left(\mathcal{x}\right)}|\\ \le |{\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}-{\pi }_{𝔥}\\ +{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)\overline{{Ȿ}_{𝔥}}\overline{{𝔗}_{𝔪}}+{\delta }_{𝔥}\overline{{Ȿ}_{𝔥}}\theta \overline{{𝔎}_{𝔥}}|\\ +|{\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}-{\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}\overline{{Ȿ}_{𝔥}}\overline{{𝔗}_{𝔪}}\\ +\left(\tau +\gamma +{\delta }_{𝔥}\right)\overline{{𝔗}_{{Ȿ}_{𝔥}}}|\\ +|{\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}-{\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}\overline{{Ȿ}_{𝔥}}\overline{{𝔗}_{𝔪}}\\ +\left(\varrho +\gamma +{\delta }_{𝔥}\right)\overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}|\\ +|\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\\ -\gamma \left(\overline{{𝔗}_{{Ȿ}_{𝔥}}}+\overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)+\theta \overline{{𝔎}_{𝔥}}+{\delta }_{𝔥}\overline{{𝔎}_{𝔥}}-\tau \overline{{𝔗}_{{Ȿ}_{𝔥}}}-\varrho \overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}|\\ +|{\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}-{\pi }_{𝔪}+{\alpha }_{𝔪}{\chi }_{𝔪}\\ \left(\overline{{𝔗}_{{Ȿ}_{𝔥}}}+\overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)\overline{{Ȿ}_{𝔪}}+{\delta }_{𝔪}\overline{{Ȿ}_{𝔪}}|+|{\alpha }_{𝔪}{\chi }_{𝔪}\\ \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left(\overline{{𝔗}_{{Ȿ}_{𝔥}}}+\overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)\overline{{Ȿ}_{𝔪}}+{\delta }_{𝔪}\overline{{𝔗}_{𝔪}}|.\\ \le \left[2{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)\mathcal{P}+{\delta }_{𝔥}\right]|{Ȿ}_{𝔥}-\overline{{Ȿ}_{𝔥}}|+\left(2\tau +2\gamma +{\delta }_{𝔥}\right)\\ |{𝔗}_{{Ȿ}_{𝔥}}-\overline{{𝔗}_{{Ȿ}_{𝔥}}}|+\left(2\varrho +2\gamma +{\delta }_{𝔥}\right)\\ |{𝔗}_{{\mathcal{A}}_{𝔥}}-\overline{{𝔗}_{{\mathcal{A}}_{𝔥}}}|+\left(2\theta {+\delta }_{𝔥}\right)|{𝔎}_{𝔥}-\overline{{𝔎}_{𝔥}}|+\left[2{\alpha }_{𝔪}\left(\mathcal{Q}+\mathcal{R}\right){\chi }_{𝔪}+{\delta }_{𝔥}\right]\\ |{Ȿ}_{𝔪}-\overline{{Ȿ}_{𝔪}}|+{\delta }_{𝔪}|{𝔗}_{\text{m}}-\overline{{𝔗}_{\text{m}}}|\mathcal{K}\le \left|\right|\mathcal{x}-\overline{\mathcal{x}}\left|\right|.\end{array}$

Where $\mathcal{K}=max\left[\begin{array}{c}2{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)\mathcal{P}+{\delta }_{𝔥},\\ 2\left(\tau +\gamma \right)+{\delta }_{𝔥},2\left(\varrho +\gamma \right)+{\delta }_{𝔥},2\theta {+\delta }_{𝔥},\\ 2{\alpha }_{𝔪}\left(\mathcal{Q}+\mathcal{R}\right){\chi }_{𝔪}+{\delta }_{𝔥},{\delta }_{𝔪}\end{array}\right]$

Basically, since $𝕄\left(\mathcal{x}\right)$ satisfies the Lipschitz requirement. Model (Equation 6) has a singular solution based on Lemma 1.

4.2 Positivity solution

Since system (Equation 6) deals with mosquitoes and populace, all components of system are positive. Following is our discussion:

Theorem 4.2.

Let (Ȿ_𝔥, 𝔗_Ȿ𝔥, ${𝔗}_{{\mathcal{A}}_{𝔥}}$, 𝔎_𝔥, Ȿ_𝔪, 𝔗_𝔪) > 0 be represent the system (Equation 6) solution for the primary points ${Ȿ}_{𝔥}\left(0\right),{𝔗}_{{Ȿ}_{𝔥}}\left(0\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(0\right),{𝔎}_{𝔥}\left(0\right),{Ȿ}_{𝔪}\left(0\right),{𝔗}_{𝔪}\left(0\right)$ and represents an immutable set

$\Phi =\left\{{Ȿ}_{𝔥},{𝔗}_{{Ȿ}_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}},{𝔎}_{𝔥},{Ȿ}_{𝔪},{𝔗}_{𝔪}\in {ℝ}_{+}^6:{\mathcal{N}}_{𝔥}=\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}},{\mathcal{N}}_{𝔪}=\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}\right\}$, then, all elements of the closed set

Φ is traveling in ${ℝ}_{+}^6$ space is positive invariant.

Proof

The given equation is used to construct the Lyapunov function:

$\begin{array}{l}𝕃\left(ȶ\right)=\left({𝕃}_1\left(ȶ\right),{𝕃}_2\left(ȶ\right)\right)=\left({Ȿ}_{𝔥}+{𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}+{𝔎}_{𝔥},{Ȿ}_{𝔪}+{𝔗}_{𝔪}\right)\end{array}$

The function 𝕃(ȶ) satisfies

$\begin{array}{r}\dot{𝕃\left(ȶ\right)}=(\dot{{𝕃}_1\left(ȶ\right),}\dot{{𝕃}_2\left(ȶ\right),})=\left(\dot{{Ȿ}_{𝔥}}+\dot{{𝔗}_{{Ȿ}_{𝔥}}}+\dot{{𝔗}_{{\mathcal{A}}_{𝔥}}}+\dot{{𝔎}_{𝔥}},\dot{{Ȿ}_{𝔪}}+\dot{{𝔗}_{𝔪}}\right)\\ =\left({\pi }_{𝔥}-{\delta }_{𝔥}{Ȿ}_{𝔥}-{𝔗}_{{Ȿ}_{𝔥}}{\delta }_{𝔥}-{𝔗}_{{\mathcal{A}}_{𝔥}}{\delta }_{𝔥}{-\delta }_{𝔥}{𝔎}_{𝔥},{\pi }_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}\right)\\ =\left({\pi }_{𝔥}-{\delta }_{𝔥}{𝕃}_1,{\pi }_{𝔪}-{\delta }_{𝔪}{𝕃}_2\right)& (7)\end{array}$

Therefore, it is simple to demonstrate Equation 7 as regards:

$\begin{array}{ll}\{\begin{array}{c}\dot{{𝕃}_1\left(ȶ\right)}={\pi }_{𝔥}-{\delta }_{𝔥}{𝕃}_1\le 0,for{𝕃}_1\ge \frac{{\pi }_{𝔥}}{{\delta }_{𝔥}}\\ \dot{{𝕃}_2\left(ȶ\right)}={\pi }_{𝔪}-{\delta }_{𝔪}{𝕃}_2\le 0,for{𝕃}_2\ge \frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}\end{array}& (8)\end{array}$

Inferring $\dot{𝕃\left(ȶ\right)}\le 0$ from the above equations, which indicates that f is positively stable collection. On the other hand, by solving system (Equation 6)

$\begin{array}{l}0\le \left({𝕃}_1\left(ȶ\right),{𝕃}_2\left(ȶ\right)\right)<\left(\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}}+{𝕃}_1\left(0\right)e^{-{\delta }_{𝔥}t},\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}+{𝕃}_2\left(0\right)e^{-{\delta }_{𝔪}t}\right)\end{array}$

Where 𝕃₁(0) and 𝕃₂(0) are the primary states of 𝕃₁(ȶ) and 𝕃₂(ȶ) respectively. Therefore, t → ∞, $0\le \left({𝕃}_1\left(ȶ\right),{𝕃}_2\left(ȶ\right)\right)\le \left(\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}},\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}\right)$ and we can conclusion that Φ is a desirable set.

This establishes the theorem.

4.3 Basic reproduction value ${\mathcal{R}}_0$

Let ${\mathcal{C}}_{\mathcal{f}}=\left({Ȿ}_{𝔥}^{*},{𝔗}_{{Ȿ}_{𝔥}}^{*},{𝔗}_{{\mathcal{A}}_{𝔥}}^{*},{𝔎}_{𝔥}^{*},{Ȿ}_{𝔪}^{*},{𝔗}_{𝔪}^{*}\right)$ be the contagious free equilibrium of Equation 6. We have ${\mathcal{C}}_{\mathcal{f}}=\left(\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}},0,0,\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}},0\right)$. The algorithm of the next iteration matrix is utilized to estimate ${\mathcal{R}}_0$. Obviously, the infected compartments are 𝔗_Ȿ𝔥, ${𝔗}_{{\mathcal{A}}_{𝔥}}$ and 𝔗_𝔪 as a consequence of Equation 6. There are

$\begin{array}{l}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{Ȿ}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{\mathcal{A}}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{𝔪}={\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}& (9)\end{array}$

Then we derive

$\begin{array}{l}\mathcal{F}=\left(\begin{array}{ccc}0& 0& {\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}\\ 0& 0& {\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}\\ {\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& {\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& 0\end{array}\right),\\ \mathcal{V}=\left(\begin{array}{ccc}\tau +\gamma +{\delta }_{𝔥}& 0& 0\\ 0& \varrho +\gamma +{\delta }_{𝔥}& 0\\ 0& 0& {\delta }_{𝔪}\end{array}\right)\end{array}$

The basic reproduction value is given by

$\begin{array}{l}{\mathcal{R}}_0=\rho \left({\mathcal{FV}}^{-1}\right)=\frac{{\pi }_{𝔥}{\pi }_{𝔪}{\chi }_{𝔪}{{\alpha }_{𝔪}}^2\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)}{{\delta }_{𝔥}{\left({\delta }_{𝔪}\right)}^2\left(\tau +\gamma +{\delta }_{𝔥}\right)\left(\varrho +\gamma +{\delta }_{𝔥}\right)}\end{array}$

Where $\rho \left({\mathcal{FV}}^{-1}\right)$ denote the spectral radius. Surmise that ${\mathcal{C}}_{\mathcal{p}}=\left({Ȿ}_{𝔥}^{**},{𝔗}_{{Ȿ}_{𝔥}}^{**},{𝔗}_{{\mathcal{A}}_{𝔥}}^{**},{𝔎}_{𝔥}^{**},{Ȿ}_{𝔪}^{**},{𝔗}_{𝔪}^{**}\right)$ represents the endemic equilibrium for Equation 6. So that

$\begin{array}{l}{Ȿ}_{𝔥}^{**}=\frac{{\pi }_{𝔥}\left(\theta +{\delta }_{𝔥}\right)u_1{u}_2}{\left({\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){𝔗}_{𝔪}+{\delta }_{𝔥}\right)u_1{u}_2\left(\theta +{\delta }_{𝔥}\right)-\theta {\alpha }_{𝔪}{{𝔗}_{𝔪}u}_4}\\ {𝔗}_{{Ȿ}_{𝔥}}^{**}=\frac{{\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}}{u_1}\\ {𝔗}_{{\mathcal{A}}_{𝔥}}^{**}=\frac{{\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}}{u_2}\\ {𝔎}_{𝔥}^{**}=\frac{{\alpha }_{𝔪}{Ȿ}_{𝔥}{𝔗}_{𝔪}{u}_4}{{\left(\theta +{\delta }_{𝔥}\right)u}_1{u}_2}\\ {Ȿ}_{𝔪}^{**}=\frac{u_1{u}_2{\pi }_{𝔪}}{{{\alpha }_{𝔪}}^2{\chi }_{𝔪}{Ȿ}_{𝔥}{𝔗}_{𝔪}{u}_3+{\delta }_{𝔪}{u}_1{u}_2}\\ {𝔗}_{𝔪}^{**}=\frac{{\left(\theta +{\delta }_{𝔥}\right)u}_3{\alpha }_{𝔪}{\pi }_{𝔥}{Ȿ}_{𝔪}-{\delta }_{𝔥}{u}_1{u}_2{\delta }_{𝔪}}{\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)u_5-\theta u_4}\end{array}$

Where $u_1=\tau +\gamma +{\delta }_{𝔥},u_2=\varrho +\gamma +{\delta }_{𝔥},u_3=u_1{\chi }_{{\mathcal{A}}_{𝔥}}+u_2{\chi }_{{Ȿ}_{𝔥}},u_4=u_1{\chi }_{{\mathcal{A}}_{𝔥}}\left(\tau +\gamma \right)+u_2{\chi }_{{Ȿ}_{𝔥}}\left(\varrho +\gamma \right),$ and u₅ = u₁u₂ (θ + δ_𝔥) δ_𝔪.

4.4 Local stability

In this part, we are covering the analysis of firmness conditions of contagious free equilibrium ${\mathcal{C}}_{\mathcal{f}}$ and contagious persistence equilibrium ${\mathcal{C}}_{\mathcal{p}}$ points. A steady state analysis of this equilibrium results in the following Theorem 4.3 and Theorem 4.4.

The obtained Jacobian matrix is:

$\begin{array}{ll}\Im =\left(\begin{array}{cccccc}{-\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){𝔗}_{𝔪}-{\delta }_{𝔥}& 0& 0& \theta & 0& {-\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}\\ {\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{𝔗}_{𝔪}& -\left(\tau +\gamma +{\delta }_{𝔥}\right)& 0& 0& 0& {\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}\\ {\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{𝔗}_{𝔪}& 0& -\left(\varrho +\gamma +{\delta }_{𝔥}\right)& 0& 0& {\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}\\ 0& \tau +\gamma & \varrho +\gamma & -\left(\theta +{\delta }_{𝔥}\right)& 0& 0\\ 0& {-\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& -{\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& 0& -\left({\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)+{\delta }_{𝔪}\right)& 0\\ 0& {\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& {\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}& 0& {\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)& -{\delta }_{𝔪}\end{array}\right)& (10)\end{array}$

Theorem 4.3.

If ${\mathcal{R}}_0<1,$ the non - contagious equilibrium ${\mathcal{C}}_{\mathcal{f}}$ is locally stable.

Proof

The structure (Equation 6) in the Jacobian matrix of ${\mathcal{C}}_{\mathcal{f}}$ follows

$\begin{array}{ll}\Im \left({\mathcal{C}}_{\mathcal{f}}\right)=\left(\begin{array}{cccccc}-{\delta }_{𝔥}& 0& 0& \theta & 0& {-\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}}\\ 0& -\left(\tau +\gamma +{\delta }_{𝔥}\right)& 0& 0& 0& {\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}}\\ 0& 0& -\left(\varrho +\gamma +{\delta }_{𝔥}\right)& 0& 0& {\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}}\\ 0& \tau +\gamma & \varrho +\gamma & -\left(\theta +{\delta }_{𝔥}\right)& 0& 0\\ 0& {-\alpha }_{𝔪}{\chi }_{𝔪}\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}& {-\alpha }_{𝔪}{\chi }_{𝔪}\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}& 0& {\delta }_{𝔪}& 0\\ 0& {\alpha }_{𝔪}{\chi }_{𝔪}\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}& {\alpha }_{𝔪}{\chi }_{𝔪}\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}}& 0& 0& -{\delta }_{𝔪}\end{array}\right)& (11)\end{array}$

To determine the eigenvalue from the above-described matrix $det\left(\Im \left({\mathcal{C}}_{\mathcal{f}}\right)-\lambda 𝔗\right)=0$

We obtain the Eigen values λ₁ = −δ_𝔥, λ₂ = −(τ+γ+δ_𝔥), λ₃ = δ_𝔪, λ₄ = θ+δ_𝔥 and the characteristic relation is

$\begin{array}{r}{\lambda }^2+\left(\varrho +\gamma +{\delta }_{𝔥}+{\delta }_{𝔪}\right)\lambda +{\delta }_{𝔪}\left(\tau +\gamma +{\delta }_{𝔥}\right)\left(\varrho +\gamma +{\delta }_{𝔥}\right)\\ \left[1-{\mathcal{R}}_0\right]=0\end{array}$

When ${\mathcal{R}}_0<1,$ it is obvious that λ₅ < 1 and λ₆ < 1, all the Eigen values satisfy the condition $|arg({\lambda }_i)|>\frac{\eta \pi }{2},i=1,2,\dots ,6.$ the without contagious equilibrium ${\mathcal{C}}_{\mathcal{f}}$ is locally asymptotically stable.

Theorem 4.4.

If ${\mathcal{R}}_0>1$, the equilibrium point ${\mathcal{C}}_{\mathcal{p}}$ is locally stable, then system (Equation 6) has ubiquitous contagion.

Proof

Jacobian matrix evaluated in static equilibrium:

$\begin{array}{ll}det\left(\Im \left({\mathcal{C}}_{\mathcal{p}}\right)-\lambda 𝔗\right)=0& (12)\end{array}$

We obtain the Eigen values are λ₁ = (θ+δ_𝔥), ${\lambda }_2=\left({\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}^{**}+{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}\right)+{\delta }_{𝔪}\right)$,

${\lambda }_3={\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){𝔗}_{𝔪}^{**}+{\delta }_{𝔥}$ and the characteristic relation

$\begin{array}{ll}{\lambda }^3+a_1{\lambda }^2+a_2\lambda +a_3=0& (13)\end{array}$

Where

$\begin{array}{l}{a}_1=\left(\tau +\varrho +2\left(\gamma +{\delta }_{𝔥}\right)\right)\\ a_2=\left(\varrho +\gamma +{\delta }_{𝔥}+{\delta }_{𝔪}\right)\left[\left(\tau +\gamma +{\delta }_{𝔥}\right)-{{\alpha }_{𝔪}}^2\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){\chi }_{𝔪}{Ȿ}_{𝔥}^{**}{Ȿ}_{𝔪}^{**}\right]\\ a_3=\left(\tau +\gamma +{\delta }_{𝔥}\right)\left[\left(\varrho +\gamma +{\delta }_{𝔥}\right){\delta }_{𝔪}-{{\alpha }_{𝔪}}^2\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){\chi }_{𝔪}{Ȿ}_{𝔥}^{**}{Ȿ}_{𝔪}^{**}\right]\end{array}$

By using Routh-Hurwitz Criteria (22, 23), if the following provisions are handling

$\begin{array}{l}{a}_1>0,a_2>0,a_3>0\text{ and }{a}_1{a}_2-a_3>0.\end{array}$

Then ${\mathcal{C}}_{\mathcal{p}}$ is approximately stable locally. The evidence is conclusive.

4.5 Global stability

Theorem 4.5.

If ${\mathcal{R}}_0<1,$ the point of without contagious equilibrium ${\mathcal{C}}_{\mathcal{f}}$ is global stability on Φ.

Proof

Create a Lyapunov function 𝕍₁(ȶ),

$\begin{array}{ll}{𝕍}_1\left(ȶ\right)=\left({Ȿ}_{𝔥}-{Ȿ}_{𝔥}^{*}ln{Ȿ}_{𝔥}\right)+{𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}+{𝔎}_{𝔥}+\left({Ȿ}_{𝔪}-{Ȿ}_{𝔪}^{*}ln{Ȿ}_{𝔪}\right)+{𝔗}_{𝔪}& (14)\end{array}$

Calculating the fractional order derivatives of 𝕍₁(ȶ) in the solution direction of Equation 6, from Lemma 2, we obtain

$\begin{array}{l}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_1\left(ȶ\right)\le \left(1-\frac{{Ȿ}_{𝔥}^{*}}{{Ȿ}_{𝔥}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔥}\\ +\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{Ȿ}_{𝔥}}+\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{\mathcal{A}}_{𝔥}}+\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔎}_{𝔥}\\ +\left(1-\frac{{Ȿ}_{𝔪}^{*}}{{Ȿ}_{𝔪}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔪}+\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{𝔪}\\ \le \left(1-\frac{{Ȿ}_{𝔥}^{*}}{{Ȿ}_{𝔥}}\right)\left({\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\right)\\ +\left({\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\right)\\ +\left({\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\right)+(\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\\ -\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}})\\ +\left(1-\frac{{Ȿ}_{𝔪}^{*}}{{Ȿ}_{𝔪}}\right)\left({\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\right)\\ +\left({\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}\right)\\ \le {\pi }_{𝔥}\left(1-\frac{{Ȿ}_{𝔥}^{*}}{{Ȿ}_{𝔥}}\right)+{\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}^{*}{𝔗}_{𝔪}+{\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}^{*}{𝔗}_{𝔪}\\ -{\delta }_{𝔥}{Ȿ}_{𝔥}+{\delta }_{𝔥}{Ȿ}_{𝔥}^{*}+\theta {𝔎}_{𝔥}\left(1-\frac{{Ȿ}_{𝔥}^{*}}{{Ȿ}_{𝔥}}\right)-{\delta }_{𝔥}{𝔗}_{{Ȿ}_{𝔥}}\\ -{\delta }_{𝔥}{𝔗}_{{\mathcal{A}}_{𝔥}}-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+{\pi }_{𝔪}\left(1-\frac{{Ȿ}_{𝔪}^{*}}{{Ȿ}_{𝔪}}\right)+{\alpha }_{𝔪}{\chi }_{𝔪}\\ \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}^{*}-{\delta }_{𝔪}\left(1-\frac{{Ȿ}_{𝔪}^{*}}{{Ȿ}_{𝔪}}\right)-{\delta }_{𝔪}{𝔗}_{𝔪}\end{array}$

Substituting the reaction of without contagious free ${\mathcal{C}}_{\mathcal{f}}=\left(\frac{{\pi }_{𝔥}}{{\delta }_{𝔥}},0,0,\frac{{\pi }_{𝔪}}{{\delta }_{𝔪}},0\right)$, we obtain:

$\begin{array}{ll}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_1\left(ȶ\right)\le \left({\pi }_{𝔥}+{\delta }_{𝔥}\right)\left(2-\frac{{Ȿ}_{𝔥}^{*}}{{Ȿ}_{𝔥}}-\frac{{Ȿ}_{𝔥}}{{Ȿ}_{𝔥}^{*}}\right)+\left({\pi }_{𝔪}-{\delta }_{𝔪}\right)\left(2-\frac{{Ȿ}_{𝔪}^{*}}{{Ȿ}_{𝔪}}-\frac{{Ȿ}_{𝔪}}{{Ȿ}_{𝔪}^{*}}\right)& (15)\end{array}$

It is clear that each term in Equation 15 must be negative. We have $\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_1\left(ȶ\right)\le 0$ due to LaSalle's invariance principle (24), the function $\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_1\left(ȶ\right)$ is required to be negative finite.

The maximally invariant sets ${Ȿ}_{𝔥}={Ȿ}_{𝔥}^{*}$, ${Ȿ}_{𝔪}={Ȿ}_{𝔪}^{*}$ which is singleton ${\mathcal{C}}_{\mathcal{f}}=\left({Ȿ}_{𝔥}^{*},{𝔗}_{{Ȿ}_{𝔥}}^{*},{𝔗}_{{\mathcal{A}}_{𝔥}}^{*},{𝔎}_{𝔥}^{*},{Ȿ}_{𝔪}^{*},{𝔗}_{𝔪}^{*}\right)$ contains the limit set for each solution. This demonstrates ${\mathcal{C}}_{\mathcal{f}}$ is globally asymptotically stable on Φ.

Theorem 4.6.

When ${\mathcal{R}}_0>1,$ the positive contagious equalization level of system (Equation 6) arises and is globally stable on Φ.

Proof

Let's create a lyapunov function of the following form

$\begin{array}{l}{𝕍}_2\left(ȶ\right)=\left({Ȿ}_{𝔥}-{Ȿ}_{𝔥}^{**}ln{Ȿ}_{𝔥}\right)+\left({𝔗}_{{Ȿ}_{𝔥}}-{𝔗}_{{Ȿ}_{𝔥}}^{**}ln{𝔗}_{{Ȿ}_{𝔥}}\right)\\ +\left({𝔗}_{{\mathcal{A}}_{𝔥}}-{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}ln{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\\ +\left({𝔎}_{𝔥}-{𝔎}_{𝔥}^{**}ln{𝔎}_{𝔥}\right)+\left({Ȿ}_{𝔪}-{Ȿ}_{𝔪}^{**}ln{Ȿ}_{𝔪}\right)+\left({𝔗}_{𝔪}-{𝔗}_{𝔪}^{**}ln{𝔗}_{𝔪}\right)\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_2\left(ȶ\right)\le \left(1-\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔥}\\ +\left(1-\frac{{𝔗}_{{Ȿ}_{𝔥}}^{**}}{{𝔗}_{{Ȿ}_{𝔥}}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{Ȿ}_{𝔥}}+\left(1-\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{\mathcal{A}}_{𝔥}}+\\ \left(1-\frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔎}_{𝔥}\\ +\left(1-\frac{{Ȿ}_{𝔪}^{**}}{{Ȿ}_{𝔪}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔪}+\left(1-\frac{{𝔗}_{𝔪}^{**}}{{𝔗}_{𝔪}}\right)\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{𝔪}\\ \le \left(1-\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}\right)\left({\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\right)\\ +\left(1-\frac{{𝔗}_{{Ȿ}_{𝔥}}^{**}}{{𝔗}_{{Ȿ}_{𝔥}}}\right)\left({\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{{Ȿ}_{𝔥}}\right)\\ +\left(1-\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)\left({\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\right)+\left(1-\frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}\right)\\ \left(\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\right)\\ +\left(1-\frac{{Ȿ}_{𝔪}^{**}}{{Ȿ}_{𝔪}}\right)\left({\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\right)+\left(1-\frac{{𝔗}_{𝔪}^{**}}{{𝔗}_{𝔪}}\right)\\ \left({\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}\right)\\ \le {\pi }_{𝔥}\left(1-\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}\right)+{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}^{**}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}\left(1-\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}\right)\\ +\theta {𝔎}_{𝔥}\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}+{𝔗}_{{Ȿ}_{𝔥}}^{**}\left(\tau +\gamma \right)-{\delta }_{𝔥}{𝔗}_{{Ȿ}_{𝔥}}\left(1-\frac{{𝔗}_{{Ȿ}_{𝔥}}^{**}}{{𝔗}_{{Ȿ}_{𝔥}}}\right)-{\alpha }_{𝔪}{Ȿ}_{𝔥}\\ \left({\chi }_{{Ȿ}_{𝔥}}\frac{{𝔗}_{{Ȿ}_{𝔥}}^{**}}{{𝔗}_{{Ȿ}_{𝔥}}}+{\chi }_{{\mathcal{A}}_{𝔥}}\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right){𝔗}_{𝔪}+\left(\varrho +\gamma \right){𝔗}_{{\mathcal{A}}_{𝔥}}^{**}\\ -{\delta }_{𝔥}{𝔗}_{{\mathcal{A}}_{𝔥}}\left(1-\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right)-\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\\ \frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}+\theta {𝔎}_{𝔥}^{**}-{\delta }_{𝔥}{𝔎}_{𝔥}\left(1-\frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}\right)+\left(\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\right)\frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}\\ +{\pi }_{𝔪}\left(1-\frac{{Ȿ}_{𝔪}^{**}}{{Ȿ}_{𝔪}}\right)-{\delta }_{𝔪}\left(1-\frac{{Ȿ}_{𝔪}^{**}}{{Ȿ}_{𝔪}}\right)-{\delta }_{𝔪}\left(1-\frac{{𝔗}_{𝔪}^{**}}{{𝔗}_{𝔪}}\right)\\ +{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\left({Ȿ}_{𝔪}^{**}+{Ȿ}_{𝔪}\frac{{𝔗}_{𝔪}^{**}}{{𝔗}_{𝔪}}\right)\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_2\left(ȶ\right)\le \left({\pi }_{𝔥}-{\delta }_{𝔥}+\theta {𝔎}_{𝔥}\right)\\ \left(2-\frac{{Ȿ}_{𝔥}^{**}}{{Ȿ}_{𝔥}}-\frac{{Ȿ}_{𝔥}}{{Ȿ}_{𝔥}^{**}}\right)+\left(\tau +\gamma -{\delta }_{𝔥}\right)\left(2-\frac{{𝔗}_{{Ȿ}_{𝔥}}^{**}}{{𝔗}_{{Ȿ}_{𝔥}}}-\frac{{𝔗}_{{Ȿ}_{𝔥}}}{{𝔗}_{{Ȿ}_{𝔥}}^{**}}\right)\\ +\left(\varrho +\gamma -{\delta }_{𝔥}\right)\\ \left(2-\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}{{𝔗}_{{\mathcal{A}}_{𝔥}}}-\frac{{𝔗}_{{\mathcal{A}}_{𝔥}}}{{𝔗}_{{\mathcal{A}}_{𝔥}}^{**}}\right)+\left(\gamma +{\delta }_{𝔥}-\theta \right)\left(2-\frac{{𝔎}_{𝔥}^{**}}{{𝔎}_{𝔥}}-\frac{{𝔎}_{𝔥}}{{𝔎}_{𝔥}^{**}}\right)\\ +\left({\pi }_{𝔪}-{\delta }_{𝔪}\right)\left(2-\frac{{Ȿ}_{𝔪}^{**}}{{Ȿ}_{𝔪}}\frac{{Ȿ}_{𝔪}}{{Ȿ}_{𝔪}^{**}}\right)+{\delta }_{𝔪}\left(2-\frac{{𝔗}_{𝔪}^{**}}{{𝔗}_{𝔪}}-\frac{{𝔗}_{𝔪}}{{𝔗}_{𝔪}^{**}}\right)& (16)\end{array}$

Hence, the condition in Equation 16 ensures

$\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝕍}_2\left(ȶ\right)\le 0$ for all $\left({Ȿ}_{𝔥}^{**},{𝔗}_{{Ȿ}_{𝔥}}^{**},{𝔗}_{{\mathcal{A}}_{𝔥}}^{**},{𝔎}_{𝔥}^{**},S_{𝔪}^{**},{𝔗}_{𝔪}^{**}\right)\in \Phi$ and strict the quality holds for $S_{𝔥}={Ȿ}_{𝔥}={Ȿ}_{𝔥}^{**},{𝔗}_{S_{𝔥}}={𝔗}_{S_{𝔥}}^{**},$ ${𝔗}_{{\mathcal{A}}_{𝔥}}={𝔗}_{{\mathcal{A}}_{𝔥}}^{**},{𝔎}_{𝔥}={𝔎}_{𝔥}^{**},{Ȿ}_{𝔪}={Ȿ}_{𝔪}^{**}$ and ${𝔗}_{𝔪}={𝔗}_{𝔪}^{**}.$ therefore the equilibrium point ${\mathcal{C}}_{\mathcal{p}}$ becomes globally stable on Φ.

5 Optimum control approach

In this portion, we will discuss how to optimize the problem and analyze the performance of the control function. Consolidation of optimal controlling problem a dynamics of control system can be described as system (Equation 6).

$\begin{array}{l}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔥}={\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)S_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{S}_{𝔥}+\theta {𝔎}_{𝔥}-{𝕌}_1\left(ȶ\right)S_{𝔥}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{Ȿ}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{S}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{S_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{{\mathcal{A}}_{𝔥}}={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{Ȿ}_{𝔪}={\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}-{𝕌}_2\left(ȶ\right){Ȿ}_{𝔪}\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{𝔗}_{𝔪}={\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}-{𝕌}_2\left(ȶ\right){𝔗}_{𝔪}& (17)\end{array}$

Where

𝕌₁− Self-precaution (long sleeved pants and shorts, increase immune system, consultation at

the neatest health care) minimizes the susceptible individuals.

𝕌₂− Use of chemical insecticide sprays destroying the susceptible and infected mosquito cases

The optimal solution being minimized could be expressed as:

$\begin{array}{ll}ℂ\left({𝕌}_1,{𝕌}_2\right)={\displaystyle {\int }_0^{{ȶ}_f}}\left(\mathcal{a}{Ȿ}_{𝔥}+\mathcal{b}{𝔗}_{{Ȿ}_{𝔥}}+\mathcal{c}{𝔗}_{{\mathcal{A}}_{𝔥}}+\mathcal{d}{Ȿ}_{𝔪}+\mathcal{e}{𝔗}_{𝔪}+\mathcal{f}{{𝕌}_1}^2+\mathcal{g}{{𝕌}_2}^2\right)\mathcal{d}ȶ& (18)\end{array}$

To reduce the cost of two controls 𝕌₁ and 𝕌₂ the objective is reduced ${Ȿ}_{𝔥},{𝔗}_{S_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}}$ and S_𝔪, 𝔗_𝔪.

Therefore, we need to obtain optimal controls ${𝕌}_1^{*}$ and ${𝕌}_2^{*}$

$\begin{array}{ll}ℂ\left({𝕌}_1^{*},{𝕌}_2^{*}\right)={min}_{{}_{{𝕌}_1,{𝕌}_2}}\left\{ℂ\left({𝕌}_1,{𝕌}_2\right)|{𝕌}_1,{𝕌}_2\in \Phi \right\}& (19)\end{array}$

A set of constraints $\Phi =\left\{\left({𝕌}_1,{𝕌}_2\right)|{𝕌}_{\mathcal{i}}:\left[0,{ȶ}_f\right]\to \left[0,\infty \right)lebesquequantifiable\mathcal{i}=1,2\right\}$.

The expense of minimizing ${Ȿ}_{𝔥},{𝔗}_{{Ȿ}_{𝔥}},{𝔗}_{{\mathcal{A}}_{𝔥}}$, Ȿ_𝔪 and 𝔗_𝔪 is represented by the term $\mathcal{a}{Ȿ}_{𝔥},\mathcal{b}{𝔗}_{{Ȿ}_{𝔥}},\mathcal{c}{𝔗}_{{\mathcal{A}}_{𝔥}},\mathcal{d}{Ȿ}_{𝔪}$ and $\mathcal{e}{𝔗}_{𝔪}$ respectively. Likewise, $f{𝕌}_1^2,\mathcal{g}{𝕌}_2^2$ represents the cost for controls 𝕌₁, 𝕌₂. The most prevalent PMP can be used to find the adequacy condition required for the control system to be satisfied. Equations 17, 19 can be transformed into the following point-wise Hamiltonian ℍ for (𝕌₁, 𝕌₂) regression problem using the aforesaid principle.

$\begin{array}{l}ℍ=\left\{\mathcal{a}{Ȿ}_{𝔥}+\mathcal{b}{𝔗}_{S_{𝔥}}+\mathcal{c}{𝔗}_{{\mathcal{A}}_{𝔥}}+\mathcal{d}{S}_{𝔪}+\mathcal{e}{𝔗}_{𝔪}+f{{𝕌}_1}^2+\mathcal{g}{{𝕌}_2}^2\right\}\\ +{\lambda }_{{Ȿ}_{𝔥}}\left\{{\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{S_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}-{𝕌}_1\left(ȶ\right){Ȿ}_{𝔥}\right\}\\ +{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left\{{\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{S}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{S_{𝔥}}\right\}\\ +{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left\{{\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\right\}\\ +{\lambda }_{{Ȿ}_{𝔪}}\left\{{\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{S_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{S}_{𝔪}-{𝕌}_2\left(ȶ\right)S_{𝔪}\right\}\\ +{\lambda }_{{𝔗}_{𝔪}}\left\{{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-{\delta }_{𝔪}{𝔗}_{𝔪}-{𝕌}_2\left(ȶ\right){𝔗}_{𝔪}\right\}& (20)\end{array}$

Where ${\lambda }_{{Ȿ}_{𝔥}},{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}},{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}},{\lambda }_{{\mathcal{S}}_{\mathcal{m}}}$ and λ_𝔗𝔪 are the ad-joint variable or co-state variable.

$\begin{array}{l}\frac{d{\lambda }_{{Ȿ}_{𝔥}}}{dȶ}=\frac{\partial ℍ}{\partial {Ȿ}_{𝔥}}=\mathcal{a}+{\lambda }_{{Ȿ}_{𝔥}}\left\{-{\alpha }_{𝔪}\left({\chi }_{S_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){𝔗}_{𝔪}-{\delta }_{𝔥}\right\}\\ -\left({\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}+{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right){\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){𝔗}_{𝔪}\\ \frac{d{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}}{dȶ}=\frac{\partial ℍ}{\partial {𝔗}_{{Ȿ}_{𝔥}}}=\mathcal{b}-{\lambda }_{{Ȿ}_{𝔥}}{𝕌}_1\left(ȶ\right)+{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left\{-\left(\tau +\gamma +{\delta }_{𝔥}\right)\right\}\\ +{\lambda }_{{𝔗}_{𝔪}}{\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}\\ \frac{d{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}}{dȶ}=\frac{\partial ℍ}{\partial {𝔗}_{{\mathcal{A}}_{𝔥}}}=\mathcal{c}-{\lambda }_{{Ȿ}_{𝔥}}{𝕌}_1\left(ȶ\right)\\ +{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left\{-\left(\varrho +\gamma +{\delta }_{𝔥}\right)\right\}+{\lambda }_{{𝔗}_{𝔪}}{\alpha }_{𝔪}{\chi }_{𝔪}{Ȿ}_{𝔪}\\ \frac{d{\lambda }_{{Ȿ}_{𝔪}}}{dȶ}=\frac{\partial ℍ}{\partial {Ȿ}_{𝔪}}\\ =\mathcal{d}+{\lambda }_{{Ȿ}_{𝔪}}\left\{-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-{\delta }_{𝔪}-{𝕌}_2\left(ȶ\right)\right\}\\ +{\lambda }_{{𝔗}_{𝔪}}\left({\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\right)\\ \frac{d{\lambda }_{{𝔗}_{𝔪}}}{dȶ}=\frac{\partial ℍ}{\partial {𝔗}_{𝔪}}=\mathcal{e}+{\lambda }_{{𝔗}_{𝔪}}\{-{\delta }_{𝔪}-{𝕌}_2\left(ȶ\right)\}\\ -\left({\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}+{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\right){\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)S_{𝔥}+{\lambda }_{S_{𝔥}}\\ \left\{-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔥}\right\}& (21)\end{array}$

The conditions for transversality are

${\lambda }_{{Ȿ}_{𝔥}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left(ȶ\right)={\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left({ȶ}_f\right)={\lambda }_{{Ȿ}_{𝔪}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{𝔪}}\left({ȶ}_f\right)=0.$

For $0<{𝕌}_{\mathcal{i}}<1,\mathcal{i}=1,2.$ From the interior of controls, we have:

$\begin{array}{l}\frac{\partial ℍ}{\partial {𝕌}_1}=2f{𝕌}_1-{\lambda }_{{Ȿ}_{𝔥}}\left({𝔗}_{S_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\\ \frac{\partial ℍ}{\partial {𝕌}_2}=2\mathcal{g}{𝕌}_2-{\lambda }_{{Ȿ}_{𝔪}}{Ȿ}_{𝔪}-{\lambda }_{{𝔗}_{𝔪}}{𝔗}_{𝔪}& (22)\end{array}$

From where:

$\begin{array}{l}{𝕌}_1=\frac{{\lambda }_{{Ȿ}_{𝔥}}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)}{2f}\\ {𝕌}_2=\frac{{\lambda }_{{Ȿ}_{𝔪}}{Ȿ}_{𝔪}+{\lambda }_{{𝔗}_{𝔪}}{𝔗}_{𝔪}}{2\mathcal{g}}& (23)\end{array}$

5.1 Utilization of optimal solutions

Theorem 5.1. $\left({𝕌}_1^{*},{𝕌}_2^{*}\right)$ is a control factor can reduce over 𝕌 provided by

$\begin{array}{l}{𝕌}_1^{*}=max\left\{0,min\left\{1,{\frac{1}{2f}\lambda }_{{Ȿ}_{𝔥}}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)\right\}\right\}\\ {𝕌}_2^{*}=max\left\{0,min\left\{1,\frac{1}{2\mathcal{g}}{\lambda }_{{Ȿ}_{𝔪}}{Ȿ}_{𝔪}+{\lambda }_{{𝔗}_{𝔪}}{𝔗}_{𝔪}\right\}\right\}& (24)\end{array}$

Where ${\lambda }_{{Ȿ}_{𝔥}},{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}},{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}},{\lambda }_{{Ȿ}_{𝔪}}$ and λ_𝔗𝔪 are co-state variable that satisfy the condition (Equations 17–24) in addition, the transversality characteristic that follows

$\begin{array}{l}{\lambda }_{{Ȿ}_{𝔥}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left({ȶ}_f\right)={\lambda }_{{Ȿ}_{𝔪}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{𝔪}}\left({ȶ}_f\right)=0.\end{array}$

$\begin{array}{l}{𝕌}_1^{*}=\{\begin{array}{c}0if{𝕌}_1\le 0,\\ {𝕌}_{1}if0<{𝕌}_1<1,\\ 1if{𝕌}_1\ge 0.\end{array}\end{array}$

And

$\begin{array}{ll}{𝕌}_2^{*}=\{\begin{array}{c}0if{𝕌}_2\le 0,\\ {𝕌}_{2}if0<{𝕌}_2<1,\\ 1if{𝕌}_2\ge 0.\end{array}& (25)\end{array}$

Proof

To demonstrate the survival of optimal control solutions, the configuration of the Lipschitz criterion of the system and the convexity of the integral in Equation 21 are related and state variable that constrains 𝕌₁ and 𝕌₂ to the boundary of the state solution. So we employ PMP and get the following:

$\begin{array}{l}\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{\lambda }_{{Ȿ}_{𝔥}}\left(ȶ\right)=\frac{\partial ℍ}{\partial {Ȿ}_{𝔥}};\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left(ȶ\right)=\frac{\partial ℍ}{\partial {𝔗}_{{Ȿ}_{𝔥}}};\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left(ȶ\right)=\frac{\partial ℍ}{\partial {𝔗}_{{\mathcal{A}}_{𝔥}}};\\ \begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{\lambda }_{{Ȿ}_{𝔪}}\left(ȶ\right)=\frac{\partial ℍ}{\partial {Ȿ}_{𝔪}};\begin{array}{r}\mathcal{A}ℬ\mathcal{C}\\ 0\end{array}{𝔇}_{ȶ}^{\eta }{\lambda }_{{𝔗}_{𝔪}}\left(ȶ\right)=\frac{\partial ℍ}{\partial {𝔗}_{𝔪}};& (26)\end{array}$

with,

${\lambda }_{{Ȿ}_{𝔥}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{{Ȿ}_{𝔥}}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{{\mathcal{A}}_{𝔥}}}\left({ȶ}_f\right)={\lambda }_{{Ȿ}_{𝔪}}\left({ȶ}_f\right)={\lambda }_{{𝔗}_{𝔪}}\left({ȶ}_f\right)=0.$

The Hamilton can be differentiated with regard to achieve the conditional optimum:

$\begin{array}{ll}\frac{\partial ℍ}{\partial {𝕌}_1}=0,\frac{\partial ℍ}{\partial {𝕌}_2}=0.& (27)\end{array}$

The ad-joint system (Equations 20, 21) derived from Equation 17, the optimum system (Equation 23) is accessible from Equation 24. The optimal method is the constrained system (Equation 17) and its initial state is ad-joint the system includes (Equation 20), and condition for intersection.

6 Adams-Bash forth method

Here, we formulate the system of Equation 6 a recently invented numerical approach, the Adams-Bash forth method (24). The framework (Equation 6) can be used to test the essential theorem from fractional calculus,

$\begin{array}{l}{Ȿ}_{𝔥}\left(ȶ\right)={Ȿ}_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_1(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),{𝔗}_{S_{𝔥}}\left(\varpi \right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(\varpi \right),\\ {𝔎}_{𝔥}\left(\varpi \right),{Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (28)\end{array}$

$\begin{array}{l}{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right)={𝔗}_{{Ȿ}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_2(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),{𝔗}_{{Ȿ}_{𝔥}}\left(\varpi \right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(\varpi \right),{𝔎}_{𝔥}\\ \left(\varpi \right),{Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (29)\end{array}$

$\begin{array}{l}{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right)={𝔗}_{{\mathcal{A}}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_3(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),\\ T_{{Ȿ}_h}\left(\varpi \right),T_{{\mathcal{A}}_h}\left(\varpi \right),K_h\left(\varpi \right),\\ {Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (30)\end{array}$

$\begin{array}{l}{𝔎}_{𝔥}\left(ȶ\right)={𝔎}_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_4(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_4(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),{𝔗}_{S_{𝔥}}\left(\varpi \right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(\varpi \right),\\ {𝔎}_{𝔥}\left(\varpi \right),{Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (31)\end{array}$

$\begin{array}{l}{Ȿ}_{𝔪}\left(ȶ\right)={Ȿ}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_5(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),\\ {𝔗}_{S_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_5(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),{𝔗}_{S_{𝔥}}\left(\varpi \right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(\varpi \right),\\ {𝔎}_{𝔥}\left(\varpi \right),{Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (32)\end{array}$

$\begin{array}{l}{𝔗}_{𝔪}\left(ȶ\right)={𝔗}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{A}\mathcal{B}\mathcal{C}\left(\eta \right)}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right))+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{ȶ}}{𝕂}_6(\varpi ,{Ȿ}_{𝔥}\left(\varpi \right),{𝔗}_{{Ȿ}_{𝔥}}\left(\varpi \right),\\ {𝔗}_{{\mathcal{A}}_{𝔥}}\left(\varpi \right),{𝔎}_{𝔥}\left(\varpi \right),{Ȿ}_{𝔪}\left(\varpi \right),{𝔗}_{𝔪}\left(\varpi \right)){\left(ȶ-\varpi \right)}^{\eta -1}\mathcal{d}\varpi & (33)\end{array}$

Where,

$\begin{array}{ll}\{\begin{array}{l}{𝕂}_1\left(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ ={\pi }_{𝔥}-{\alpha }_{𝔪}\left({\chi }_{{Ȿ}_{𝔥}}+{\chi }_{{\mathcal{A}}_{𝔥}}\right)S_{𝔥}{𝔗}_{𝔪}-{\delta }_{𝔥}{Ȿ}_{𝔥}+\theta {𝔎}_{𝔥}\\ {𝕂}_2\left(ȶ,S_{𝔥}\left(ȶ\right),{𝔗}_{S_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),S_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ ={\alpha }_{𝔪}{\chi }_{{Ȿ}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\tau +\gamma +{\delta }_{𝔥}\right){𝔗}_{S_{𝔥}}\\ {𝕂}_3\left(ȶ,S_{𝔥}\left(ȶ\right),{𝔗}_{S_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),S_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ ={\alpha }_{𝔪}{\chi }_{{\mathcal{A}}_{𝔥}}{Ȿ}_{𝔥}{𝔗}_{𝔪}-\left(\varrho +\gamma +{\delta }_{𝔥}\right){𝔗}_{{\mathcal{A}}_{𝔥}}\\ {𝕂}_4\left(ȶ,S_{𝔥}\left(ȶ\right),{𝔗}_{S_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),S_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ =\gamma \left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)-\theta {𝔎}_{𝔥}-{\delta }_{𝔥}{𝔎}_{𝔥}+\tau {𝔗}_{{Ȿ}_{𝔥}}+\varrho {𝔗}_{{\mathcal{A}}_{𝔥}}\\ {𝕂}_5\left(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ ={\pi }_{𝔪}-{\alpha }_{𝔪}{\chi }_{𝔪}\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right){Ȿ}_{𝔪}-\Omega {Ȿ}_{𝔪}-{\delta }_{𝔪}{Ȿ}_{𝔪}\\ {𝕂}_6\left(ȶ,S_{𝔥}\left(ȶ\right),{𝔗}_{S_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),S_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)\right)\\ ={\alpha }_{𝔪}{\chi }_{𝔪}{\left({𝔗}_{{Ȿ}_{𝔥}}+{𝔗}_{{\mathcal{A}}_{𝔥}}\right)Ȿ}_{\mathcal{m}}-\left(\Omega +{\delta }_{𝔪}\right){𝔗}_{𝔪}\end{array}& (34)\end{array}$

The following structure is obtained at time t_{𝔫+ 1},

$\begin{array}{l}{Ȿ}_{𝔥}\left({ȶ}_{𝔫+1}\right)={Ȿ}_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_1\\ \left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (35)\end{array}$

$\begin{array}{l}{𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫+1}\right)={𝔗}_{{Ȿ}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_2\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (36)\end{array}$

$\begin{array}{l}{𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫+1}\right)={𝔗}_{{\mathcal{A}}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_3\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (37)\end{array}$

$\begin{array}{l}{𝔎}_{𝔥}\left({ȶ}_{𝔫+1}\right)={𝔎}_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_4\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_4(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (38)\end{array}$

$\begin{array}{l}{Ȿ}_{𝔪}\left({ȶ}_{𝔫+1}\right)={Ȿ}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_5\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_5(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (39)\end{array}$

$\begin{array}{l}{𝔗}_{𝔪}\left({ȶ}_{𝔫+1}\right)={𝔗}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_6\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (40)\end{array}$

While, at t_𝔫 we have

$\begin{array}{l}{Ȿ}_{𝔥}\left({ȶ}_{𝔫}\right)=S_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_1({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},\\ {𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}})+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (41)\end{array}$

$\begin{array}{l}{𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫}\right)={𝔗}_{{Ȿ}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_2({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},\\ {𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}})\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (42)\end{array}$

$\begin{array}{l}{𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫}\right)={𝔗}_{{\mathcal{A}}_{𝔥}}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}{𝕂}_3({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},\\ {𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}})\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ0\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (43)\end{array}$

$\begin{array}{l}{𝔎}_{𝔥}\left({ȶ}_{𝔫}\right)={𝔎}_{𝔥}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_4\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_4(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (44)\end{array}$

$\begin{array}{l}{Ȿ}_{𝔪}\left({ȶ}_{𝔫}\right)={Ȿ}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_5\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_5(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (45)\end{array}$

$\begin{array}{l}{𝔗}_{𝔪}\left({ȶ}_{𝔫}\right)={𝔗}_{𝔪}\left(0\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ {𝕂}_6\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},{Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (46)\end{array}$

By subtracting Ȿ_𝔥(ȶ_𝔫) from Ȿ_𝔥(ȶ_𝔫+1), 𝔗_Ȿ𝔥(ȶ_𝔫) from 𝔗_Ȿ𝔥(ȶ_𝔫+1), ${𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫}\right)$ from ${𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫+1}\right)$, 𝔎_𝔥(ȶ_𝔫) from 𝔎_𝔥(ȶ_𝔫+1), Ȿ_𝔪(ȶ_𝔫) from S_𝔪(ȶ_𝔫+1) and 𝔗_𝔪(ȶ_𝔫) from 𝔗_𝔪(ȶ_𝔫+1), we get the following

$\begin{array}{l}{Ȿ}_{𝔥}\left({ȶ}_{𝔫+1}\right)=S_{𝔥}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_1\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_1\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (47)\end{array}$

$\begin{array}{l}{𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫+1}\right)={𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_2\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_2\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (48)\end{array}$

$\begin{array}{l}{𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫+1}\right)={𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_3\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_3\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (49)\end{array}$

$\begin{array}{l}{𝔎}_{𝔥}\left({ȶ}_{𝔫+1}\right)={𝔎}_{𝔥}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_4\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)\\ -{𝕂}_4\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_4\left(ȶ,{Ȿ}_{𝔥}(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_4(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (50)\end{array}$

$\begin{array}{l}{Ȿ}_{𝔪}\left({ȶ}_{𝔫+1}\right)=S_{𝔪}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_5\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)\\ -{𝕂}_5\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (51)\end{array}$

$\begin{array}{l}{𝔗}_{𝔪}\left({ȶ}_{𝔫+1}\right)={𝔗}_{𝔪}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_6\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_6\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]+\\ \frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ-\\ +\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),{𝔎}_{𝔥}\left(ȶ\right),\\ {Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫}-t\right)}^{\eta -1}\mathcal{d}ȶ& (52)\end{array}$

The Equations 47–52 become

$\begin{array}{l}{Ȿ}_{𝔥}\left({ȶ}_{𝔫+1}\right)=S_{𝔥}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_1({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},\\ {Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}(-\\ {𝕂}_1\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +{𝔸}_{\eta ,1}^1-{𝔸}_{\eta ,2}^1& (53)\end{array}$

$\begin{array}{l}{𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫+1}\right)={𝔗}_{{Ȿ}_{𝔥}}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_2({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},\\ {Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}(-\\ {𝕂}_2({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},\\ {Ȿ}_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}(\end{array}\right]\\ +{𝔸}_{\eta ,1}^2-{𝔸}_{\eta ,2}^2& (54)\end{array}$

$\begin{array}{l}{𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫+1}\right)+{𝔗}_{{\mathcal{A}}_{𝔥}}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_3\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_3\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +{𝔸}_{\eta ,1}^3-{𝔸}_{\eta ,2}^3& (55)\end{array}$

$\begin{array}{l}{𝔎}_{𝔥}\left({ȶ}_{𝔫+1}\right)={𝔎}_{𝔥}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_4\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_4\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +{𝔸}_{\eta ,1}^4-{𝔸}_{\eta ,2}^4& (56)\end{array}$

$\begin{array}{l}{Ȿ}_{𝔪}\left({ȶ}_{𝔫+1}\right)=S_{𝔪}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_5\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_5\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +{𝔸}_{\eta ,1}^5-{𝔸}_{\eta ,2}^5& (57)\end{array}$

$\begin{array}{l}{𝔗}_{𝔪}\left({ȶ}_{𝔫+1}\right)={𝔗}_{𝔪}\left({ȶ}_{𝔫}\right)+\frac{1-\eta }{\mathcal{ABC}\left(\eta \right)}\\ \left[\begin{array}{c}{𝕂}_6\left({ȶ}_{𝔫},{Ȿ}_{{𝔥}_{𝔫}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫}}},{𝔎}_{{𝔥}_{𝔫}},{Ȿ}_{{𝔪}_{𝔫}},{𝔗}_{{𝔪}_{𝔫}}\right)-\\ {𝕂}_6\left({ȶ}_{𝔫-1},{Ȿ}_{{𝔥}_{𝔫-1}},{𝔗}_{{Ȿ}_{{𝔥}_{𝔫-1}}},{𝔗}_{{\mathcal{A}}_{{𝔥}_{𝔫-1}}},{𝔎}_{{𝔥}_{𝔫-1}},S_{{𝔪}_{𝔫-1}},{𝔗}_{{𝔪}_{𝔫-1}}\right)\end{array}\right]\\ +{𝔸}_{\eta ,1}^6-{𝔸}_{\eta ,2}^6& (58)\end{array}$

Where

$\begin{array}{c}{𝔸}_{\eta ,1}^1=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (59)\end{array}$

$\begin{array}{c}{𝔸}_{\eta ,1}^2=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_2(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (60)\end{array}$

$\begin{array}{c}{𝔸}_{\eta ,1}^3=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_3(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (61)\end{array}$

$\begin{array}{c}{𝔸}_{\eta ,1}^4=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_4(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (62)\end{array}$

$\begin{array}{c}{𝔸}_{\eta ,1}^5=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_5(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (63)\end{array}$

$\begin{array}{c}{𝔸}_{\eta ,1}^6=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫+1}}}{𝕂}_6(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){\left({ȶ}_{𝔫+1}-t\right)}^{\eta -1}\mathcal{d}ȶ& (64)\end{array}$

and

$\begin{array}{c}{𝔸}_{\eta ,2}^1=\frac{\eta }{\mathcal{ABC}\left(\eta \right)\lceil \eta }{\displaystyle {\int }_0^{{ȶ}_{𝔫}}}{𝕂}_1(ȶ,{Ȿ}_{𝔥}\left(ȶ\right),{𝔗}_{{Ȿ}_{𝔥}}\left(ȶ\right),{𝔗}_{{\mathcal{A}}_{𝔥}}\left(ȶ\right),\\ {𝔎}_{𝔥}\left(ȶ\right),{Ȿ}_{𝔪}\left(ȶ\right),{𝔗}_{𝔪}\left(ȶ\right)){({�}_{}}^{}\end{array}$