Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory

Maysaa Al Qurashi; Saima Rashid; Sobia Sultana; Fahd Jarad; Abdullah M. Alsharif; Maysaa Al Qurashi; Saima Rashid; Sobia Sultana; Fahd Jarad; Abdullah M. Alsharif

doi:10.3934/math.2022697

AIMS Mathematics

2022, Volume 7, Issue 7: 12587-12619. doi: 10.3934/math.2022697

Previous Article Next Article

Research article Special Issues

Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory

1.
Department of Mathematics, King Saud University, P. O. Box 22452, Riyadh 11495, Saudi Arabia
2.
Department of Mathematics, Government College University, Faisalabad, Pakistan
3.
Department of Mathematics, Imam Mohammad Ibn Saud Islamic University, Riyadh 12211, Saudi Arabia
4.
Department of Mathematics, Cankaya University, Ankara, Turkey
5.
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan
6.
Department of Mathematics, King Abdulaziz University, Jeddah, Saudi Arabia
7.
Department of Mathematics and Statistics, College of Science, Taif University, P. O. Box 11099, Taif 21944, Saudi Arabia

Received: 26 February 2022 Revised: 14 April 2022 Accepted: 19 April 2022 Published: 28 April 2022
MSC : 46S40, 47H10, 54H25

In this research, the $\bar{\mathbf{q}}$ -homotopy analysis transform method ( $\bar{\mathbf{q}}$ -HATM) is employed to identify fractional-order Whitham–Broer–Kaup equation (WBKE) solutions. The WBKE is extensively employed to examine tsunami waves. With the aid of Caputo and Atangana-Baleanu fractional derivative operators, to obtain the analytical findings of WBKE, the predicted algorithm employs a combination of $\bar{\mathbf{q}}$ -HAM and the Aboodh transform. The fractional operators are applied in this work to show how important they are in generalizing the frameworks connected with kernels of singularity and non-singularity. To demonstrate the applicability of the suggested methodology, various relevant problems are solved. Graphical and tabular results are used to display and assess the findings of the suggested approach. In addition, the findings of our recommended approach were analyzed in relation to existing methods. The projected approach has fewer processing requirements and a better accuracy rate. Ultimately, the obtained results reveal that the improved strategy is both trustworthy and meticulous when it comes to assessing the influence of nonlinear systems of both integer and fractional order.

Keywords:

Citation: Maysaa Al Qurashi, Saima Rashid, Sobia Sultana, Fahd Jarad, Abdullah M. Alsharif. Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory[J]. AIMS Mathematics, 2022, 7(7): 12587-12619. doi: 10.3934/math.2022697

Related Papers:

[1]	Nadiyah Hussain Alharthi, Abdon Atangana, Badr S. Alkahtani . Numerical analysis of some partial differential equations with fractal-fractional derivative. AIMS Mathematics, 2023, 8(1): 2240-2256. doi: 10.3934/math.2023116
[2]	Abdon Atangana, Ali Akgül . Analysis of a derivative with two variable orders. AIMS Mathematics, 2022, 7(5): 7274-7293. doi: 10.3934/math.2022406
[3]	Hasib Khan, Jehad Alzabut, Anwar Shah, Sina Etemad, Shahram Rezapour, Choonkil Park . A study on the fractal-fractional tobacco smoking model. AIMS Mathematics, 2022, 7(8): 13887-13909. doi: 10.3934/math.2022767
[4]	Khaled M. Saad, Manal Alqhtani . Numerical simulation of the fractal-fractional reaction diffusion equations with general nonlinear. AIMS Mathematics, 2021, 6(4): 3788-3804. doi: 10.3934/math.2021225
[5]	Amir Ali, Abid Ullah Khan, Obaid Algahtani, Sayed Saifullah . Semi-analytical and numerical computation of fractal-fractional sine-Gordon equation with non-singular kernels. AIMS Mathematics, 2022, 7(8): 14975-14990. doi: 10.3934/math.2022820
[6]	Abdon Atangana, Seda İğret Araz . Extension of Chaplygin's existence and uniqueness method for fractal-fractional nonlinear differential equations. AIMS Mathematics, 2024, 9(3): 5763-5793. doi: 10.3934/math.2024280
[7]	Rahat Zarin, Amir Khan, Pushpendra Kumar, Usa Wannasingha Humphries . Fractional-order dynamics of Chagas-HIV epidemic model with different fractional operators. AIMS Mathematics, 2022, 7(10): 18897-18924. doi: 10.3934/math.20221041
[8]	Muhammad Farman, Ali Akgül, Sameh Askar, Thongchai Botmart, Aqeel Ahmad, Hijaz Ahmad . Modeling and analysis of fractional order Zika model. AIMS Mathematics, 2022, 7(3): 3912-3938. doi: 10.3934/math.2022216
[9]	Muhammad Farman, Ali Akgül, Kottakkaran Sooppy Nisar, Dilshad Ahmad, Aqeel Ahmad, Sarfaraz Kamangar, C Ahamed Saleel . Epidemiological analysis of fractional order COVID-19 model with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(1): 756-783. doi: 10.3934/math.2022046
[10]	Muhammad Farman, Aqeel Ahmad, Ali Akgül, Muhammad Umer Saleem, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Dynamical behavior of tumor-immune system with fractal-fractional operator. AIMS Mathematics, 2022, 7(5): 8751-8773. doi: 10.3934/math.2022489

Abstract

1. Introduction

Chemical kinetics deals with chemistry experiments and interprets them in terms of a mathematical model. The experiments are done on chemical reactions with the passage of time. The models are differential equations for the rates at which reactants are consumed and products are produced. Chemists are able to understand how chemical reactions take place at the molecular level by combining models with investigation. Molecules react in steps to lead to the overall stoichiometric reaction which is reaction mechanism for collection of reactions. The set of reactions specifies the path (or paths) that reactant molecules take to finally arrive at the product molecules. All species in the reaction appear in at least one step and the sum of the steps gives the overall reaction. The govern the rate of the reaction which leads directly to the mechanism of differential equations ^[1]. Many processes and phenomena in chemistry generally in sciences can be designated by first-order differential equations. These equations are the most important and most frequently used to describe natural laws. The following examples are discussed: the Bouguer-Lambert-Beer law in spectroscopy, time constants of sensors, chemical reaction kinetics, radioactive decay, relaxation in nuclear magnetic resonance, and the RC constant of an electrode ^[2]. The induced kinetic differential equations of a reaction network endowed with mass action type kinetics are a system of polynomial differential equations ^[3]. We review the basic ideas of fractional differential equations and their applications on non-linear biochemical reaction models. We apply this idea to a non-linear model of enzyme inhibitor reactions ^[4].

The fractional-order, which involves integration and transect differentiation using fractional calculus is helping to better understand the explanation of real-world problems than ordinary integer order, as well as in the modeling of real phenomena due to a characterization of memory and hereditary properties in ^[5,6]. Riemann Liouville developed the concept of fractional derivative, which is based on power law, ^[7,8] offers a novel fractional derivative that makes use of the exponential kernel. Several issues include the non-singular kernel fractional derivative, which covers the trigonometric and exponential functions, and ^[9,10,11,12] illustrates some relevant techniques for epidemic models. This virus's suggested outbreak efficiently catches the timeline for the COVID-19 disease conceptual model ^[13]. In the literature, many fractional operators are employed to solve real-world issues ^[14,15].

In this paper, section 1 is introduction and section 2 consists of some basic fractional order derivative which are helpful to solve the epidemiological model. Section 3 and 4 consists of generalized from of the model, uniqueness and stability of the model. Fractal Fractional techniques with exponential decay kernel and Mittag-Leffler kernel respectively in section 5. Results and conclusion are discussed in section 6, and 7 respectively.

2. Basic definitions

Following are the basic definitions ^[7,8,14,15] used for analysis and solution of the problem.

Definition 1: Sumudu transform for any function $\phi \left(t\right)$ over a set is given as,

$A = \left\{\phi \left(t\right):there\ exist\ \mathrm{\Lambda }, {\tau }_{1}, {\tau }_{2} > 0, \left|\phi \left(t\right)\right| < \mathrm{\Lambda }exp\left(\frac{\left|t\right|}{{\tau }_{i}}\right), if\ t\in {\left(-1\right)}^{j}\times \left[0, \left.\infty \right)\right.\right\}$

is defined by

$F\left(u\right) = ST\left[\phi \left(t\right)\right] = \underset{0}{\overset{\infty }{\int }}exp\left(-t\right)\phi \left(ut\right)dt, \ \ \ \ \ \ u\in \left({-\tau }_{1}, {\tau }_{2}\right).$

Definition 2: For a function $g\left(t\right)\in {W}_{2}^{1}\left(\mathrm{0, 1}\right), b > a\ and\ \sigma \in \left(\mathrm{0, 1}\right],$ the definition of Atangana–Baleanu derivative in the Caputo sense is given by

${_0^{ABC}}{D}_{t}^{\sigma }g\left(t\right) = \frac{AB\left(\sigma \right)}{1-\sigma }\underset{0}{\overset{t}{\int }}\frac{d}{d\tau }g\left(\tau \right){M}_{\sigma }\left[-\frac{\sigma }{1-\sigma }{\left(t-\tau \right)}^{\sigma }\right]d\tau , \ \ \ \ \ \ n-1 < \sigma < n$

where

$AB\left(\sigma \right) = 1-\sigma +\frac{\sigma }{\mathrm{\Gamma }\left(\sigma \right)}.$

By using Sumudu transform (ST) for (1), we obtain

$ST\left[{_0^{ABC}}{D}_{t}^{\sigma }g\left(t\right)\right]\left(s\right) = \frac{q\left(\sigma \right)}{1-\sigma }\left\{\sigma \mathrm{\Gamma }\left(\sigma +1\right){M}_{\sigma }\left(-\frac{1}{1-\sigma }{V}^{\sigma }\right)\right\}\times \left[ST\left(g\left(t\right)\right)-g\left(0\right)\right].$

Definition 3: For a function $g\left(t\right)\in {W}_{2}^{1}\left(\mathrm{0, 1}\right), b > a\ and\ {\alpha }_{1}\in \left(\mathrm{0, 1}\right],$ the definition of Atangana–Baleanu derivative in the Caputo sense is given by

${_0^{ABC}}{D}_{t}^{{\alpha }_{1}}g\left(t\right) = \frac{AB\left({\alpha }_{1}\right)}{1-{\alpha }_{1}}\underset{0}{\overset{t}{\int }}\frac{d}{d\tau }g\left(\tau \right){E}_{{\alpha }_{1}}[-\frac{{\alpha }_{1}}{1-{\alpha }_{1}}{(t-\tau )}^{{\alpha }_{1}}]d\tau ,$

where

$AB\left({\alpha }_{1}\right) = 1-{\alpha }_{1}+\frac{{\alpha }_{1}}{\mathrm{\Gamma }\left({\alpha }_{1}\right)}.$

Definition 4: Suppose that $g\left(t\right)$ is continuous on an open interval $(a, b),$ then the fractal-fractional integral of $g\left(t\right)$ of order ${\alpha }_{1}$ having Mittag-Leffler type kernel and given by

${{}^{FFM}J}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}\left(g\left(t\right)\right) = \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{s}^{{\alpha }_{2}-1}g\left(s\right){\left(t-s\right)}^{{\alpha }_{1}}ds+\frac{{\alpha }_{2}\left(1-{\alpha }_{1}\right){t}^{{\alpha }_{2}-1}g\left(t\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}$

3. Hires problem with fractional operator

Robertson introduces this chemical process in ^[19,20]. Schafer pioneered the following chemical reactions method in 1975 ^[19,20]. It represents the high irradiance response (HIRES) of photomorphogenesis based on phytochrome. A stiff system of eight non-linear ordinary differential equations is used to create the following mathematical model.

$\begin{array}{l} {{y}_{1}}^{'} = -{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}, \\{{y}_{2}}^{'} = {M}_{1}{y}_{1}-{M}_{5}{y}_{2}, \\ {{y}_{3}}^{'} = -{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}, \\{{y}_{4}}^{'} = {M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}, \\{{y}_{5}}^{'} = -{M}_{10}{y}_{5}+{M}_{2}{y}_{6}+{M}_{2}{y}_{7}, \\{{y}_{6}}^{'} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}, \\{{y}_{7}}^{'} = {M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}, \\{{y}_{8}}^{'} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}. \end{array}$

(1)

Here ${M}_{1} = 1.7, {M}_{2} = 0.43, {M}_{3} = 8.32,$ ${M}_{4} = 0.0007, {M}_{5} = 8.75,$ ${M}_{6} = 10.03, {M}_{7} = 0.035,$ ${M}_{8} = 1.71, {M}_{9} = 1.12, {M}_{10} = 1.745,$ ${M}_{11} = 280, {M}_{12} = 0.69, {M}_{13} = 1.81.$ The initial values can be represented by $y = {\left(\mathrm{1, 0}, \mathrm{0, 0}, \mathrm{0, 0}, \mathrm{0, 0.0057}\right)}^{T}$ . By using Atangana-Baleanu in Caputo sense for system (1), we get

