An Accurate Approach to Simulate the Fractional Delay Differential Equations

: The fractional Legendre polynomials (FLPs) that we present as an effective method for solving fractional delay differential equations (FDDEs) are used in this work. The Liouville–Caputo sense is used to characterize fractional derivatives. This method uses the spectral collocation technique based on FLPs. The proposed method converts FDDEs into a set of algebraic equations. We lay out a study of the convergence analysis and ﬁgure out the upper bound on error for the approximate solution. Examples are provided to demonstrate the precision of the suggested approach.


Introduction
In the hope of simulating a variety of physical events, it is crucial to use differential equations based on fractional-order operators [1,2].Some applications of these equations include earthquake simulation, damping principles, fluid dynamics, biology, physics, and engineering [3][4][5][6][7].Numerous techniques have been employed to locate appropriate solutions to these equations depending on their physical structure.The majority of the current systems are flawed and unworkable considering the actual physical existence of these problems.The numerical findings show the full reliability of the proposed algorithms [8][9][10].Since explicit exact solutions for FDEs are still lacking, approximation and numerical techniques such as the spectral collocation method (SCM), FEM, FVM, and many others have been used to solve many fractional models [11][12][13][14][15].The paper's primary goal is to use the mentioned algorithm ( [16,17]) to obtain the numerical solution of the following FDDE using SCM depending on the fractional Legendre polynomials [18]: where where the order of the Liouville-Caputo fractional is α.The delay function is a function g that satisfies condition a ≤ g(x) ≤ x, x ∈ [a, b] and g(x) ∈ C[a, b].Also, we shall calculate an upper bound for the estimated solution's resulting error and the convergence will be studied.In ( [19,20]), the authors used spline functions to study both the error and stability analysis for first-order delay differential equations.We plan to discretize Equation (1) using the Legendre collocation method to produce an algebraic system of equations (linear/nonlinear) and solve it.
The manuscript is structured as follows: Section 2 presents the definitions and approximate formula for the fractional derivative.Section 3 provides some thoughts on the recently revealed fractional-order Legendre polynomials.The polynomials approximation and error estimation are discussed in Section 4. The implementation of the Legendre spectral method for solving FDDE is introduced in Section 5. Section 6 also discusses numerical simulation.We provide the conclusions and some observations in Section 7.

Preliminaries and Notations
Definition 1.The following formula defines the fractional derivative of order α of Liouville-Caputo ( [15,21]): where Additionally, a generalization of Taylor's formula is introduced based on the fractional derivative of Liouville-Caputo [22].Assume that
We shall introduce the fractional-order Legendre polynomials in this section.By changing the variables t = x ν and ν > 0 on L * k (t), these polynomials will be obtained.We represent FLPs, L * k (x ν ) by P ν k (x).We note that the recurrence formula for P ν k (x) may be deduced from the recurrence relation of L * k (t) [18]: The analytic form of the FLPs, P ν k (x) of degree νk is given by: Note that P ν k (0) = (−1) k and P ν k (1) = 1.
Proof.The formula can be obtained directly by taking t = x ν , in the orthogonality condition

Approximating Polynomials and Estimating Errors
Let a non-negative, integrable, real-valued function spanning the range (0, 1) be denoted by w ν (x) = x ν−1 , ν > 0. In the first (m + 1)-terms of FLPs, the function u(x) ∈ L 2 w ν , which is a square-integrable function in [0, 1], can be represented and approximated as follows: The error estimate for this approximation in w ν -norm is provided by the following theorem.
Theorem 1.The error bound can be calculated using the following formula if u ν m is the best Proof.As an approximate representation of u(x) and denoted by y(x), consider the generalized Taylor's polynomial as follows: where it is known that the following error bound exists and: this together with the well-known Stirling formula (n! = √ 2π n n+1/2 e −n+θ/12n , for some (0 < θ < 1) lead to the desired result.
The primary formula for calculating an approximation of the Liouville-Caputo fractional derivative is given in the following theorem.Theorem 2. Let u(x) be as described in (6) and α > 0, then: Proof.We can obtain this formula directly by using the linearity properties of the Liouville-Caputo operation and then employing Equation (3) in Formula ( 5) with some simplification. where , j = 0, 1, . . . .

