Study of a Nonlinear System of Fractional Differential Equations with Deviated Arguments Via Adomian Decomposition Method.

This paper studies a system of nonlinear fractional differential equations (FDEs) with deviated arguments. Many linear and nonlinear problems are faced in the real-life. Generally, linear problems are solved quickly, but some difficulties appear while solving nonlinear problems. Our purpose is to approximate those solutions numerically via the Adomian decomposition method (ADM). Here, our main goal is to apply the ADM to solve higher-order nonlinear system of FDEs with deviated arguments. We prove the existence and uniqueness of the solution using Banach contraction principle. Moreover, we plot the figures of ADM solutions using MATLAB.

Introduction

“Science is a differential equation,” said Alan Turing, and Paul Ormerod “Baseball players or cricketers do not need to be able to solve explicitly the nonlinear differential equations which govern the flight of the ball. They just catch it.” Nonlinear differential equations describe many real-world physical phenomena. To understand the nature of these phenomena, we must first solve differential equations. “In order to solve this differential equation, you look at it until a solution occurs to you,” George Polya explained. In 1980, George Adomian proposed a new iterative scheme known as ADM. This method provides analytical solutions to both linear and nonlinear differential equations.

Classical differential equations cannot adequately describe more and more phenomena as science and technology advance. Various physical processes, for example, have memory and heritability properties that the classical local differential operators cannot adequately represent. Many great mathematicians such as Euler, Liouville, Riemann, Caputo, and Letnikov developed a new excellent tool to describe these nonlocal processes (fractional differential equations described by nonlocal operators) [21, 23]. Fractional Calculus was born in 1695. G.F.A. de L Hospital thought about what happens if the order is 1/2. And in 1697, G.W. Leibniz used fractional derivatives of order 1/2 and made some remarks on it. However, in 1819, S.F. Lacroix mentioned the derivative of arbitrary order in his text on differential and integral calculus. N.H. Abel gave the first application of fractional calculus in 1823.

Fractional calculus has been included in recent studies on CoViD-19 [3, 17, 20]. Many studies have been conducted in the field of mathematics, and it has been shown that differential equations using fractional operators are effective in demonstrating epidemic models linked to many infectious illnesses [7, 12]. Two leading implementations of fractional calculus are in epidemiological and biomathematical models [6, 18, 19]. Ahmad et al. [2] performed simulations of a fractional model for CoViD-19 transmission, taking into account various values of the non-integer order derivative and came to the conclusion that the value of = 0.97 best matched the actual data. Furthermore, Zhang et al. [29] created a non-integer order model for the dynamics of CoViD-19. The authors investigated the stability of the system and reproduction number.

FDEs with deviated arguments have many applications in science and engineering, including fractals theory, chemistry, biology, physics, neural network, weather prediction model, etc. [19, 26]. Brauer et al. [5] presented logistic equations, which are particularly applicable to epidemic systems. Xu et al. [26] explored the effect of numerous time delays on fractional-order neural network bifurcation. There are various methods to solve fractional differential equations: the homotopy perturbation method [9, 13, 14, 22], the Adomian decomposition method [1, 8, 15, 16, 25, 27], the polynomial least square method [4], and so on [24].

In [10], Duan et al. provided a review of ADM and its application to FDEs. Evans and Raslan [11] applied the Adomian method to solve a particular ordinary delay differential equations in which the delay is located in the linear or nonlinear part, where the history function is not necessary. In [16], Li and Pang provided an application of ADM to a nonlinear system. In [25], Saeed and Rahman studied the ADM for solving the system of delay differential equation. In [27], Ziada studied the nonlinear system of fractional differential equations via ADM, and the fractional order rabies model was solved as an application. However, in [28] Ziada studied the analytical and numerical solutions of a multi-term nonlinear differential equation with deviated arguments.

Motivated by the works of Ziada [27, 28] as well as Saeed and Rahman [25], we construct the system (1)–(2). We extend the work of [25] for higher order nonlinear FDEs with deviated arguments.

Delay differential equations are far more complicated than traditional ordinary differential equations, they explain many processes found in several fields such as biology, medicine, chemistry, economics, engineering and physics. Systems of FDEs have many applications in engineering and science, including electrical networks, control theory, fractals theory, viscoelasticity, optical and neural network systems. This paper aims to discuss the approximate solution of a nonlinear system of FDEs with deviated arguments via ADM (an algorithm that uses a decomposition technique). Here, our main goal is to apply the ADM to solve higher-order nonlinear systems of FDEs with deviated arguments. This method has numerous advantages. It is very simple to use and can solve a wide range of nonlinear systems, such as ordinary and partial differential equations, fractional delay differential equations, and so on. It avoids the Picard method’s time-consuming integrations. It decomposes the solution into a series with easily computed components. It has the advantage of converging to the exact solution.

The remaining paper is designed as “Formulation of the Problem” section introduces the formulation of the problem and ADM’s iterative scheme. “Basic Definitions” section contains some useful basic definitions. In the next two sections, we prove the existence, uniqueness and convergence of the solution of the system (1)–(2). Numerical examples have been provided in “Numerical Examples” section. The conclusion is added in the last section.

Formulation of the Problem