$\begin{array}{l} {_0^{ABC}}{D}_{t}^{\alpha }{y}_{1} = -{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{2} = {M}_{1}{y}_{1}-{M}_{5}{y}_{2}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{3} = -{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{4} = {M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{5} = -{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{6} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{7} = {M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}, \\ {_0^{ABC}}{D}_{t}^{\alpha }{y}_{8} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}. \end{array}$

(2)

Here ${}_{O}{}^{ABC}{D}_{t}^{\alpha }$ is the Atanagana-Baleanue Caputo sense fractional derivative with $0 < \alpha \le 1$ .

With given initial conditions

${y}_{i}\left(0\right)\ge 0, i = \mathrm{1, 2}, 3, \dots , 8$

(3)

Theorem 3.1: The solution of the proposed fractional-order model (1) along initial conditions is unique and bounded in R₊⁸.

Proof: In (1), we can get its existence and uniqueness on the time interval (0, ∞). Afterwards, we need to show that the non-negative region R₊⁸ is a positively invariant region. For this

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{1 = 0} = {M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{2 = 0} = {M}_{1}{y}_{1}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{3 = 0} = {M}_{2}{y}_{4}+{M}_{7}{y}_{5}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{4 = 0} = {M}_{3}{y}_{2}+{M}_{8}{y}_{3}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{5 = 0} = {M}_{2}{y}_{6}-{M}_{2}{y}_{7}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{6 = 0} = {M}_{12}{y}_{4}+{M}_{8}{y}_{5}+{M}_{12}{y}_{7}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{7 = 0} = {M}_{11}{y}_{6}{y}_{8}\ge 0,$

${_0^{ABC}}{D}_{t}^{\alpha }{y}_{1}|{y}_{8 = 0} = {M}_{13}{y}_{7}\ge 0$

If ( ${y}_{1}$ (0)), ( ${y}_{2}$ (0)), ( ${y}_{3}$ (0)), ( ${y}_{4}$ (0)), ( ${y}_{5}$ (0)), ( ${y}_{6}$ (0)), ( ${y}_{7}$ (0)), ( ${y}_{8}$ (0)) $\epsilon$ ${R}_{+}^{8}$ , then from above expression, the solution cannot escape from the hyperplane. Also on each hyperplane bounding the non-negative orthant, the vector field points into ${R}_{+}^{8}$ , i.e., the domain ${R}_{+}^{8}$ is a positively invariant set.

Now, with the help of Sumudu transform definition, we get

$\begin{array}{l} Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{1}\left(t\right)-{y}_{1}\left(0\right)\right\} = ST\left[-{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{2}\left(t\right)-{y}_{2}\left(0\right)\right\} = ST\left[{M}_{1}{y}_{1}-{M}_{5}{y}_{2}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{3}\left(t\right)-{y}_{3}\left(0\right)\right\} = ST\left[-{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{4}\left(t\right)-{y}_{4}\left(0\right)\right\} = ST\left[{M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{5}\left(t\right)-{y}_{5}\left(0\right)\right\} = ST\left[-{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{6}\left(t\right)-{y}_{6}\left(0\right)\right\} = ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{7}\left(t\right)-{y}_{7}\left(0\right)\right\} = ST\left[{M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}\right], \\ Q{E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)ST\left\{{y}_{8}\left(t\right)-{y}_{8}\left(0\right)\right\} = ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}\right]. \end{array}$

(4)

Where $Q = \frac{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right)}{1-\alpha }$

Rearranging, we get

$\begin{array}{l} ST\left({y}_{1}\left(t\right)\right) = {y}_{1}\left(0\right)+H\times ST\left[-{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}\right], \\ ST\left({y}_{2}\left(t\right)\right) = {y}_{2}\left(0\right)+H\times ST\left[{M}_{1}{y}_{1}-{M}_{5}{y}_{2}\right], \\ ST\left({y}_{3}\left(t\right)\right) = {y}_{3}\left(0\right)+H\times ST\left[-{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}\right], \\ ST\left({y}_{4}\left(t\right)\right) = {y}_{4}\left(0\right)+H\times ST\left[{M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}\right], \\ ST\left({y}_{5}\left(t\right)\right) = {y}_{5}\left(0\right)+H\times ST\left[-{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}\right], \\ ST\left({y}_{6}\left(t\right)\right) = {y}_{6}\left(0\right)+H\times ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}\right], \\ ST\left({y}_{7}\left(t\right)\right) = {y}_{7}\left(0\right)+H\times ST\left[{M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}\right], \\ ST\left({y}_{8}\left(t\right)\right) = {y}_{8}\left(0\right)+H\times ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}\right]. \end{array}$

(5)

Using the inverse Sumudu transform on both sides of the system (5), we obtain

$\begin{array}{l} {y}_{1}\left(t\right) = {y}_{1}\left(0\right)+{ST}^{-1}\left[H\times ST\left[-{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}\right]\right], \\ {y}_{2}\left(t\right) = {y}_{2}\left(0\right)+{ST}^{-1}\left[H\times ST\left[{M}_{1}{y}_{1}-{M}_{5}{y}_{2}\right]\right], \\ {y}_{3}\left(t\right) = {y}_{3}\left(0\right)+{ST}^{-1}\left[H\times ST\left[-{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}\right]\right], \\ {y}_{4}\left(t\right) = {y}_{4}\left(0\right)+{ST}^{-1}\left[H\times ST\left[{M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}\right]\right], \\ {y}_{5}\left(t\right) = {y}_{5}\left(0\right)+{ST}^{-1}\left[H\times ST\left[-{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}\right]\right], \\ {y}_{6}\left(t\right) = {y}_{6}\left(0\right)+{ST}^{-1}\left[H\times ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}\right]\right], \\ {y}_{7}\left(t\right) = {y}_{7}\left(0\right)+{ST}^{-1}\left[H\times ST\left[{M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}\right]\right], \\ {y}_{8}\left(t\right) = {y}_{8}\left(0\right)+{ST}^{-1}\left[H\times ST\left[-{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}\right]\right]. \end{array}$

(6)

We next obtain the following recursive formula.

$\begin{array}{l} {{y}_{1}}_{\left(n+1\right)}\left(t\right) = {{y}_{1}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{-{M}_{1}{{y}_{1}}_{\left(n\right)}+{M}_{2}{{y}_{2}}_{\left(n\right)}+{M}_{3}{{y}_{3}}_{\left(n\right)}+{M}_{4}\right\}\right], \\ {{y}_{2}}_{\left(n+1\right)}\left(t\right) = {{y}_{2}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{{M}_{1}{{y}_{1}}_{\left(n\right)}-{M}_{5}{{y}_{2}}_{\left(n\right)}\right\}\right], \\ {{y}_{3}}_{\left(n+1\right)}\left(t\right) = {{y}_{3}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{-{M}_{6}{{y}_{3}}_{\left(n\right)}+{M}_{2}{{y}_{4}}_{\left(n\right)}+{M}_{7}{{y}_{5}}_{\left(n\right)}\right\}\right], \\ {{y}_{4}}_{\left(n+1\right)}\left(t\right) = {{y}_{4}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{{M}_{3}{{y}_{2}}_{\left(n\right)}+{M}_{8}{{y}_{3}}_{\left(n\right)}-{M}_{9}{{y}_{4}}_{\left(n\right)}\right\}\right], \\ {{y}_{5}}_{\left(n+1\right)}\left(t\right) = {{y}_{5}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{-{M}_{10}{{y}_{5}}_{\left(n\right)}+{M}_{2}{{y}_{6}}_{\left(n\right)}-{M}_{2}{{y}_{7}}_{\left(n\right)}\right\}\right], \\ {{y}_{6}}_{\left(n+1\right)}\left(t\right) = {{y}_{6}}_{\left(n\right)}\left(0\right)+{ST}^{-1}[H\times ST\{-{M}_{11}{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}+{M}_{12}{{y}_{4}}_{\left(n\right)}+{M}_{8}{{y}_{5}}_{\left(n\right)}-{M}_{2}{{y}_{6}}_{\left(n\right)}+\\ {M}_{12}{{y}_{7}}_{\left(n\right)}\}], \\ {{y}_{7}}_{\left(n+1\right)}\left(t\right) = {{y}_{7}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{{M}_{11}{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{M}_{13}{{y}_{7}}_{\left(n\right)}\right\}\right], \\ {{y}_{8}}_{\left(n+1\right)}\left(t\right) = {{y}_{8}}_{\left(n\right)}\left(0\right)+{ST}^{-1}\left[H\times ST\left\{-{M}_{11}{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}+{M}_{13}{{y}_{7}}_{\left(n\right)}\right\}\right]. \end{array}$

(7)

Where $H = \frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}$

And the solution of system is provided by

${y}_{1}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{1}}_{\left(n\right)}\left(t\right), \ \ \ \ {y}_{2}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{2}}_{\left(n\right)}\left(t\right), \ \ \ \ {y}_{3}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{3}}_{\left(n\right)}\left(t\right),$

${y}_{4}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{4}}_{\left(n\right)}\left(t\right), \ \ \ \ {y}_{5}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{5}}_{\left(n\right)}\left(t\right), \ \ \ \ {y}_{6}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{6}}_{\left(n\right)}\left(t\right),$

${y}_{7}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{7}}_{\left(n\right)}\left(t\right), \ \ \ \ {y}_{8}\left(t\right) = \underset{n\to \infty }{\mathrm{lim}}{{y}_{8}}_{\left(n\right)}\left(t\right).$

4. Fixed-point theorem for stability analysis of iteration method

Theorem 4.1: Define K be a self-map is given by

$\begin{array}{l} K\left[ {{y_1}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_1}_{\left( {n + 1} \right)}\left( t \right) = {y_1}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ - {M_1}{y_1}_{\left( n \right)} + \\ {M_2}{y_2}_{\left( n \right)} + {M_3}{y_3}_{\left( n \right)} + {M_4}\} ], \\ K\left[ {{y_2}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_2}_{\left( {n + 1} \right)}\left( t \right) = {y_2}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ {M_1}{y_1}_{\left( n \right)}-\\ {M_5}{y_2}_{\left( n \right)}\} ], \\ K\left[ {{y_3}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_3}_{\left( {n + 1} \right)}\left( t \right) = {y_3}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ - {M_6}{y_3}_{\left( n \right)} + \\ {M_2}{y_4}_{\left( n \right)} + {M_7}{y_5}_{\left( n \right)}\} ], \\ K\left[ {{y_4}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_4}_{\left( {n + 1} \right)}\left( t \right) = {y_4}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ {M_3}{y_2}_{\left( n \right)} + {M_8}{y_3}_{\left( n \right)} - \\ {M_9}{y_4}_{\left( n \right)}\} ], \\ K\left[ {{y_5}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_5}_{\left( {n + 1} \right)}\left( t \right) = {y_5}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ - {M_{10}}{y_5}_{\left( n \right)} + \\ {M_2}{y_{6(n)}} - {M_2}{y_7}_{\left( n \right)}\} ], \\ K\left[ {{y_6}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_6}_{\left( {n + 1} \right)}\left( t \right) = {y_6}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ - {M_{11}}{y_6}_{\left( n \right)}{y_8}_{\left( n \right)} + \\ {M_{12}}{y_4}_{\left( n \right)} + {M_8}{y_5}_{\left( n \right)} - {M_2}{y_6}_{\left( n \right)} + {M_{12}}{y_7}_{\left( n \right)}\} ], \\ K\left[ {{y_7}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_7}_{\left( {n + 1} \right)}\left( t \right) = {y_7}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ {M_{11}}{y_6}_{\left( n \right)}{y_8}_{\left( n \right)} - \\ {M_{13}}{y_7}_{\left( n \right)}\} ], \\ K\left[ {{y_8}_{\left( {n + 1} \right)}\left( t \right)} \right] = {y_8}_{\left( {n + 1} \right)}\left( t \right) = {y_8}_{\left( n \right)}\left( 0 \right) + S{T^{ - 1}}[\frac{{1 - \alpha }}{{M\left( \alpha \right)\alpha \Gamma \left( {\alpha + 1} \right){E_\alpha }\left( { - \frac{1}{{1 - \alpha }}{P^\alpha }} \right)}} \times ST\{ - {M_{11}}{y_6}_{\left( n \right)}{y_8}_{\left( n \right)} + \\ {M_{13}}{y_7}_{\left( n \right)}\} ]. \end{array}$

(8)

Proof: By using triangular inequality with the definition of norms, we get