Proof.
The proof can be found in [24].
From ( 9) in (2), and use the properties of the Formula (5), we can obtain the following equations: Using Equation (11) and α equations of B.Cs gives (m + 1) algebraic equations in the unknown c i , i = 0, 1, . . ., m, which can be solved.Consequently, u(x) for (1) can be found.

Numerical Results
Example 1.Consider the following differential equation for a linear fractional delay: with the initial condition: Using the provided method and m = 2, we arrive at the following approximation of the solution: Using Equation (11) with α = 0.5, we have: with p = 0, 1 where x p are roots of P ν 2 (x) and x 0 = 0.0446582, x 1 = 0.622008, at ν = 0.5.Using Equations ( 14) and ( 15) we get: Now, solving Equations ( 16) and ( 17), we get: Thus, using Equation ( 15) we can get: u(x) Example 2. Consider the following differential equation for a linear fractional delay: The exact solution of Equation ( 18) Using the provided method and m = 4, we arrive at the following approximation of the solution: Using Equation (11), we have: with p = 0, 1, 2, 3 where x p are roots of P ν 4 (x) and x 0 = 0.004821, x 1 = 0.108961, x 2 = 0.448887, x 3 = 0.865957, at ν = 0.5.
The approximate solution and the exact one are given in Figure 1 at α = 1 and ν = 0.5.But, the approximate solution with α(= 0.75, 0.85) is given in Figure 2.These graphs demonstrate that the exact solution and numerical solution are in very good agreement, and the approximate solution exhibits the same behavior.Example 3. Consider the following non-linear fractional delay differential equation: where h(x) = Γ(2ν+1) Γ(2ν−0.5)x 2ν−1.5 − (x − 0.5) 2ν − x 6ν , with Using the provided method and m = 3, we arrive at the following approximation of the solution: Using Equation ( 11), we have: with p = 0, 1, 2 where x p are roots of P ν 3 (x).Using Equations ( 24) and ( 25), we obtain: Solving ( 26)-( 28), we obtain: Thus, using Equation ( 25) at ν = 1.25, we can obtain: Example 4 ((Ikeda system-One delay)).Consider the Ikeda delay system with one delay, which is described in ( [26,27]): where κ, λ, µ, and τ are constants; with Using the provided method and m = 10, we arrive at the approximation of the solution as in Equation (25).Using Equation (11), we have with p = 0(1)8 where t p are roots of P ν 9 (t).Using Equations ( 25) and (30), we can obtain: Now, solving the system of Equations ( 31)-(32), we obtain the approximate solution of this model (29).
To conclude this numerical analysis with the simulation; we define the residual error function (REF) as: The absolute relative error drops to zero in all situations when the residual is minimal (REF(t, m) → 0), indicating that the solution tends to the exact one.Since the precise answer is unknown in the situation of α in fractional order, we shall employ this form of error.Finally, REF includes other forms; for further information, see [28].
Here, we define the REF to conduct a full numerical study with the simulation: Using the provided method, we present a numerical simulation of this model (34) using the suggested approach in Figures 6-8

Conclusions and Remarks
This work's primary focus is on presenting fractional Legendre polynomials with states and demonstrating some of their characteristics.In addition, the Legendre SCM is used in this article to solve fractional delay differential equations.The suggested problem is reduced to a system of algebraic equations, which can be solved using an appropriate numerical method, leveraging the properties of the new fractional Legendre polynomials.We examine the convergence analysis, the generated formula's accuracy, and the approximation of the solution.The numerical simulations show that this approach to applying the Liouville-Caputo derivative is an effective way to solve FDDE.In addition, Only a few shifted fractional Legendre polynomials are required to produce a good outcome.By including new terms in the series, overall error rates can be reduced.

Figure 1 .
Figure 1.The numerical solution and the exact solution of Example 2.

Figure 2 .
Figure 2. The approximate solution with different values of α of Example 2.

Figure 5 .
Figure 5.The numerical solution of Example 4, with m = 15; with different values of τ (a) and α (b).