Consider the following higher-order nonlinear system of FDEs with deviated argumentswith initial conditionswhere denote the Caputo derivative of order   Here, we use Caputo fractional derivative amongst a variety of definitions for fractional order derivatives as it is suitable for describing various phenomena, since the initial values of the function and its integer order derivatives have to be specified. are nonlinear operators that satisfy Lipschitz condition with Lipschitz constant , such as are unknown functions, are given continuous functions and are given constants. In order to solve the problem (1) with (2) by using the ADM, performing the fractional integral to both sides of (2), we haveAdomian’s method defines the solution by seriesSo that, the components will be determined recursively. Moreover, the method defines the nonlinear term by the Adomian polynomialswhere are Adomian polynomials that can be generated for all forms of nonlineraity aswhere is a parameter.

In view of (5) and (6), (4) becomesTo determine the components First we identify the zero component by the terms and where represent the non-homogeneous parts of Thus, the recurrence relation isWe can approximate the solution by the truncated series

Basic Definitions

Definition 3.1

[21, 23] (a) Caputo fractional derivativewhere(b)(c)

Existence and Uniqueness of Solution

Define the operator where is the Banach space the space of all continuous functions on equipped with the norm

Theorem 4.1

If satisfy the Lipschitz condition (3) and where then the system (1)–(2) has a unique solution on

Proof

Define the operator asLet thenSince therefore the mapping is contraction. By Banach contraction principle, there exists a unique solution This completes the proof.

Proof of Convergence

Theorem 5.1

The series solution (5) of the system (1)–(2) using ADM converges if and where

Define the sequence of partial sum as

Since,So, we haveFurther, we prove that is a Cauchy sequence in Let be two arbitrary partial sums such that thenLet then

Using triangle inequality, we haveSince and then Hence,Since therefore as Hence, is a Cauchy sequence in and thus the series (5) converges. The proof is completed.

Numerical Examples

Example 1

Consider the following nonlinear systemsubject to the initial conditionswhere

On applying ADM to system (12), we obtain the following schemewhere and represent the Adomian polynomials of nonlinear terms and respectively.

Using the relations (13)–(15), the first four terms of the series solutions areFigure 1A–C show ADM solution of and at different values of respectively.

Fig. 1

A–C, ADM Sol. of

Example 2

Consider the following systemsubject to the initial conditions

On applying ADM to system (19), we have the following recursive relationswhere and represent the Adomian polynomials of nonlinear terms and respectively.

Using the relations (20)–(22), the first few terms of the series solution areFigure 2 shows the ADM solution of

Fig. 2

ADM Sol. of

Example 3

Consider the following systemsubject to the initial conditionswhere

Applying ADM to system (26), we get the following schemewhere and represent the Adomian polynomials of nonlinear terms and respectively.

Using the relations (27)–(28), the first few terms of the series solution areFigure 3A, B show the ADM solution of

Fig. 3

A, B, ADM Sol. of

Example 4

Consider the following nonlinear systemsubject to the initial conditions Applying ADM to system (31) leads to the following schemewhere and represent the Adomian polynomials of nonlinear terms and respectively.

Using the relations (32)–(33), the first few terms of the series solution areFigure 4 shows the ADM solution of

Fig. 4

ADM Sol. of

In general, finding the exact solution of every differential equation is difficult, particularly the higher-order non-linear fractional differential equations with deviated arguments. The question is how one can check the accuracy of the method in the absence of the exact solution. Therefore, we give the following example with the known exact solution and show the accuracy of the proposed method.

Example 5

Consider the following nonlinear systemsubject to the initial conditionswhere , which has the exact solution and at (Fig. 5).

Fig. 5

A, B, ADM and exact sol. of

On applying ADM to system (36), we obtain the following schemewhere and represent the Adomian polynomial of nonlinear terms and respectively. Using the relations (37)–(38), the first few terms of the series solution areTables  1 and  2 show the exact and numerical values of and . Figure  5 A, B show the solutions of and , respectively.

Table 1

The exact and numerical values of the solution of example 5

t	Exact solution	ADM solution	Error
0	1.0000	1.0000	0
0.1	0.9950	0.9950	0.0000
0.2	0.9801	0.9801	0.0000
0.3	0.9553	0.9553	0.0000
0.4	0.9211	0.9211	0.0000
0.5	0.8776	0.8776	0.0000
0.6	0.8253	0.8254	0.0001
0.7	0.7648	0.7650	0.0002
0.8	0.6967	0.6971	0.0004
0.9	0.6216	0.6223	0.0007
1	0.5403	0.5417	0.0014

Table 2

The exact and numerical values of the solution of example 5

t	Exact solution	ADM solution	Error
0	0	0	0
0.1	0.0998	0.0998	0.0000
0.2	0.1987	0.1987	0.0000
0.3	0.2955	0.2955	0.0000
0.4	0.3894	0.3893	0.0001
0.5	0.4794	0.4792	0.0003
0.6	0.5646	0.5640	0.0006
0.7	0.6442	0.6428	0.0014
0.8	0.7174	0.7147	0.0027
0.9	0.7833	0.7785	0.0048
1	0.8415	0.8333	0.0081

Conclusion

The focus of this paper is to approximate the solution of a nonlinear system of FDEs with deviated arguments using a simple method. Using Banach contraction principle, we prove the existence and uniqueness of the solution. As we know that in real-life, many linear and nonlinear problems occur in the form of a differential equation. Some difficulties occur while solving nonlinear FDEs. Therefore, we apply the ADM method to solve the higher-order nonlinear system of FDEs with deviated arguments and provide some numerical examples to show the effectuality of the method. We plot the figures of ADM solutions using MATLAB.