$\begin{array}{l} ||K\left[{{y}_{1}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{1}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{1}}_{\left(n\right)}\left(t\right)-{{y}_{1}}_{\left(m\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times\\ ST\left\{-{M}_{1}||{{y}_{1}}_{\left(n\right)}-{{y}_{1}}_{\left(m\right)}||+{M}_{2}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||+{M}_{3}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||+{M}_{4}\right\}], \\ ||K\left[{{y}_{2}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{2}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{2}}_{\left(n\right)}\left(t\right)-{{y}_{2}}_{\left(m\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times \\ ST\left\{{M}_{1}||{{y}_{1}}_{\left(n\right)}-{{y}_{1}}_{\left(m\right)}||-{M}_{5}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||\right\}], \\ ||K\left[{{y}_{3}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{3}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{3}}_{\left(n\right)}\left(t\right)-{{y}_{3}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times\\ ST\left\{-{M}_{6}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||+{M}_{2}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||+{M}_{7}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||\right\}], \\ ||K\left[{{y}_{4}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{4}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{4}}_{\left(n\right)}\left(t\right)-{{y}_{4}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times \\ ST\left\{{M}_{3}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||+{M}_{8}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||-{M}_{9}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||\right\}], \\ ||K\left[{{y}_{5}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{5}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{5}}_{\left(n\right)}\left(t\right)-{{y}_{5}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times \\ ST\left\{-{M}_{10}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||+{M}_{2}||{{y}_{6}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}||-{M}_{2}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\right\}], \\ ||K\left[{{y}_{6}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{6}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{6}}_{\left(n\right)}\left(t\right)-{{y}_{6}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times\\ ST\{-{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||+{M}_{12}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||+{M}_{8}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||- \\ {M}_{2}||{{y}_{6}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}||+{M}_{12}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\}], \\ ||K\left[{{y}_{7}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{7}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{7}}_{\left(n\right)}\left(t\right)-{{y}_{7}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times\\ ST\left\{{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||-{M}_{13}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\right\}], \\ ||K\left[{{y}_{8}}_{\left(n\right)}\left(t\right)\right]-K\left[{{y}_{8}}_{\left(m\right)}\left(t\right)\right]||\le ||{{y}_{8}}_{\left(n\right)}\left(t\right)-{{y}_{8}}_{\left(n\right)}\left(t\right)||+{ST}^{-1}[\frac{1-\alpha }{M\left(\alpha \right)\alpha \mathrm{\Gamma }\left(\alpha +1\right){E}_{\alpha }\left(-\frac{1}{1-\alpha }{P}^{\alpha }\right)}\times\\ ST\left\{-{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||+{M}_{13}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\right\}]. \end{array}$

Hence satisfied given conditions.

$\begin{array}{l} \theta = \left(\mathrm{0, 0}, \mathrm{0, 0}, \mathrm{0, 0}, \mathrm{0, 0}\right), \theta = \\\left\{\begin{array}{c}||{{y}_{1}}_{\left(n\right)}\left(t\right)-{{y}_{1}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{1}}_{\left(n\right)}\left(t\right)+{{y}_{1}}_{\left(m\right)}\left(t\right)\right)||-{M}_{1}||{{y}_{1}}_{\left(n\right)}-{{y}_{1}}_{\left(m\right)}||+{M}_{2}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||\\ +{M}_{3}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||+{M}_{4}\\ ||{{y}_{2}}_{\left(n\right)}\left(t\right)-{{y}_{2}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{2}}_{\left(n\right)}\left(t\right)+{{y}_{2}}_{\left(m\right)}\left(t\right)\right)||+{M}_{1}||{{y}_{1}}_{\left(n\right)}-{{y}_{1}}_{\left(m\right)}||-{M}_{5}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||\\ ||{{y}_{3}}_{\left(n\right)}\left(t\right)-{{y}_{3}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{3}}_{\left(n\right)}\left(t\right)+{{y}_{3}}_{\left(m\right)}\left(t\right)\right)||-{M}_{6}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||+{M}_{2}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||\\ +{M}_{7}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||\\ ||{{y}_{4}}_{\left(n\right)}\left(t\right)-{{y}_{4}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{4}}_{\left(n\right)}\left(t\right)+{{y}_{4}}_{\left(m\right)}\left(t\right)\right)||+{M}_{3}||{{y}_{2}}_{\left(n\right)}-{{y}_{2}}_{\left(m\right)}||+{M}_{8}||{{y}_{3}}_{\left(n\right)}-{{y}_{3}}_{\left(m\right)}||\\ -{M}_{9}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||\\ ||{{y}_{5}}_{\left(n\right)}\left(t\right)-{{y}_{5}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{5}}_{\left(n\right)}\left(t\right)+{{y}_{5}}_{\left(m\right)}\left(t\right)\right)||-{M}_{10}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||+{M}_{2}||{{y}_{6}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}||\\ -{M}_{2}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\\ ||{{y}_{6}}_{\left(n\right)}\left(t\right)-{{y}_{6}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{6}}_{\left(n\right)}\left(t\right)+{{y}_{6}}_{\left(m\right)}\left(t\right)\right)||-{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||\\ +{M}_{12}||{{y}_{4}}_{\left(n\right)}-{{y}_{4}}_{\left(m\right)}||+{M}_{8}||{{y}_{5}}_{\left(n\right)}-{{y}_{5}}_{\left(m\right)}||-{M}_{2}||{{y}_{6}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}||+{M}_{12}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\\ ||{{y}_{7}}_{\left(n\right)}\left(t\right)-{{y}_{7}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{7}}_{\left(n\right)}\left(t\right)+{{y}_{7}}_{\left(m\right)}\left(t\right)\right)||+{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||\\ -{M}_{13}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\\ ||{{y}_{8}}_{\left(n\right)}\left(t\right)-{{y}_{8}}_{\left(m\right)}\left(t\right)||\times ||-\left({{y}_{8}}_{\left(n\right)}\left(t\right)+{{y}_{8}}_{\left(m\right)}\left(t\right)\right)||-{M}_{11}||{{y}_{6}}_{\left(n\right)}{{y}_{8}}_{\left(n\right)}-{{y}_{6}}_{\left(m\right)}{{y}_{8}}_{\left(m\right)}||\\ +{M}_{13}||{{y}_{7}}_{\left(n\right)}-{{y}_{7}}_{\left(m\right)}||\end{array}\right. \end{array}$

Hence the system is stable.

Theorem 4.2: Unique singular solution with the iterative method for the special solution of system $\left(2\right).$

Proof: Considering the Hilbert space $H = {L}^{2}\left(\left(p, q\right)\times \left(0, T\right)\right)$ which can be defined as

$h:\left(p, q\right)\times \left(0, T\right)\to \mathbb{R}, \iint ghdgdh < \infty .$

For this purpose, we consider the following operator

$\theta \left(\mathrm{0, 0}, \mathrm{0, 0}, \mathrm{0, 0}, \mathrm{0, 0}\right), \theta = \left\{\begin{array}{c}-{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}, \\ {M}_{1}{y}_{1}-{M}_{5}{y}_{2}, \\ -{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}, \\ {M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}, \\ -{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}, \\ -{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}, \\ {M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}, \\ -{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}.\end{array}\right.$

By using inner product, we get

$\begin{array}{l} T(({{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}, {{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}, {{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}, {{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}, {{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}, {{y}_{6}}_{\left(61\right)}-\\ \ \ \ \ \ \ \ \ {{y}_{6}}_{\left(62\right)}, {{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}, {{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}), \left({V}_{1}, {V}_{2}, {V}_{3}, {V}_{4}, {V}_{5}, {V}_{6}, {V}_{7}, {V}_{8}\right)). \end{array}$

Where

$({{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}, {{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}, {{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)},$ ${{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}, {{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)},$ ${{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)},$ ${{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)},$ ${{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}),$ are the special solutions of the system. Taking into account the inner function and the norm, we have

$\begin{array}{l} \left\{ { - {M_1}\left( {{y_1}_{\left( {11} \right)} - {y_1}_{\left( {12} \right)}} \right) + {M_2}\left( {{y_2}_{\left( {21} \right)} - {y_2}_{\left( {22} \right)}} \right) + {M_3}\left( {{y_3}_{\left( {31} \right)} - {y_3}_{\left( {32} \right)}} \right) + {M_4}, {V_1}} \right\} \le {M_1}||{y_1}_{\left( {11} \right)} - \\ {y_1}_{\left( {12} \right)}||||{V_1}|| + {M_2}||{y_2}_{\left( {21} \right)} - {y_2}_{\left( {22} \right)}||||{V_1}|| + {M_3}||{y_3}_{\left( {31} \right)} - {y_3}_{\left( {32} \right)}||||{V_1}|| + {M_4}||{V_1}||, \\ \left\{ {{M_1}\left( {{y_1}_{\left( {11} \right)} - {y_1}_{\left( {12} \right)}} \right) - {M_5}\left( {{y_2}_{\left( {21} \right)} - {y_2}_{\left( {22} \right)}} \right), {V_2}} \right\} \le {M_1}||{y_1}_{\left( {11} \right)} - {y_1}_{\left( {12} \right)}||||{V_2}|| + {M_5}||{y_2}_{\left( {21} \right)} - \\ {y_2}_{\left( {22} \right)}||||{V_2}||, \\ \left\{ { - {M_6}\left( {{y_3}_{\left( {31} \right)} - {y_3}_{\left( {32} \right)}} \right) + {M_2}\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right) + {M_7}\left( {{y_5}_{\left( {51} \right)} - {y_5}_{\left( {52} \right)}} \right), {V_3}} \right\} \le {M_6}||({y_3}_{\left( {31} \right)} - \\ {y_3}_{\left( {32} \right)})||||{V_3}|| + {M_2}||\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right)||||{V_3}|| + {M_7}||\left( {{y_5}_{\left( {51} \right)} - {y_5}_{\left( {52} \right)}} \right)||||{V_3}||, \\ \left\{ {{M_3}\left( {{y_2}_{\left( {21} \right)} - {y_2}_{\left( {22} \right)}} \right) + {M_8}\left( {{y_3}_{\left( {31} \right)} - {y_3}_{\left( {32} \right)}} \right) - {M_9}\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right), {V_4}} \right\} \le {M_3}||({y_2}_{\left( {21} \right)} - \\ {y_2}_{\left( {22} \right)})||||{V_4}|| + {M_8}||\left( {{y_3}_{\left( {31} \right)} - {y_3}_{\left( {32} \right)}} \right)||||{V_4}|| + {M_9}||\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right)||||{V_4}||, \\ \left\{ { - {M_{10}}\left( {{y_5}_{\left( {51} \right)} - {y_5}_{\left( {52} \right)}} \right) + {M_2}\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right) - {M_2}\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right), {V_5}} \right\} \le {M_{10}}||({y_5}_{\left( {51} \right)} - \\ {y_5}_{\left( {52} \right)})||||{V_5}|| + {M_2}||\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right)||||{V_5}|| + {M_2}||\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right)||||{V_5}||, \\ \{ - {M_{11}}\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right)\left( {{y_8}_{\left( {81} \right)} - {y_8}_{\left( {82} \right)}} \right) + {M_{12}}\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right) + {M_8}\left( {{y_5}_{\left( {51} \right)} - {y_5}_{\left( {52} \right)}} \right) - \\ {M_2}\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right) + {M_{12}}\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right), {V_6}\} \le {M_{11}}||\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right)||||({y_8}_{\left( {81} \right)} - \\ {y_8}_{\left( {82} \right)})||||{V_6}|| + {M_{12}}||\left( {{y_4}_{\left( {41} \right)} - {y_4}_{\left( {42} \right)}} \right)||||{V_6}|| + {M_8}||\left( {{y_5}_{\left( {51} \right)} - {y_5}_{\left( {52} \right)}} \right)||||{V_6}|| + {M_2}||({y_6}_{\left( {61} \right)} - \\ {y_6}_{\left( {62} \right)})||||{V_6}|| + {M_{12}}||\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right)||||{V_6}||, \\ \left\{ {{M_{11}}\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right)\left( {{y_8}_{\left( {81} \right)} - {y_8}_{\left( {82} \right)}} \right) - {M_{13}}\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right), {V_7}} \right\} \le {M_{11}}||({y_6}_{(61)} - \\ {y_6}_{\left( {62} \right)})||||\left( {{y_8}_{\left( {81} \right)} - {y_8}_{\left( {82} \right)}} \right)||||{V_7}|| + {M_{13}}||\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right)||||{V_7}||, \\ \left\{ { - {M_{11}}\left( {{y_6}_{\left( {61} \right)} - {y_6}_{\left( {62} \right)}} \right)\left( {{y_8}_{\left( {81} \right)} - {y_8}_{\left( {82} \right)}} \right) + {M_{13}}\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right), {V_8}} \right\} \le {M_{11}}||({y_6}_{(61)} - \\ {y_6}_{\left( {62} \right)})||||\left( {{y_8}_{\left( {81} \right)} - {y_8}_{\left( {82} \right)}} \right)||||{V_8}|| + {M_{13}}||\left( {{y_7}_{\left( {71} \right)} - {y_7}_{\left( {72} \right)}} \right)||||{V_8}||. \end{array}$

In the case for large number ${e}_{1}, {e}_{2}, {e}_{3}, {e}_{4}, {e}_{5}, {e}_{6}, {e}_{7}and{e}_{8},$ both solutions happen to be converged to the exact solution. Employing the topology concept, we can obtain eight positive very small parameters $\left({\chi }_{{e}_{1}}, {\chi }_{{e}_{2}}, {\chi }_{{e}_{3}}, {\chi }_{{e}_{4}}, {\chi }_{{e}_{5}}, {\chi }_{{e}_{6}}, {\chi }_{{e}_{7}}and{\chi }_{{e}_{8}}\right).$

$||{y}_{1}-{{y}_{1}}_{\left(11\right)}||, ||{y}_{1}-{{y}_{1}}_{\left(12\right)}||\le \frac{{\chi }_{{e}_{1}}}{\varpi }, ||{y}_{2}-{{y}_{2}}_{\left(21\right)}||, ||{y}_{2}-{{y}_{2}}_{\left(22\right)}||\le \frac{{\chi }_{{e}_{2}}}{\varsigma },$

$||{y}_{3}-{{y}_{3}}_{\left(31\right)}||, ||{y}_{3}-{{y}_{3}}_{\left(32\right)}||\le \frac{{\chi }_{{e}_{3}}}{\upsilon }, ||{y}_{4}-{{y}_{4}}_{\left(41\right)}||, ||{y}_{4}-{{y}_{4}}_{\left(42\right)}||\le \frac{{\chi }_{{e}_{4}}}{\kappa },$

$||{y}_{5}-{{y}_{5}}_{\left(51\right)}||, ||{y}_{5}-{{y}_{5}}_{\left(52\right)}||\le \frac{{\chi }_{{e}_{5}}}{\varrho }, ||{y}_{6}-{{y}_{6}}_{\left(61\right)}||, ||{y}_{6}-{{y}_{6}}_{\left(62\right)}||\le \frac{{\chi }_{{e}_{6}}}{\zeta },$

$||{y}_{7}-{{y}_{7}}_{\left(71\right)}||, ||{y}_{7}-{{y}_{7}}_{\left(72\right)}||\le \frac{{\chi }_{{e}_{7}}}{\nu }, ||{y}_{8}-{{y}_{8}}_{\left(81\right)}||, ||{y}_{8}-{{y}_{8}}_{\left(82\right)}||\le \frac{{\chi }_{{e}_{8}}}{\varepsilon }.$

Where

$\varpi = 8\left({M}_{1}||{{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}||+{M}_{2}||{{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}||+{M}_{3}||{{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}||+{M}_{4}\right)||{V}_{1}||$

$\varsigma = 8\left({M}_{1}||{{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}||+{M}_{5}||{{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}||\right)||{V}_{2}||$

$\upsilon = 8\left({M}_{6}||\left({{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}\right)||+{M}_{2}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||+{M}_{7}||\left({{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}\right)||\right)||{V}_{3}||$

$\kappa = 8\left({M}_{3}||\left({{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}\right)||+{M}_{8}||\left({{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}\right)||+{M}_{9}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||\right)||{V}_{4}||$

$\varrho = 8\left({M}_{10}||\left({{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}\right)||+{M}_{2}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||+{M}_{2}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||\right)||{V}_{5}||$

$\begin{array}{l} \zeta = 8({M}_{11}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)||+{M}_{12}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||+{M}_{8}||({{y}_{5}}_{\left(51\right)}-\\ {{y}_{5}}_{\left(52\right)})||+{M}_{2}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||+{M}_{12}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||)||{V}_{6}|| \end{array}$

$\nu = 8\left({M}_{11}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)||+{M}_{13}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||\right)||{V}_{7}||$

$\varepsilon = 8\left({M}_{11}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)||+{M}_{13}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||\right)||{V}_{8}||$

But, it is obvious that

$\left({M}_{1}||{{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}||+{M}_{2}||{{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}||+{M}_{3}||{{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}||+{M}_{4}\right)\ne 0$

$\left({M}_{1}||{{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}||+{M}_{5}||{{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}||\right)\ne 0$

$\left({M}_{6}||\left({{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}\right)||+{M}_{2}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||+{M}_{7}||\left({{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}\right)||\right)\ne 0$

$\left({M}_{3}||\left({{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}\right)||+{M}_{8}||\left({{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}\right)||+{M}_{9}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||\right)\ne 0$

$\left({M}_{10}||\left({{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}\right)||+{M}_{2}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||+{M}_{2}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||\right)\ne 0$

$\begin{array}{l} ({M}_{11}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)||+{M}_{12}||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)||+{M}_{8}||({{y}_{5}}_{\left(51\right)}-\\ {{y}_{5}}_{\left(52\right)})||+{M}_{2}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||+{M}_{12}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||)\ne 0 \end{array}$

$\left({M}_{11}||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)||||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)||+{M}_{13}||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)||\right)\ne 0$

where $||{V}_{1}||, ||{V}_{2}||, ||{V}_{3}||, ||{V}_{4}||, ||{V}_{5}||, ||{V}_{6}||, ||{V}_{7}||, ||{V}_{8}||\ne 0$ .

Therefore, we have

$||{{y}_{1}}_{\left(11\right)}-{{y}_{1}}_{\left(12\right)}|| = 0, ||{{y}_{2}}_{\left(21\right)}-{{y}_{2}}_{\left(22\right)}|| = 0, ||{{y}_{3}}_{\left(31\right)}-{{y}_{3}}_{\left(32\right)}|| = 0,$

$||\left({{y}_{4}}_{\left(41\right)}-{{y}_{4}}_{\left(42\right)}\right)|| = 0, ||\left({{y}_{5}}_{\left(51\right)}-{{y}_{5}}_{\left(52\right)}\right)|| = 0, ||\left({{y}_{6}}_{\left(61\right)}-{{y}_{6}}_{\left(62\right)}\right)|| = 0,$

$||\left({{y}_{7}}_{\left(71\right)}-{{y}_{7}}_{\left(72\right)}\right)|| = 0, ||\left({{y}_{8}}_{\left(81\right)}-{{y}_{8}}_{\left(82\right)}\right)|| = 0.$

Which yields that

$\begin{array}{*{20}{c}} { {{y}_{1}}_{\left(11\right)} = {{y}_{1}}_{\left(12\right)}, {{y}_{2}}_{\left(21\right)} = {{y}_{2}}_{\left(22\right)}, {{y}_{3}}_{\left(31\right)} = {{y}_{3}}_{\left(32\right)}, {{y}_{4}}_{\left(41\right)} = {{y}_{4}}_{\left(42\right)}, {{y}_{5}}_{\left(51\right)} = {{y}_{5}}_{\left(52\right)},{{y}_{6}}_{\left(61\right)} = }\\ {{{y}_{6}}_{\left(62\right)}, {{y}_{7}}_{\left(71\right)} = {{y}_{7}}_{\left(72\right)}, {{y}_{8}}_{\left(81\right)} = {{y}_{8}}_{\left(82\right)} } \end{array}$

This completes the proof of uniqueness.

An operator $B:Z\to Z$ can be defined as:

$B\left(\varphi \right)\left(t\right) = \varphi \left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}(1-{\alpha }_{1})}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, \varphi \left(t\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(t, \varphi \left(t\right)\right)d\lambda$

(10)

If $£ \left(t, \varphi \left(t\right)\right)$ satisfies the Lipschitz condition and the following extension then

● For every $\varphi \in Z$ there exists constants ${L}_{£ } > 0$ and ${M}_{£ }$ such that

$\left|£ \left(t, \varphi \left(t\right)\right)|\le {L}_{£ }\right|\varphi \left(t\right)|+{M}_{£ }$

(11)

● For every $\varphi, \overline{\varphi }\in Z$ , there exists a constant ${M}_{£ } > 0$ such that

$\left|£ \left(t, \varphi \left(t\right)\right)-£ \left(t, \overline{\varphi }\left(t\right)\right)\left|\right|\le {M}_{£ }\right|\varphi \left(t\right)-\overline{\varphi }\left(t\right)|$

(12)

Theorem 4.2: If the condition of (11) holds then for the function $£ :[0, T]\times Z\to R$ there exists at least one solution for the (1).

Proof: Since $£$ in (10) is continuous function, so $B$ is also a continuous. Assume $M = \left\{\varphi \in \left|\left|\varphi \right|\right|\le R, R > 0\right\}$ , then for $\varphi \in Z,$ we have

$\begin{array}{l} B\left(\varphi \right)\left(t\right) = {max}_{t\in \left[0, T\right]}|\varphi \left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, \varphi \left(t\right)\right)\\ \;\;\;\;\;\;\;\;\;\;+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(t, \varphi \left(t\right)\right)d\lambda \end{array}$

$\begin{array}{l} \le |\varphi \left(0\right)+\frac{\mathsf{μ} {T}^{\mathsf{μ} -1}(1-{\alpha }_{1})}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}\left({L}_{£ }\left|\right|\varphi \left(t\right)\left|\right|+{M}_{£ }\right)\\ \;\;\;\;\;\;\;\;\;\;+{max}_{t\in [0, T]}\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(t, \varphi \left(t\right)\right)d\lambda | \end{array}$

$\begin{array}{l} \le \varphi \left(0\right)+\frac{\mathsf{μ} {T}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}\left({L}_{£ }\left|\left|\varphi \left(t\right)\right|\right|+{M}_{£ }\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{\left({L}_{£ }\left|\left|\varphi \left(t\right)\right|\right|+{M}_{£ }\right)T}^{\mathsf{μ} +{\alpha }_{1-1}}M\left(\mathsf{μ} , {\alpha }_{1}\right)\\ \le R. \end{array}$

Hence, $B$ is uniformly bounded, and $M\left(\mathsf{μ}, {\alpha }_{1}\right)$ is a beta function. For equicontinuity of $B$ , we take ${t}_{1} < {t}_{2}\le T,$ then consider

$B\left(\varphi \right)\left({t}_{2}\right)-B\left(\varphi \right)\left({t}_{1}\right) = |\frac{\mathsf{μ} {{t}_{2}}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left({t}_{2}, \varphi \left({t}_{2}\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}$

$\begin{array}{l} \underset{0}{\overset{{t}_{2}}{\int }}{\lambda }^{\mathsf{μ} -1}{\left({t}_{2}-\lambda \right)}^{\mathsf{μ} -1}£ \left(t, \varphi \left(t\right)\right)d\lambda -\frac{\mathsf{μ} {{t}_{1}}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left({t}_{1}, \varphi \left({t}_{1}\right)\right)\\ \ \ \ \ \ \ \ \ \ \ \ \ +\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{{t}_{2}}{\int }}{\lambda }^{\mathsf{μ} -1}{\left({t}_{1}-\lambda \right)}^{\mathsf{μ} -1}£ \left(t, \varphi \left(t\right)\right)d\lambda | \end{array}$

$\le \frac{\mathsf{μ} {{t}_{2}}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}\left({L}_{£ }\left|\varphi \left(t\right)\right|+{M}_{£ }\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{\left({L}_{£ }\left|\varphi \left(t\right)\right|+{M}_{£ }\right){t}_{2}}^{\mathsf{μ} +{\alpha }_{1-1}}M\left(\mathsf{μ} , {\alpha }_{1}\right)$

$\frac{\mathsf{μ} {{t}_{1}}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}\left({L}_{£ }\left|\varphi \left(t\right)\right|+{M}_{£ }\right)-\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{\left({L}_{£ }\left|\varphi \left(t\right)\right|+{M}_{£ }\right){t}_{1}}^{\mathsf{μ} +{\alpha }_{1-1}}M\left(\mathsf{μ} , {\alpha }_{1}\right)$

If ${t}_{1}\to {t}_{2}$ then $\left|\left|B\left(\varphi \right)\left({t}_{2}\right)-B\left(\varphi \right)\left({t}_{1}\right)\to 0\right|\right|$ Consequently $\left|\left|B\left(\varphi \right)\left({t}_{2}\right)-B\left(\varphi \right)\left({t}_{1}\right)\to 0\right|\right|, as{t}_{1}\to {t}_{2}.$ Hence $B$ is equicontinous. Thus, by Arzela-Ascoli theorem $B$ is completely continuous. Consequently, by the result of Schauder's fixed point, it has at least one solution.

Theorem 4.3: If $\eta = \frac{\mathsf{μ} {T}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{T}^{\mathsf{μ} +{\alpha }_{1}-1}M\left(\mathsf{μ}, {\alpha }_{1}\right){M}_{£ } < 1$ and the condition (12) holds, then $\eta$ has a unique solution.

Proof: For $\varphi, \overline{\varphi }\in Z$ , we have

$\left|B\left(\varphi \right)-B\left(\overline{\varphi }\right)\right| = {max}_{t\in \left[0, T\right]}|\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}\left[£ \left(t, \varphi \left(t\right)\right)|-£ \left(t, \overline{\varphi }\left(t\right)\right)\right]$

$+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}\left[£ \left(t, \varphi \left(t\right)\right)|-£ \left(t, \overline{\varphi }\left(t\right)\right)\right]d\lambda |$

$\le [\frac{\mathsf{μ} {T}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{T}^{\mathsf{μ} +{\alpha }_{1}-1}M\left(\mathsf{μ} , {\alpha }_{1}\right)]\left|\left|B\left(\varphi \right)-B\left(\overline{\varphi }\right)\right|\right|$

$\le \eta \left|\left|B\left(\varphi \right)-B\left(\overline{\varphi }\right)\right|\right|.$

Hence, $B$ is a contraction. So, by the principle of Banach contraction, it has a unique solution.

Ulam-Hyres stability

The proposed model is Ulam-Hyres stable if there exists ${B}_{\mathsf{μ}, {\alpha }_{1}}\ge 0$ such that for every $\varepsilon > 0$ and for every $\varphi \in (L\left[0, T\right], R)$ satisfies the following inequality ${{}^{FFM}J}_{0, t}^{\mathsf{μ}, {\alpha }_{1}}\left(\varphi \left(t\right)\right)-£ \left(t, \varphi \left(t\right)\right)\le \varepsilon, t\in \left[0, T\right]$ such that $\left|\varphi \left(t\right)-£ \left(t\right)\right|\le {B}_{\mathsf{μ}, {\alpha }_{1}}\varepsilon, t\in \left[0, T\right]$ .

Suppose a perturbation $\omega \in L\left[0, T\right], R$ then $\omega \left(0\right) = 0$ and

● For every $\varepsilon > 0 \exists \omega \left(t\right)\le \varepsilon |$

● ${{}_{0}{}^{FFM}J}_{t}^{\mathsf{μ}, {\alpha }_{1}}\left(\varphi \left(t\right)\right) = £ \left(t, \varphi \left(t\right)\right)+\omega \left(t\right).$

Lemma 4.4: The solution of the perturbed model ${{}_{0}{}^{FFM}J}_{t}^{\mathsf{μ}, {\alpha }_{1}}\left(\varphi \left(t\right)\right) = £ \left(t, \varphi \left(t\right)\right)+\omega \left(t\right), \varphi \left(0\right) = {\varphi }_{0}$ fulfills the relation

$\begin{array}{l} B\left(t\right)-[\varphi \left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}(1-{\alpha }_{1})}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, \varphi \left(t\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(\lambda , \varphi \left(\lambda \right)\right)d\lambda \\ \ \ \ \ \ \ \ \ \le {{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}\varepsilon \end{array}$

Where ${{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}\varepsilon = \frac{\mathsf{μ} {T}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{T}^{\mathsf{μ} +{\alpha }_{1}-1}M\left(\mathsf{μ}, {\alpha }_{1}\right).$

Lemma 4.5: By using condition (12) with lemma (4.4), proposed model is Ulam-Hyres stable if $\eta < 1.$

Proof: Suppose ${\alpha }_{1}\in Z$ be a solution and $\varphi \in Z$ be any solution of (1), then

$\begin{array}{l} |\varphi \left(t\right)-{\alpha }_{1}\left(t\right)| = |\varphi \left(t\right)-[{\alpha }_{1}\left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, {\alpha }_{1}\left(t\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t}{\lambda }^{\mathsf{μ} -1}{(1-}\\ \lambda )^{\mathsf{μ} -1}£ \left(\lambda , {\alpha }_{1}\left(\lambda \right)\right)d\lambda ]| \end{array}$

$\begin{array}{l} \le \left|\varphi \left(t\right)-\left[\varphi \left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, \varphi \left(t\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(\lambda , \varphi \left(\lambda \right)\right)d\lambda \right]\right|\\ \ \ \ \ \ \ \ \ \ \ \ \ +|\varphi \left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, \varphi \left(t\right)\right)\\ \ \ \ \ \ \ \ \ \ \ \ \ +\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\underset{0}{\overset{t}{\int }}{\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(\lambda , \varphi \left(\lambda \right)\right)d\lambda | \end{array}$

$- |{\alpha }_{1}\left(0\right)+\frac{\mathsf{μ} {t}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}£ \left(t, {\alpha }_{1}\left(t\right)\right)+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\lambda }^{\mathsf{μ} -1}{\left(1-\lambda \right)}^{\mathsf{μ} -1}£ \left(\lambda , {\alpha }_{1}\left(\lambda \right)\right)d\lambda |$

$\le {{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}\varepsilon +(\frac{\mathsf{μ} {T}^{\mathsf{μ} -1}\left(1-{\alpha }_{1}\right)}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)}+\frac{\mathsf{μ} {\alpha }_{1}}{\mathrm{A}\mathrm{B}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}{T}^{\mathsf{μ} +{\alpha }_{1}-1}){L}_{£ }|\varphi \left(t\right)-{\alpha }_{1}\left(t\right)|$

$\le {{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}\varepsilon +\eta \left|\varphi \left(t\right)-{\alpha }_{1}\left(t\right)\right|.$

Consequently,

$\left|\right|\varphi -{\alpha }_{1}\left|\right|\le {{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}\varepsilon +\eta \left|\left|\varphi \left(t\right)-{\alpha }_{1}\left(t\right)\right|\right|.$

So, we can write it as

$\left|\left|\varphi -{\alpha }_{1}\right|\right|\le {B}_{{\alpha }_{1}, \mathsf{μ} }\varepsilon ,$

Where ${B}_{{\alpha }_{1}, \mathsf{μ} }\varepsilon = \frac{{{\alpha }_{1}}_{{\alpha }_{1}, \mathsf{μ} }^{*}}{1-\eta }.$ Hence the solution is Ulam-Hyres stable.

5. Fractal-fractional operator for hires problem

In this section, we present the Hires problem model $\left(1\right)$ using fractal-fractional Atangana-Baleanu derivative. We have

$\begin{array}{l} {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{1} = -{M}_{1}{y}_{1}+{M}_{2}{y}_{2}+{M}_{3}{y}_{3}+{M}_{4}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{2} = {M}_{1}{y}_{1}-{M}_{5}{y}_{2}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{3} = -{M}_{6}{y}_{3}+{M}_{2}{y}_{4}+{M}_{7}{y}_{5}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{4} = {M}_{3}{y}_{2}+{M}_{8}{y}_{3}-{M}_{9}{y}_{4}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{5} = -{M}_{10}{y}_{5}+{M}_{2}{y}_{6}-{M}_{2}{y}_{7}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{6} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{12}{y}_{4}+{M}_{8}{y}_{5}-{M}_{2}{y}_{6}+{M}_{12}{y}_{7}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{7} = {M}_{11}{y}_{6}{y}_{8}-{M}_{13}{y}_{7}, \\ {}^{FF}{D}_{0, t}^{{\alpha }_{1}, {\alpha }_{2}}{y}_{8} = -{M}_{11}{y}_{6}{y}_{8}+{M}_{13}{y}_{7}. \end{array}$

(13)

With initial conditions

$\begin{array}{l} {y}_{1}\left(0\right) = {{y}_{1}}_{\left(0\right)}, {y}_{2}\left(0\right) = {{y}_{2}}_{\left(0\right)}, {y}_{3}\left(0\right) = {{y}_{3}}_{\left(0\right)}, {y}_{4}\left(0\right) = {{y}_{4}}_{\left(0\right)}, \\ {y}_{5}\left(0\right) = {{y}_{5}}_{\left(0\right)}, {y}_{6}\left(0\right) = {{y}_{6}}_{\left(0\right)}, {y}_{7}\left(0\right) = {{y}_{7}}_{\left(0\right)}, {y}_{8}\left(0\right) = {{y}_{8}}_{\left(0\right)}. \end{array}$

We present the numerical algorithm for the fractal-fractional Hires problem model (13). The following is obtained by integrating the system (13).

$\begin{array}{l} {y}_{1}\left(t\right)-{y}_{1}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{1}{y}_{1}\left(t\right)+{M}_{2}{y}_{2}\left(t\right)+{M}_{3}{y}_{3}\left(t\right)+{M}_{4}\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\left\{-{M}_{1}{y}_{1}\left(\tau \right)+{M}_{2}{y}_{2}\left(\tau \right)+{M}_{3}{y}_{3}\left(\tau \right)+{M}_{4}\right\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{2}\left(t\right)-{y}_{2}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{1}{y}_{1}\left(t\right)-{M}_{5}{y}_{2}\left(t\right)\right\}+\frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\{\{{M}_{1}{y}_{1}\left(\tau \right)-\\ {M}_{5}{y}_{2}\left(\tau \right)\}\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{3}\left(t\right)-{y}_{3}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{6}{y}_{3}\left(t\right)+{M}_{2}{y}_{4}\left(t\right)+{M}_{7}{y}_{5}\left(t\right)\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\left\{-{M}_{6}{y}_{3}\left(\tau \right)+{M}_{2}{y}_{4}\left(\tau \right)+{M}_{7}{y}_{5}\left(\tau \right)\right\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{4}\left(t\right)-{y}_{4}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{3}{y}_{2}\left(t\right)+{M}_{8}{y}_{3}\left(t\right)-{M}_{9}{y}_{4}\left(t\right)\right\}+\frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\{{M}_{3}{y}_{2}\left(\tau \right)+\\ {M}_{8}{y}_{3}\left(\tau \right)-{M}_{9}{y}_{4}\left(\tau \right)\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{5}\left(t\right)-{y}_{5}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{10}{y}_{5}\left(t\right)+{M}_{2}{y}_{6}\left(t\right)-{M}_{2}{y}_{7}\left(t\right)\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\left\{\left\{-{M}_{10}{y}_{5}\left(\tau \right)+{M}_{2}{y}_{6}\left(\tau \right)-{M}_{2}{y}_{7}\left(\tau \right)\right\}\right\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{6}\left(t\right)-{y}_{6}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{-{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)+{M}_{12}{y}_{4}\left(t\right)+{M}_{8}{y}_{5}\left(t\right)-{M}_{2}{y}_{6}\left(t\right)+\\ {M}_{12}{y}_{7}\left(t\right)\}+\frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\{-{M}_{11}{y}_{6}\left(\tau \right){y}_{8}\left(\tau \right)+{M}_{12}{y}_{4}\left(\tau \right)+{M}_{8}{y}_{5}\left(\tau \right)-{M}_{2}{y}_{6}\left(\tau \right)+\\ {M}_{12}{y}_{7}\left(\tau \right)\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{7}\left(t\right)-{y}_{7}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)-{M}_{13}{y}_{7}\left(t\right)\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\left\{{M}_{11}{y}_{6}\left(\tau \right){y}_{8}\left(\tau \right)-{M}_{13}{y}_{7}\left(\tau \right)\right\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{8}\left(t\right)-{y}_{8}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)+{M}_{13}{y}_{7}\left(t\right)\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} {\tau }^{{\alpha }_{2}-1}\left\{\left\{-{M}_{11}{y}_{6}\left(\tau \right){y}_{8}\left(\tau \right)+{M}_{13}{y}_{7}\left(\tau \right)\right\}\right\}{\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \end{array}$

(14)

Let

$k\left(t, {y}_{1}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{1}{y}_{1}\left(t\right)+{M}_{2}{y}_{2}\left(t\right)+{M}_{3}{y}_{3}\left(t\right)+{M}_{4}\right\},$

$k\left(t, {y}_{2}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{1}{y}_{1}\left(t\right)-{M}_{5}{y}_{2}\left(t\right)\right\},$

$k\left(t, {y}_{3}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{6}{y}_{3}\left(t\right)+{M}_{2}{y}_{4}\left(t\right)+{M}_{7}{y}_{5}\left(t\right)\right\},$

$k\left(t, {y}_{4}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{3}{y}_{2}\left(t\right)+{M}_{8}{y}_{3}\left(t\right)-{M}_{9}{y}_{4}\left(t\right)\right\},$

$k\left(t, {y}_{5}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{10}{y}_{5}\left(t\right)+{M}_{2}{y}_{6}\left(t\right)-{M}_{2}{y}_{7}\left(t\right)\right\},$

$k\left(t, {y}_{6}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)+{M}_{12}{y}_{4}\left(t\right)+{M}_{8}{y}_{5}\left(t\right)-{M}_{2}{y}_{6}\left(t\right)+{M}_{12}{y}_{7}\left(t\right)\right\},$

$k\left(t, {y}_{7}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)-{M}_{13}{y}_{7}\left(t\right)\right\},$

$k\left(t, {y}_{8}\left(t\right)\right) = {\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)+{M}_{13}{y}_{7}\left(t\right)\right\}.$

Then system (14) becomes

$\begin{array}{l} {y}_{1}\left(t\right)-{y}_{1}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{1}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{1}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{2}\left(t\right)-{y}_{2}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{2}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{2}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{3}\left(t\right)-{y}_{3}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{3}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{3}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{4}\left(t\right)-{y}_{4}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{4}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{4}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{5}\left(t\right)-{y}_{5}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{5}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{5}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{6}\left(t\right)-{y}_{6}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{6}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{6}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{7}\left(t\right)-{y}_{7}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{7}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{7}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{8}\left(t\right)-{y}_{8}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left(t, {y}_{8}\left(t\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{t} k\left(\tau , {y}_{8}\left(\tau \right)\right){\left(t-\tau \right)}^{{\alpha }_{1}-1}d\tau , \end{array}$

(15)

At ${t}_{n+1} = \left(n+1\right)\Delta t,$ we have

$\begin{array}{l} {y}_{1}\left({t}_{n+1}\right)-{y}_{1}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{1}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{1}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{2}\left({t}_{n+1}\right)-{y}_{2}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{2}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{2}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{3}\left({t}_{n+1}\right)-{y}_{3}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{3}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{3}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{4}\left({t}_{n+1}\right)-{y}_{4}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{4}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{4}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{5}\left({t}_{n+1}\right)-{y}_{5}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{5}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{5}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{6}\left({t}_{n+1}\right)-{y}_{6}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{6}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{6}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{7}\left({t}_{n+1}\right)-{y}_{7}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{7}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{7}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{8}\left({t}_{n+1}\right)-{y}_{8}\left(0\right) = \frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{8}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\int_{0}^{{t}_{n+1}} k\left(\tau , {y}_{8}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau . \end{array}$

(16)

Also, we have

$\begin{array}{l} {y}_{1}\left({t}_{n+1}\right) = {y}_{1}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{1}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{1}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{2}\left({t}_{n+1}\right) = {y}_{2}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{2}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{2}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{3}\left({t}_{n+1}\right) = {y}_{3}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{3}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{3}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{4}\left({t}_{n+1}\right) = {y}_{4}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{4}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{4}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{5}\left({t}_{n+1}\right) = {y}_{5}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{5}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{5}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{6}\left({t}_{n+1}\right) = {y}_{6}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{6}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{6}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{7}\left({t}_{n+1}\right) = {y}_{7}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{7}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{7}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \\ {y}_{8}\left({t}_{n+1}\right) = {y}_{8}\left(0\right)+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{8}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 0}^{n}\int_{{t}_{j}}^{{t}_{j+1}} k\left(\tau , {y}_{8}\left(\tau \right)\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau . \end{array}$

(17)

In general, approximating the function $k\left(\tau, y\left(\tau \right)\right),$ using the Newton polynomial, we have

$\begin{array}{l} {P}_{n}\left(\tau \right) = \frac{k\left({t}_{n}, y\left({t}_{n}\right)\right)}{{t}_{n}-{t}_{n-1}}\left(\tau -{t}_{n-1}\right)+\frac{k\left({t}_{n-1}, y\left({t}_{n-1}\right)\right)}{{t}_{n}-{t}_{n-1}}\left(\tau -{t}_{n}\right) = \frac{k\left({t}_{n}, y\left({t}_{n}\right)\right)}{h}\left(\tau -{t}_{n-1}\right)\frac{k\left({t}_{n-1}, y\left({t}_{n-1}\right)\right)}{h}\left(\tau -\\ {t}_{n}\right). \end{array}$

(18)

Using Eq (18) into system (17) we have

(19)

Rearranging the above equation, we have

(20)

Writing further system (20) we have

$\begin{array}{l} {{y}_{1}}^{n+1} = {{y}_{1}}^{0}+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{1}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 2}^{n}k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)\int_{{t}_{j}}^{{t}_{j+1}} {\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau +\\ \frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 2}^{n}\frac{k\left({t}_{j-1}, {{y}_{1}}^{j-1}\right)-k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)}{\Delta t}\int_{{t}_{j}}^{{t}_{j+1}} \left(\tau -{t}_{j-2}\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau +\\ \frac{{\alpha }_{1}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}\right)}\sum _{j = 2}^{n}\frac{k\left({t}_{j}, {{y}_{1}}^{j}\right)-2k\left({t}_{j-1}, {{y}_{1}}^{j-1}\right)+k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)}{{2\left(\Delta t\right)}^{2}}\int_{{t}_{j}}^{{t}_{j+1}} \left(\tau -{t}_{j-2}\right)\left(\tau -{t}_{j-1}\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau , \end{array}$

(21)

Now, calculating the integrals in system (21) we get

$\int_{{t}_{j}}^{{t}_{j+1}} {\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau = \frac{{\left(\Delta t\right)}^{{\alpha }_{1}}}{{\alpha }_{1}}\left[{\left(n-j+1\right)}^{{\alpha }_{1}}-{\left(n-j\right)}^{{\alpha }_{1}}\right],$

$\begin{array}{l} \int_{{t}_{j}}^{{t}_{j+1}} \left(\tau -{t}_{j-2}\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau = \frac{{\left(\Delta t\right)}^{{\alpha }_{1}+1}}{{\alpha }_{1}\left({\alpha }_{1}+1\right)}[{\left(n-j+1\right)}^{{\alpha }_{1}}\left(n-j+3+2{\alpha }_{1}\right)-{\left(n-j+1\right)}^{{\alpha }_{1}}(n-\\ j+3+3{\alpha }_{1})], \end{array}$

$\begin{array}{l} \int_{{t}_{j}}^{{t}_{j+1}} \left(\tau -{t}_{j-2}\right)\left(\tau -{t}_{j-1}\right){\left({t}_{n+1}-\tau \right)}^{{\alpha }_{1}-1}d\tau = \frac{{\left(\Delta t\right)}^{{\alpha }_{1}+2}}{{\alpha }_{1}\left({\alpha }_{1}+1\right)\left({\alpha }_{1}+2\right)}[{\left(n-j+1\right)}^{{\alpha }_{1}}\left\{2{\left(n-j\right)}^{2}+\\ \left(3{\alpha }_{1}+10\right)\left(n-j\right)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\right\}-{\left(n-j\right)}^{{\alpha }_{1}}\{2{\left(n-j\right)}^{2}+\left(5{\alpha }_{1}+10\right)\left(n-j\right)+6{{\alpha }_{1}}^{2}+\\ 18{\alpha }_{1}+12\}]. \end{array}$

Inserting them into system (21) we get

$\begin{array}{l} {{y}_{1}}^{n+1} = {{y}_{1}}^{0}+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}k\left({t}_{n}, {y}_{1}\left({t}_{n}\right)\right)+\frac{{\alpha }_{1}{\left(\Delta t\right)}^{{\alpha }_{1}}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+1\right)}\sum _{j = 2}^{n}k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)\left[{\left(n-j+1\right)}^{{\alpha }_{1}}-{\left(n-j\right)}^{{\alpha }_{1}}\right]+\\ \frac{{\alpha }_{1}{\left(\Delta t\right)}^{{\alpha }_{1}}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+2\right)}\sum _{j = 2}^{n}\left[k\left({t}_{j-1}, {{y}_{1}}^{j-1}\right)-k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)\right][{(n-j+1)}^{{\alpha }_{1}}\left(n-j+3+2{\alpha }_{1}\right)-(n-j+\\ 1)^{{\alpha }_{1}}\left(n-j+3+3{\alpha }_{1}\right)]+\frac{{\alpha }_{1}{\left(\Delta t\right)}^{{\alpha }_{1}}}{2\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+3\right)}\sum _{j = 2}^{n}\left[k\left({t}_{j}, {{y}_{1}}^{j}\right)-2k\left({t}_{j-1}, {{y}_{1}}^{j-1}\right)+k\left({t}_{j-2}, {{y}_{1}}^{j-2}\right)\right][(n-\\ j+1)^{{\alpha }_{1}}\left\{2{\left(n-j\right)}^{2}+\left(3{\alpha }_{1}+10\right)\left(n-j\right)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\right\}-{\left(n-j\right)}^{{\alpha }_{1}}\{2{\left(n-j\right)}^{2}+(5{\alpha }_{1}+\\ 10)\left(n-j\right)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}] \end{array}$

(22)

Finally, we have the following approximation:

$\begin{array}{l} {{y}_{1}}^{n+1} = {{y}_{1}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{-{M}_{1}{y}_{1}(t)+{M}_{2}{y}_{2}(t)+{M}_{3}{y}_{3}(t)+{M}_{4}\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j-2}+{M}_{2}{{y}_{2}}^{j-2}+{M}_{3}{{y}_{3}}^{j-2}+{M}_{4}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j-1}+{M}_{2}{{y}_{2}}^{j-1}+{M}_{3}{{y}_{3}}^{j-1}+{M}_{4}\}-{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j-2}+{M}_{2}{{y}_{2}}^{j-2}+\\ {M}_{3}{{y}_{3}}^{j-2}+{M}_{4}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j}+{M}_{2}{{y}_{2}}^{j}+{M}_{3}{{y}_{3}}^{j}+{M}_{4}\}-2{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j-1}+{M}_{2}{{y}_{2}}^{j-1}+\\ {M}_{3}{{y}_{3}}^{j-1}+{M}_{4}\}+{t}^{{\alpha }_{2}-1}\{-{M}_{1}{{y}_{1}}^{j-2}+{M}_{2}{{y}_{2}}^{j-2}+{M}_{3}{{y}_{3}}^{j-2}+{M}_{4}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+\\ (3{\alpha }_{1}+10)(n-j)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+6{{\alpha }_{1}}^{2}+\\ 18{\alpha }_{1}+12\}], \\ {{y}_{2}}^{n+1} = {{y}_{2}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{{M}_{1}{y}_{1}(t)-{M}_{5}{y}_{2}(t)\}+\frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j-2}-\\ {M}_{5}{{y}_{2}}^{j-2}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j-1}-{M}_{5}{{y}_{2}}^{j-1}\}-\\ {t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j-2}-{M}_{5}{{y}_{2}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+\\ 3{\alpha }_{1})]+\frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j}-{M}_{5}{{y}_{2}}^{j}\}-2{t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j-1}-{M}_{5}{{y}_{2}}^{j-1}\}+\\ {t}^{{\alpha }_{2}-1}\{{M}_{1}{{y}_{1}}^{j-2}-{M}_{5}{{y}_{2}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(3{\alpha }_{1}+10)(n-j)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\}-\\ {(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}], \\ {{y}_{3}}^{n+1} = {{y}_{3}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{-{M}_{6}{y}_{3}(t)+{M}_{2}{y}_{4}(t)+{M}_{7}{y}_{5}(t)\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j-2}+{M}_{2}{{y}_{4}}^{j-2}+{M}_{7}{{y}_{5}}^{j-2}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j-1}+{M}_{2}{{y}_{4}}^{j-1}+{M}_{7}{{y}_{5}}^{j-1}\}-{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j-2}+{M}_{2}{{y}_{4}}^{j-2}+\\ {M}_{7}{{y}_{5}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j}+{M}_{2}{{y}_{4}}^{j}+\\ {M}_{7}{{y}_{5}}^{j}\}-2{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j-1}+{M}_{2}{{y}_{4}}^{j-1}+{M}_{7}{{y}_{5}}^{j-1}\}+{t}^{{\alpha }_{2}-1}\{-{M}_{6}{{y}_{3}}^{j-2}+{M}_{2}{{y}_{4}}^{j-2}+\\ {M}_{7}{{y}_{5}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(3{\alpha }_{1}+10)(n-\\ j)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}], \\ {{y}_{4}}^{n+1} = {{y}_{4}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{{M}_{3}{y}_{2}(t)+{M}_{8}{y}_{3}(t)-{M}_{9}{y}_{4}(t)\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j-2}+{M}_{8}{{y}_{3}}^{j-2}-{M}_{9}{{y}_{4}}^{j-2}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j-1}+{M}_{8}{{y}_{3}}^{j-1}-{M}_{9}{{y}_{4}}^{j-1}\}-{t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j-2}+{M}_{8}{{y}_{3}}^{j-2}-\\ {M}_{9}{{y}_{4}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j}+{M}_{8}{{y}_{3}}^{j}-{M}_{9}{{y}_{4}}^{j}\}-2{t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j-1}+{M}_{8}{{y}_{3}}^{j-1}-{M}_{9}{{y}_{4}}^{j-1}\}+\\ {t}^{{\alpha }_{2}-1}\{{M}_{3}{{y}_{2}}^{j-2}+{M}_{8}{{y}_{3}}^{j-2}-{M}_{9}{{y}_{4}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(3{\alpha }_{1}+10)(n-j)+2{{\alpha }_{1}}^{2}+\\ 9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}], \\ {{y}_{5}}^{n+1} = {{y}_{5}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{-{M}_{10}{y}_{5}(t)+{M}_{2}{y}_{6}(t)-{M}_{2}{y}_{7}(t)\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j-2}+{M}_{2}{{y}_{6}}^{j-2}-{M}_{2}{{y}_{7}}^{j-2}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j-1}+{M}_{2}{{y}_{6}}^{j-1}-{M}_{2}{{y}_{7}}^{j-1}\}-{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j-2}+{M}_{2}{{y}_{6}}^{j-2}-\\ {M}_{2}{{y}_{7}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j}+{M}_{2}{{y}_{6}}^{j}-\\ {M}_{2}{{y}_{7}}^{j}\}-2{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j-1}+{M}_{2}{{y}_{6}}^{j-1}-{M}_{2}{{y}_{7}}^{j-1}\}+{t}^{{\alpha }_{2}-1}\{-{M}_{10}{{y}_{5}}^{j-2}+\\ {M}_{2}{{y}_{6}}^{j-2}-{M}_{2}{{y}_{7}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(3{\alpha }_{1}+\\ 10)(n-j)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+\\ 10)(n-j)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+\\ 12\}], \end{array}$

(23)

$\begin{array}{l} {{y}_{6}}^{n+1} = {{y}_{6}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{-{M}_{11}{y}_{6}(t){y}_{8}(t)+{M}_{12}{y}_{4}(t)+{M}_{8}{y}_{5}(t)-{M}_{2}{y}_{6}(t)+{M}_{12}{y}_{7}(t)\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+{M}_{12}{{y}_{4}}^{j-2}+{M}_{8}{{y}_{5}}^{j-2}-{M}_{2}{{y}_{6}}^{j-2}+{M}_{12}{{y}_{7}}^{j-2}\}[(n-\\ j+1)^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}+{M}_{12}{{y}_{4}}^{j-1}+{M}_{8}{{y}_{5}}^{j-1}-\\ {M}_{2}{{y}_{6}}^{j-1}+{M}_{12}{{y}_{7}}^{j-1}\}-{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+{M}_{12}{{y}_{4}}^{j-2}+\\ {M}_{8}{{y}_{5}}^{j-2}-{M}_{2}{{y}_{6}}^{j-2}+{M}_{12}{{y}_{7}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j}{{y}_{8}}^{j}+{M}_{12}{{y}_{4}}^{j}+{M}_{8}{{y}_{5}}^{j}-{M}_{2}{{y}_{6}}^{j}+{M}_{12}{{y}_{7}}^{j}\}-\\ 2{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}+{M}_{12}{{y}_{4}}^{j-1}+{M}_{8}{{y}_{5}}^{j-1}-{M}_{2}{{y}_{6}}^{j-1}+{M}_{12}{{y}_{7}}^{j-1}\}+\\ {t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+{M}_{12}{{y}_{4}}^{j-2}+{M}_{8}{{y}_{5}}^{j-2}-{M}_{2}{{y}_{6}}^{j-2}+{M}_{12}{{y}_{7}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2(n-\\ j)^{2}+(3{\alpha }_{1}+10)(n-j)+2{{\alpha }_{1}}^{2}+9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+\\ 6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}], \end{array}$

$\begin{array}{l} {{y}_{7}}^{n+1} = {{y}_{7}}^{0}+\frac{(1-{\alpha }_{1})}{\mathrm{C}({\alpha }_{1})}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\{{M}_{11}{y}_{6}(t){y}_{8}(t)-{M}_{13}{y}_{7}(t)\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+1)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}-{M}_{13}{{y}_{7}}^{j-2}\}[{(n-j+1)}^{{\alpha }_{1}}-{(n-j)}^{{\alpha }_{1}}]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+2)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}-{M}_{13}{{y}_{7}}^{j-1}\}-{t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}-{M}_{13}{{y}_{7}}^{j-2}\}][(n-\\ j+1)^{{\alpha }_{1}}(n-j+3+2{\alpha }_{1})-{(n-j+1)}^{{\alpha }_{1}}(n-j+3+3{\alpha }_{1})]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{(\Delta t)}^{{\alpha }_{1}}}{2\mathrm{C}({\alpha }_{1})\mathrm{\Gamma }({\alpha }_{1}+3)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j}{{y}_{8}}^{j}-{M}_{13}{{y}_{7}}^{j}\}-2{t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}-{M}_{13}{{y}_{7}}^{j-1}\}+\\ {t}^{{\alpha }_{2}-1}\{{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}-{M}_{13}{{y}_{7}}^{j-2}\}][{(n-j+1)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(3{\alpha }_{1}+10)(n-j)+2{{\alpha }_{1}}^{2}+\\ 9{\alpha }_{1}+12\}-{(n-j)}^{{\alpha }_{1}}\{2{(n-j)}^{2}+(5{\alpha }_{1}+10)(n-j)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\}], \end{array}$

$\begin{array}{l} {{y}_{8}}^{n+1} = {{y}_{8}}^{0}+\frac{\left(1-{\alpha }_{1}\right)}{\mathrm{C}\left({\alpha }_{1}\right)}{\alpha }_{2}{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{y}_{6}\left(t\right){y}_{8}\left(t\right)+{M}_{13}{y}_{7}\left(t\right)\right\}+\\ \frac{{\alpha }_{1}{\alpha }_{2}{\left(\Delta t\right)}^{{\alpha }_{1}}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+1\right)}\sum _{j = 2}^{n}{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+{M}_{13}{{y}_{7}}^{j-2}\right\}\left[{\left(n-j+1\right)}^{{\alpha }_{1}}-{\left(n-j\right)}^{{\alpha }_{1}}\right]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{\left(\Delta t\right)}^{{\alpha }_{1}}}{\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+2\right)}\sum _{j = 2}^{n}[{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}+{M}_{13}{{y}_{7}}^{j-1}\right\}-{t}^{{\alpha }_{2}-1}\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+\\ {M}_{13}{{y}_{7}}^{j-2}\}]\left[{\left(n-j+1\right)}^{{\alpha }_{1}}\left(n-j+3+2{\alpha }_{1}\right)-{\left(n-j+1\right)}^{{\alpha }_{1}}\left(n-j+3+3{\alpha }_{1}\right)\right]+\\ \frac{{\alpha }_{1}{\alpha }_{2}{\left(\Delta t\right)}^{{\alpha }_{1}}}{2\mathrm{C}\left({\alpha }_{1}\right)\mathrm{\Gamma }\left({\alpha }_{1}+3\right)}\sum _{j = 2}^{n}\left[{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{{y}_{6}}^{j}{{y}_{8}}^{j}+{M}_{13}{{y}_{7}}^{j}\right\}-2{t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{{y}_{6}}^{j-1}{{y}_{8}}^{j-1}+{M}_{13}{{y}_{7}}^{j-1}\right\}+\\ {t}^{{\alpha }_{2}-1}\left\{-{M}_{11}{{y}_{6}}^{j-2}{{y}_{8}}^{j-2}+{M}_{13}{{y}_{7}}^{j-2}\right\}\right][{\left(n-j+1\right)}^{{\alpha }_{1}}\{2{\left(n-j\right)}^{2}+\left(3{\alpha }_{1}+10\right)\left(n-j\right)+2{{\alpha }_{1}}^{2}+\\ 9{\alpha }_{1}+12\}-{\left(n-j\right)}^{{\alpha }_{1}}\left\{2{\left(n-j\right)}^{2}+\left(5{\alpha }_{1}+10\right)\left(n-j\right)+6{{\alpha }_{1}}^{2}+18{\alpha }_{1}+12\right\}]. \end{array}$

6. Results and discussions

A fractional-order model is proposed for analysis and simulation, to observe the concentration of chemicals in chemistry kinematics problems with a stiff differential equation. For this purpose, we used ABC with Mittage-Lefffier law, Atangana-Tufik scheme, and fractal fractional derivative for hires problem with given initial conditions. Details of parameters values of real data are also given in ^[18,19] which will consider for simulation analysis for the proposed study. Solution of compartment shows in to with fractional fractal operator at different order. Effect of fraction order can easily be observed in simulation of the compartments having a concentration of chemical reaction with stiff differential equations. The concentration ${y}_{1}$ and ${y}_{8}$ of the chemical species start decreasing by decreasing fractional values respectively while concentration ${y}_{2}, {y}_{3}$ , ${y}_{4}$ , ${y}_{5}$ , ${y}_{6}$ and ${y}_{7}$ of the chemical species start increasing by decreasing fractional values. These concentrations of chemical species converge to our desired value according to steady state by decreasing the fractional values which shows that it provides us appropriate results at non integer value. We can get better concentration of the components by using the fractional derivative which are very important for chemical problem to check the actual behavior of the concentration of the chemical with smallest changes in derivative with respect to time. It is also very important for solutions of nonlinear problems which are commonly used researcher and scientist in kinetics chemistry.

Figure 1. Simulation of

${y}_{1}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 2. Simulation of

${y}_{2}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 3. Simulation of

${y}_{3}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 4. Simulation of

${y}_{4}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 5. Simulation of

${y}_{5}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 6. Simulation of

${y}_{6}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 7. Simulation of

${y}_{7}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

Figure 8. Simulation of

${y}_{8}\left(t\right)$ with fractal fractional derivative.

DownLoad: Full-Size Img PowerPoint

7. Conclusions

We examine the hires problems with stiff systems of nonlinear ordinary equations that rely on the concentration of chemical reaction of components in this study. The advanced techniques of fractional operator have been implemented for initial value problem arising from chemical reactions composed of large systems of stiff ordinary differential equations. The arbitrary derivative of fractional order has been taken with Atangana-Toufik scheme and fractal fractional derivative. Solutions have been obtained efficiently within limited time which shows the actual behavior of kinetic chemical reactions. Existence and uniqueness of results have been verified by fixed point theorem. Simulations are carried out for different fractional values. New chemical reactions can be done with the help of these analyses. These concepts are very important to use for real life problems like Brine tank cascade, Recycled Brine tank cascade, pond pollution, home heating and biomass transfer problem.

Acknowledgments

Research Supporting Project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia.

Conflict of interest

No conflict of interest.

References

[1]	H. I. Abdel-Gawad, M. Tantawy, D. Baleanu, Fractional KdV and Boussenisq-Burger's equations, reduction to PDE and stability approaches, Math. Method. Appl. Sci., 43 (2020), 4125–4135. http://doi.org/10.1002/mma.6178 doi: 10.1002/mma.6178
[2]	H. I. Abdel-Gawad, M. Tantawy, Traveling wave solutions of DNA-Torsional model of fractional order, Appl. Math. Inf. Sci. Lett., 6 (2018), 85–89. http://doi.org/10.18576/amisl/060205 doi: 10.18576/amisl/060205
[3]	N. I. Okposo, P. Veeresha, E. N. Okposo, Solutions for time-fractional coupled nonlinear Schrödinger equations arising in optical solitons, Chinese J. Phy., 77 (2022), 965–984. http://doi.org/10.1016/j.cjph.2021.10.014 doi: 10.1016/j.cjph.2021.10.014
[4]	P. Veeresha, M. Yavuz, C. Baishya, A computational approach for shallow water forced Korteweg-De Vries equation on critical flow over a hole with three fractional operators, IJOCTA, 11 (2021), 52–67. http://doi.org/10.11121/ijocta.2021.1177 doi: 10.11121/ijocta.2021.1177
[5]	P. Veeresha, E. Ilhan, H. M. Baskonus, Fractional approach for analysis of the model describing wind-influenced projectile motion, Phy. Scr., 96 (2021), 075209. http://doi.org/10.1088/1402-4896/abf868 doi: 10.1088/1402-4896/abf868
[6]	P. Veeresha, D. Baleanu, A unifying computational framework for fractional Gross-Pitaevskii equations, Phy. Scr., 96 (2021), 125010. http://doi.org/10.1088/1402-4896/ac28c9 doi: 10.1088/1402-4896/ac28c9
[7]	P. Veeresha, D. G. Prakasha, H. M. Baskonus, New numerical surfaces to the mathematical model of cancer chemotherapy effect in Caputo fractional derivatives, Chaos, 29 (2019), 013119. http://doi.org/10.1063/1.5074099 doi: 10.1063/1.5074099
[8]	C. Baishya, P. Veeresha, Laguerre polynomial-based operational matrix of integration for solving fractional differential equations with non-singular kernel, Proc. R. Soc. A, 477 (2021), 20210438. http://doi.org/10.1098/rspa.2021.0438 doi: 10.1098/rspa.2021.0438
[9]	S. N. Hajiseyedazizi, M. E. Samei, J. Alzabut, Y.-M. Chu, On multi-step methods for singular fractional $q$ -integro-differential equations, Open Math., 19 (2021), 1378–1405. http://doi.org/10.1515/math-2021-0093 doi: 10.1515/math-2021-0093
[10]	F. Jin, Z.-S. Qian, Y.-M. Chu, M. ur Rahman, On nonlinear evolution model for drinking behavior under Caputo-Fabrizio derivative, J. Appl. Anal. Comput., 12 (2022), 790–806. http://doi.org/10.11948/20210357 doi: 10.11948/20210357
[11]	S. Rashid, E. I. Abouelmagd, A. Khalid, F. B. Farooq, Y.-M. Chu, Some recent developments on dynamical $\hbar$ -discrete fractional type inequalities in the frame of nonsingular and nonlocal kernels, Fractals, 30 (2022), 2240110. http://doi.org/10.1142/S0218348X22401107 doi: 10.1142/S0218348X22401107
[12]	F.-Z. Wang, M. N. Khan, I. Ahmad, H. Ahmad, H. Abu-Zinadah, Y.-M. Chu, Numerical solution of traveling waves in chemical kinetics: time-fractional Fishers equations, Fractals, 30 (2022), 2240051. http://doi.org/10.1142/S0218348X22400515 doi: 10.1142/S0218348X22400515
[13]	S. Rashid, E. I. Abouelmagd, S. Sultana, Y.-M. Chu, New developments in weighted $n$ -fold type inequalities via discrete generalized h-proportional fractional operators, Fractals, 30 (2022), 2240056. http://doi.org/10.1142/S0218348X22400564 doi: 10.1142/S0218348X22400564
[14]	Y.-M. Chu, S. Bashir, M. Ramzan, M. Y. Malik, Model-based comparative study of magnetohydrodynamics unsteady hybrid nanofluid flow between two infinite parallel plates with particle shape effects, Math. Method. Appl. Sci., in press. http://doi.org/10.1002/mma.8234
[15]	S. A. Iqbal, M. G. Hafez, Y.-M. Chu, C. Park, Dynamical Analysis of nonautonomous RLC circuit with the absence and presence of Atangana-Baleanu fractional derivative, J. Appl. Anal. Comput., 12 (2022), 770–789. http://doi.org/10.11948/20210324 doi: 10.11948/20210324
[16]	C. Baishya1, S. J. Achar, P. Veeresha, D. G. Prakasha, Dynamics of a fractional epidemiological model with disease infection in both the populations, Chaos, 31 (2021), 043130. http://doi.org/10.1063/5.0028905 doi: 10.1063/5.0028905
[17]	E. Ilhan, P. Veeresha, H. M. Baskonus, Fractional approach for a mathematical model of atmospheric dynamics of CO2 gas with an efficient method, Chaos Soliton. Fract., 152 (2021), 111347. http://doi.org/10.1016/j.chaos.2021.111347 doi: 10.1016/j.chaos.2021.111347
[18]	P. Veeresha, H. M. Baskonus, W. Gao, Strong interacting internal waves in rotating ocean: Novel fractional approach, Axioms, 10 (2021), 123. http://doi.org/10.3390/axioms10020123 doi: 10.3390/axioms10020123
[19]	P. Veeresha, D. G. Prakasha, D. Kumar, An efficient technique for nonlinear time-fractional Klein–Fock–Gordon equation, Appl. Math. Comput., 364 (2020), 124637. http://doi.org/10.1016/j.amc.2019.124637 doi: 10.1016/j.amc.2019.124637
[20]	P. Veeresha, D. G. Prakasha, M. A. Qurashi, D. Baleanu, A reliable technique for fractional modified Boussinesq and approximate long wave equations, Adv. Differ. Equ., 2019 (2019), 253. http://doi.org/10.1186/s13662-019-2185-2 doi: 10.1186/s13662-019-2185-2
[21]	T.-H. Zhao, O. Castillo, H. Jahanshahi, A. Yusuf, M. O. Alassafi, F. E. Alsaadi, A fuzzy-based strategy to suppress the novel coronavirus (2019-NCOV) massive outbreak, Appl. Comput. Math., 20 (2021), 160–176.
[22]	M. Nazeer, F. Hussain, M. Ijaz Khan, Asad-ur-Rehman, E. R. El-Zahar, Y.-M. Chu, Theoretical study of MHD electro-osmotically flow of third-grade fluid in micro channel, Appl. Math. Comput., 420 (2022), 126868. http://doi.org/10.1016/j.amc.2021.126868 doi: 10.1016/j.amc.2021.126868
[23]	Y.-M. Chu, B. M. Shankaralingappa, B. J. Gireesha, F. Alzahrani, M. Ijaz Khan, S. U. Khan, Combined impact of Cattaneo-Christov double diffusion and radiative heat flux on bio-convective flow of Maxwell liquid configured by a stretched nano- material surface, Appl. Math. Comput., 419 (2022), 126883. http://doi.org/10.1016/j.amc.2021.126883 doi: 10.1016/j.amc.2021.126883
[24]	T.-H. Zhao, M. Ijaz Khan, Y.-M. Chu, Artificial neural networking (ANN) analysis for heat and entropy generation in flow of non-Newtonian fluid between two rotating disks, Math. Method. Appl. Sci., in press. http://doi.org/10.1002/mma.7310
[25]	K. Karthikeyan, P. Karthikeyan, H. M. Baskonus, K. Venkatachalam, Y.-M. Chu, Almost sectorial operators on $\Psi$ -Hilfer derivative fractional impulsive integro-differential equations, Math. Method. Appl. Sci., in press. https://doi.org/10.1002/mma.7954
[26]	Y.-M. Chu, U. Nazir, M. Sohail, M. M. Selim, J.-R. Lee, Enhancement in thermal energy and solute particles using hybrid nanoparticles by engaging activation energy and chemical reaction over a parabolic surface via finite element approach, Fractal Fract., 5 (2021), 119. http://doi.org/10.3390/fractalfract5030119 doi: 10.3390/fractalfract5030119
[27]	S. Rashid, S. Sultana, Y. Karaca, A. Khalid, Y.-M. Chu, Some further extensions considering discrete proportional fractional operators, Fractals, 30 (2022), 2240026. http://doi.org/10.1142/S0218348X22400266 doi: 10.1142/S0218348X22400266
[28]	B. S. T. Alkahtani, A. Atangana, Analysis of non-homogenous heat model with new trend of derivative with fractional order, Chaos Soliton. Fract., 89 (2016), 566–571. https://doi.org/10.1016/j.chaos.2016.03.027 doi: 10.1016/j.chaos.2016.03.027
[29]	A. Atangana, D. Baleanu, New fractional derivatives with non-local and non-singular kernel theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769.
[30]	T.-H. Zhao, M.-K. Wang, G.-J. Hai, Y.-M. Chu, Landen inequalities for Gaussian hypergeometric function, RACSAM, 116, (2022), 53. https://doi.org/10.1007/s13398-021-01197-y
[31]	T.-H. Zhao, M.-K. Wang, Y.-M. Chu, On the bounds of the perimeter of an ellipse, Acta Math. Sci., 42 (2022), 491–501. https://doi.org/10.1007/s10473-022-0204-y doi: 10.1007/s10473-022-0204-y
[32]	T.-H. Zhao, Z.-Y. He, Y.-M. Chu, Sharp bounds for the weighted Hölder mean of the zero-balanced generalized complete elliptic integrals, Comput. Methods Funct. Theory, 21 (2021), 413–426. https://doi.org/10.1007/s40315-020-00352-7 doi: 10.1007/s40315-020-00352-7
[33]	L. Wang, X. Chen, Approximate analytical solutions of time fractional Whitham–Broer–Kaup equations by a residual power series method, Entropy, 17 (2015), 6519–6533. http://doi.org/10.3390/e17096519 doi: 10.3390/e17096519
[34]	S. Rashid, K. T. Kubra, H. Jafari, S. U. Lehre, A semi-analytical approach for fractional order Boussinesq equation in a gradient unconfined aquifers, Math. Method. Appl. Sci., 45 (2022), 1033–1062. http://doi.org/10.1002/mma.7833 doi: 10.1002/mma.7833
[35]	P. Veeresha, D. G. Prakasha, H. M. Baskonus, Novel simulations to the time-fractional Fisher's equation, Math. Sci., 13 (2019), 33–42. http://doi.org/10.1007/s40096-019-0276-6 doi: 10.1007/s40096-019-0276-6
[36]	M. I. El-Bahi, K. Hilal, Lie symmetry analysis, exact solutions, and conservation laws for the generalized time-fractional KdV-Like equation, J. Funct. Space., 2021 (2021), 6628130. http://doi.org/10.1155/2021/6628130 doi: 10.1155/2021/6628130
[37]	S. C. Shiralashetti, S. Kumbinarasaiah, Laguerre wavelets collocation method for the numerical solution of the Benjamina–Bona-Mohany equations, J. Taibah Univ. Sci., 13 (2019), 9–15. http://doi.org/10.1080/16583655.2018.1515324 doi: 10.1080/16583655.2018.1515324
[38]	A. Kadem, D. Baleanu, On fractional coupled Whitham–Broer–Kaup equations, Rom. J. Phys., 56 (2011), 629–635.
[39]	G. B. Whitham, Variational methods and applications to water waves, Proc. R. Soc. Lond. A, 299 (1967), 6–25. http://doi.org/10.1098/rspa.1967.0119 doi: 10.1098/rspa.1967.0119
[40]	L. J. F. Broer, Approximate equations for long water waves, Appl. Sci. Res., 31 (1975), 377–395. http://doi.org/10.1007/BF00418048 doi: 10.1007/BF00418048
[41]	D. J. Kaup, A higher-order water-wave equation and the method for solving it, Progress of Theoretical Physics, 54 (1975), 396–408. http://doi.org/10.1143/PTP.54.396 doi: 10.1143/PTP.54.396
[42]	M. Al-Qurashi, S. Rashid, F. Jarad, M. Tahir, A. M. Alsharif, New computations for the two-mode version of the fractional Zakharov-Kuznetsov model in plasma fluid by means of the Shehu decomposition method, AIMS Mathematics, 7 (2022), 2044–2060. http://doi.org/10.3934/math.2022117 doi: 10.3934/math.2022117
[43]	S. Rashid, F. Jarad, T. M. Jawa, A study of behaviour for fractional order diabetes model via the nonsingular kernel, AIMS Mathematics, 7 (2022), 5072–5092. http://doi.org/10.3934/math.2022282 doi: 10.3934/math.2022282
[44]	S. Rashid, S. Sultana, R. Ashraf, M. K. A. Kaabar, On comparative analysis for the Black-Scholes model in the generalized fractional derivatives sense via Jafari transform, J. Funct. Space., 2021 (2021), 7767848. http://doi.org/10.1155/2021/7767848 doi: 10.1155/2021/7767848
[45]	S. T. Mohyud-Din, A. Ydrm, G. Demirli, Traveling wave solutions of Whitham–Broer–Kaup equations by homotopy perturbation method, J. King Saud Univ. Sci., 22 (2010), 173–176. http://doi.org/10.1016/j.jksus.2010.04.008 doi: 10.1016/j.jksus.2010.04.008
[46]	S. Saha Ray, A novel method for travelling wave solutions of fractional Whitham-Broer-Kaup, fractional modified Boussinesq and fractional approximate long wave equations in shallow water, Math. Method. Appl. Sci., 38 (2015), 1352–1368. http://doi.org/10.1002/mma.3151 doi: 10.1002/mma.3151
[47]	B. Tian, Y. Qiu, Exact and explicit solutions of Whitham–Broer–Kaup equations in shallow water, Pure and Applied Mathematics Journal, 5 (2016), 174–180. http://doi.org/10.11648/j.pamj.20160506.11 doi: 10.11648/j.pamj.20160506.11
[48]	P. Veeresha, D. G. Prakasha, H. M. Baskonus, An efficient technique for coupled fractional Whitham–Broer–Kaup equations describing the propagation of shallow water waves, In: 4th International conference on computational mathematics and engineering sciences (CMES-2019), Cham: Springer, 2019, 49–75. http://doi.org/10.1007/978-3-030-39112-6_4
[49]	J. Singh, D. Kumar, R. Swroop, Numerical solution of time- and space-fractional coupled Burgers' equations via homotopy algorithm, Alex. Eng. J., 55 (2016), 1753–1763. http://doi.org/10.1016/j.aej.2016.03.028 doi: 10.1016/j.aej.2016.03.028
[50]	S. J. Liao, Homotopy analysis method: A new analytical technique for nonlinear problems, Commun. Nonlinear Sci. Numer. Simul., 2 (1997), 95–100. https://doi.org/10.1016/S1007-5704(97)90047-2 doi: 10.1016/S1007-5704(97)90047-2
[51]	S. J. Liao, Homotopy analysis method: a new analytic method for nonlinear problems, Appl. Math. Mech., 19 (1998), 957–962. http://doi.org/10.1007/BF02457955 doi: 10.1007/BF02457955
[52]	K. S. Aboodh, The new integral transform "Aboodh transform", Global Journal of Pure and Applied Mathematics, 9 (2013), 35–43.
[53]	H. Bulut, D. Kumar, J. Singh, R. Swroop, H. M. Baskonus, Analytic study for a fractional model of HIV infection of CD4+T lymphocyte cells, Math. Nat. Sci., 2 (2018), 33–43. http://doi.org/10.22436/mns.02.01.04 doi: 10.22436/mns.02.01.04
[54]	P. Veeresha, D. G. Prakasha, Solution for fractional Zakharov-Kuznetsov equations by using two reliable techniques, Chinese J. Phys., 60 (2019), 313–330. http://doi.org/10.1016/j.cjph.2019.05.009 doi: 10.1016/j.cjph.2019.05.009
[55]	D. G. Prakasha, P. Veeresha, Analysis of lakes pollution model with Mittag-Leffler kernel, J. Ocean Eng. Sci., 5 (2020), 310–322. http://doi.org/10.1016/j.joes.2020.01.004 doi: 10.1016/j.joes.2020.01.004
[56]	M. H. Cherif, D. Ziane, A new numerical technique for solving systems of nonlinear fractional partial differential equations, Int. J. Anal. Appl., 15 (2017), 188–197.
[57]	M. A. Awuya, D. Subasi, Aboodh transform iterative method for solving fractional partial differential equation with Mittag–Leffler kernel, Symmetry, 13 (2021), 2055. http://doi.org/10.3390/sym13112055 doi: 10.3390/sym13112055
[58]	G. M. Mittag-Leffler, Sur La nonvelle Fonction $E_{\alpha}(x)$ , C. R. Acad. Sci. Paris, (Ser. II), 137 (1903), 554–558.
[59]	S. M. El-Sayed, D. Kaya, Exact and numerical travelling wave solutions of Whitham–Broer–Kaup equations, Appl. Math. Comput., 167 (2005), 1339–1349. http://doi.org/10.1016/j.amc.2004.08.012 doi: 10.1016/j.amc.2004.08.012
[60]	M. Rafei, H. Daniali, Application of the variational iteration method to the Whitham-Broer-Kaup equations, Comput. Math. Appl, 54 (2007), 1079–1085. http://doi.org/10.1016/j.camwa.2006.12.054 doi: 10.1016/j.camwa.2006.12.054

This article has been cited by:

1.	Hira Soomro, Nooraini Zainuddin, Hanita Daud, Joshua Sunday, Noraini Jamaludin, Abdullah Abdullah, Mulono Apriyanto, Evizal Abdul Kadir, Variable Step Block Hybrid Method for Stiff Chemical Kinetics Problems, 2022, 12, 2076-3417, 4484, 10.3390/app12094484
2.	Zeeshan Ali, Faranak Rabiei, Kamyar Hosseini, A fractal–fractional-order modified Predator–Prey mathematical model with immigrations, 2023, 207, 03784754, 466, 10.1016/j.matcom.2023.01.006
3.	Sümeyra Uçar, Analysis of hepatitis B disease with fractal–fractional Caputo derivative using real data from Turkey, 2023, 419, 03770427, 114692, 10.1016/j.cam.2022.114692
4.	Mohammad Partohaghighi, Ali Akgül, Rubayyi T. Alqahtani, New Type Modelling of the Circumscribed Self-Excited Spherical Attractor, 2022, 10, 2227-7390, 732, 10.3390/math10050732
5.	G.M. Vijayalakshmi, Roselyn Besi. P, A fractal fractional order vaccination model of COVID-19 pandemic using Adam’s moulton analysis, 2022, 8, 26667207, 100144, 10.1016/j.rico.2022.100144
6.	SHAIMAA A. M. ABDELMOHSEN, SHABIR AHMAD, MANSOUR F. YASSEN, SAEED AHMED ASIRI, ABDELBACKI M. M. ASHRAF, SAYED SAIFULLAH, FAHD JARAD, NUMERICAL ANALYSIS FOR HIDDEN CHAOTIC BEHAVIOR OF A COUPLED MEMRISTIVE DYNAMICAL SYSTEM VIA FRACTAL–FRACTIONAL OPERATOR BASED ON NEWTON POLYNOMIAL INTERPOLATION, 2023, 31, 0218-348X, 10.1142/S0218348X2340087X
7.	D.A. Tverdyi, R.I. Parovik, A.R. Hayotov, A.K. Boltaev, Распараллеливание численного алгоритма решения задачи Коши для нелинейного дифференциального уравнения дробного переменного порядка с помощью технологии OpenMP, 2023, 20796641, 87, 10.26117/2079-6641-2023-43-2-87-110
8.	Khalid Hattaf, A New Class of Generalized Fractal and Fractal-Fractional Derivatives with Non-Singular Kernels, 2023, 7, 2504-3110, 395, 10.3390/fractalfract7050395
9.	Ali Raza, Ovidiu V. Stadoleanu, Ahmed M. Abed, Ali Hasan Ali, Mohammed Sallah, Heat transfer model analysis of fractional Jeffery-type hybrid nanofluid dripping through a poured microchannel, 2024, 22, 26662027, 100656, 10.1016/j.ijft.2024.100656
10.	Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen, New solutions of time-space fractional coupled Schrödinger systems, 2023, 8, 2473-6988, 27033, 10.3934/math.20231383
11.	Muhammad Shahzad, Soma Mustafa, Sarbaz H A Khoshnaw, Inference of complex reaction mechanisms applying model reduction techniques, 2024, 99, 0031-8949, 045242, 10.1088/1402-4896/ad3291
12.	B. El Ansari, E. H. El Kinani, A. Ouhadan, Symmetry analysis of the time fractional potential-KdV equation, 2025, 44, 2238-3603, 10.1007/s40314-024-02991-1
13.	Muhammad Farman, Changjin Xu, Perwasha Abbas, Aceng Sambas, Faisal Sultan, Kottakkaran Sooppy Nisar, Stability and chemical modeling of quantifying disparities in atmospheric analysis with sustainable fractal fractional approach, 2025, 142, 10075704, 108525, 10.1016/j.cnsns.2024.108525
14.	Hira Khan, Gauhar Rahman, Muhammad Samraiz, Kamal Shah, Thabet Abdeljawad, On Generalized Fractal-Fractional Derivative and Integral Operators Associated with Generalized Mittag-Leffler Function, 2025, 24058440, e42144, 10.1016/j.heliyon.2025.e42144
15.	Emmanuel Kengne, Ahmed Lakhssassi, Dynamics of stochastic nonlinear waves in fractional complex media, 2025, 542, 03759601, 130423, 10.1016/j.physleta.2025.130423

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)