The Novel Mittag-Leffler–Galerkin Method: Application to a Riccati Differential Equation of Fractional Order

Abstract

We present a new numerical approach to solving the fractional differential Riccati equations numerically. The approach—called the Mittag-Leffler–Galerkin method—comprises the finite Mittag-Leffler function and the Galerkin method. The error analysis of the method was studied. As a result, we present two theorems by which the error can be bounded. In addition to error analysis, the residual correction method, which allows us to estimate the error and obtain new approximate solutions, is also presented. To show how the method is applied, and the efficiency of the proposed method, some test examples were considered. When the numerical results obtained were examined, it was found that while the method achieves better results than some of the known methods in the literature, it also achieves results that are similar to those of others of the known methods.

1. Introduction

Fractional differential equations (FDEs) have been used to describe real-life phenomena such as continuum mechanics [1], viscoelasticity [2], finance [3], optimal control [4], variational problems [5], hydrologic modeling [6] and fluid mechanics [7], amongst others [8,9]. Due to the difficulty of obtaining exact solutions, the importance of developing effective methods for the numerical solutions of FDEs has been recognized in recent decades. The main methods used to solve FDEs include radial basis functions [10], fractional finite volume [11], Adomian decomposition [12], operational methods [13] and other numerical approaches [14,15,16]. In this study, we propose a new numerical solution method of solving the FDRE, defined as

$$D^{\nu }x\left(t\right)+a\left(t\right)x^2\left(t\right)+b\left(t\right)x\left(t\right)=g\left(t\right),0<\nu \le 1,0\le t\le 1,$$

(1)

and the initial condition

$$x\left(0\right)=x_0.$$

(2)

Here, $x\left(t\right)$ is the unknown function, $a\left(t\right),b\left(t\right)$ and $g\left(t\right)$ are known functions defined in $[0,1]$ and continuous, and $x_0$ is a real constant.

Recently, spectral approaches have been applied to solving different types of FDEs. Esmaeili and Shamsi [17] considered a family of FDEs, and solved it by a pseudo-spectral method. Zhang et al. [18] solved the one-dimensional nonlinear space fractional Schrödinger equation, using the Crank–Nicolson–Galerkin–Legendre spectral method. Mokhtary and Ghoreishi [19] used the tau spectral method for the solutions of nonlinear fractional integrodifferential equations. Brawy et al. [20] introduced an operational approximation method, based on the spectral collocation method for the solution of fractional Schrödinger equations. Vanani and Aminataei [21] improved the algebraic formulation of fractional partial differential equations, by using matrix–vector multiplication representation, and then applied an operational approach of the tau method. Doha et al. [22] proposed a spectral method for the solution of the fractional subdiffusion equation. The approach was based on the shifted Legendre tau spectral method. Fan et al. [23] proposed a Galerkin finite element method for solving the fractional wave equation: they discretized the problem of the Crank–Nicholson scheme, and presented the stability and convergence of the numerical scheme. Saadatmandi and Dehghan [24] used a Jacobi–Gauss–Lobatto and Gauss–Radau collocation method, based on shifted Jacobi polynomials, to solve fractional Fokker–Planck equations. Kazem [25] employed an integral operational matrix, based on Jacobi polynomials, to solve FDEs.

The differential Riccati equation (DRE) is used to describe miscellaneous engineering and physical phenomena, such as the flow of rivers, the transmission line phenomenon, stochastic control, dynamic games and financial mathematics [26,27,28]. The FDRE, which is a generalization of the DRE, has many applications in science and engineering [29,30,31], so various solution strategies have been suggested. Ozturk et al. [32] used the Taylor collocation method, converting the FDRE into a system of nonlinear algebraic equations, and then solving the system. Balaji [33] applied the Legendre wavelet operational matrix method to FDRE, to obtain its approximate solution. Mokhtary and Ghoreishi [34] introduced an operational method constituted of shifted Jacobi polynomials, to solve FDREs. Kashkari and Syam [35] used the Legendre operational matrix of fractional integration, to derive a numerical solution for FDREs. According to Jafari et al., [36] adopted a modified variation iteration method for FDREs, taking into account Adomian polynomials for nonlinear terms. Bota and Caruntu [37] applied the polynomial least squares method, to find an analytical solution for FDREs. Merdan [38] applied the fractional variational iteration method, to obtain an approximate analytical solution for nonlinear FDREs. Odibat and Momani [39] applied a modification of He’s homotopy perturbation method to the quadratic FDRE. The homotopy analysis transform method, based on a combination of the homotopy analysis method and the Laplace decomposition method, was employed by Saad and Al-Shomrani [40] to solve FDREs. Haq et al. [41] applied the variation of parameters method, to obtain the analytical solutions of nonlinear quadratic FDREs. Sakar et al. [42] applied an iterative reproducing kernel Hilbert space method, to obtain the solutions of FDREs. Yuzbasi [43] studied with the Bernstein collocation method, to obtain the numerical solutions of FDREs. Li et al. [44] derived the Haar wavelet operational matrix method, to solve FDREs: they simplified the calculation of the nonlinear term using the block pulse function. Raja et al. [45] introduced a new computational intelligence technique, based on artificial neural networks and sequential quadratic programming, to solve nonlinear quadratic FDREs.

In this paper, we introduce a new method of solving FDREs. We approximate the solution by an expansion in the finite Mittag-Leffler function. By applying the Galerkin method, the FDRE is reduced to a nonlinear system of algebraic equations. Solving these equations gives the approximate solution to the problem. The rest of the paper is organized as follows: in Section 2, some necessary definitions of fractional calculus and the finite Mittag-Leffler function are presented, along with some properties of the function; the proposed method is introduced in Section 4; our analysis of the error is presented in Section 5; some theorems for the error analysis of the method, with the residual correction procedure, are presented; in Section 6, we give some test examples, to illustrate the application steps of the method; we compare the results of the proposed method to the results of some other methods; finally, in Section 7, we summarize the results.

2. Fractional Calculus and Finite Mittag–Leffler Function

In this section, we first give some fundamental definitions of fractional calculus. Then, the properties of the finite Mittag-Leffler function (MLF) and its fractional derivative are introduced.

2.1. Fractional Calculus and Mittag-Leffler Function

Definition 1

([46]). The fractional integral of order $\nu >0$ with the lower limit zero for a function x is defined as

$$\begin{array}{c}\hfill I^{\nu }x\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(\nu \right)}{\int }_0^t\frac{x\left(s\right)}{{(t-s)}^{1-\nu }}ds,\nu >0.\end{array}$$

(3)

Here, the right-hand side is defined pointwise on $[0,\infty )$, and $\mathsf{\Gamma }(.)$ is the Gamma function.

Definition 2

([46]). The Caputo derivative of order ν with the lower limit zero for a function x is defined as

$$\begin{array}{ccc}\hfill D^{\nu }x\left(t\right)& =& \frac{1}{\mathsf{\Gamma }(\lceil \nu \rceil -\nu )}{\int }_0^t\frac{x^{(\lceil \nu \rceil )}\left(s\right)}{{(t-s)}^{\nu +1-\lceil \nu \rceil }}ds\hfill \end{array}$$

(4)

$$\begin{array}{ccc}& =& I^{(\lceil \nu \rceil -\nu )}{x}^{(\lceil \nu \rceil )}\left(t\right),t>0,\nu >0,\hfill \end{array}$$

(5)

where $\lceil \nu \rceil$ is the ceiling function of ν.

The following properties [47] apply to the Caputo fractional derivative operator: we have, for constants ${\xi }_i,i=1,2,...,N$,

$$\begin{array}{c}\hfill D^{\nu }\sum _{i=1}^N{\xi }_i{x}_i\left(t\right)=\sum _{i=1}^N{\xi }_i{D}^{\nu }{x}_i\left(t\right),\end{array}$$

(6)

as well as, from [48],

$$\begin{array}{c}\hfill D^{\nu }{t}^N=\frac{\mathsf{\Gamma }(N+1)}{\mathsf{\Gamma }(N+1-\nu )}{t}^{N-\nu },N>\nu -1.\end{array}$$

(7)

In the case of $\nu$ as an integer, the Caputo differential operator will coincide with the usual differential operator.

2.2. Mittag-Leffler Function

The Mittag-Leffler function $E^{\xi ,\eta }$ is a function that is dependent on two parameters, $\xi$ and $\eta$. When the real component of $\xi$ is strictly positive, it can be described by the following series:

Definition 3

([46]). The MLF of two-parameter $\xi ,\eta$ is defined by

$$\begin{array}{c}\hfill E^{\xi ,\eta }\left(t\right)=\sum _{i=0}^{\infty }\frac{t^i}{\mathsf{\Gamma }(\xi i+\eta )},\xi >0,\eta >0,t\in \mathbb{R},\end{array}$$

where $\mathsf{\Gamma }(\xi i+\eta )$ is the gamma function. As a special case, we have $E^{1,1}\left(t\right)=e^t$.

3. Finite Mittag-Leffler Function and Its Fractional Derivative

The novel definition of the two-parameter finite Mittag-Leffler function of any integer i is:

Definition 4.

The finite MLF of two-parameter $\xi ,\eta$ can be defined as

$$\begin{array}{c}\hfill E_i^{\xi ,\eta }\left(t\right)=\sum _{k=0}^i\frac{{(-1)}^k{t}^k}{\mathsf{\Gamma }(\xi k+\eta )},\xi >0,\eta >0,t\in \mathbb{R},\end{array}$$

(8)

that is:

$$\begin{array}{c}\hfill E_i^{\xi ,\eta }\left(t\right)=\frac{{(-1)}^i{t}^i}{\mathsf{\Gamma }(\xi i+\eta )}+\frac{{(-1)}^{i-1}{t}^{i-1}}{\mathsf{\Gamma }\left(\xi \right(i-1)+\eta )}+\cdots +\frac{-t}{\mathsf{\Gamma }(\xi +\eta )}+\frac{1}{\mathsf{\Gamma }\left(\eta \right)}.\end{array}$$

Based on the above definition, we can write

$$E_i^{\xi ,\eta }\left(t\right)=\frac{{(-1)}^i{t}^i}{\mathsf{\Gamma }(\xi i+\eta )}+\mathrm{polynomials}\mathrm{of}\mathrm{degree}<\mathrm{i}.$$

(9)

Using (7), the fractional-order derivative of the Mittag-Leffler function (8) can be calculated as

$$\begin{array}{ccc}\hfill D^{\nu }{E}_i^{\xi ,\eta }\left(t\right)& =& \sum _{k=0}^i\frac{D^{\nu }{(-1)}^k{t}^k}{\mathsf{\Gamma }(\xi k+\eta )}\hfill \\ & =& \sum _{k=0}^i\frac{\mathsf{\Gamma }(k+1){(-1)}^k}{\mathsf{\Gamma }(k+1-\nu )}\frac{t^{k-\nu }}{\mathsf{\Gamma }(\xi k+\eta )},\xi >0,\eta >0,t\in \mathbb{R}.\hfill \end{array}$$

4. The Fundamental Concepts of the Mittag-Leffler–Galerkin Method

In this section, we apply the Galerkin method, Equation (13), which has been used to solve problems in structural mechanics, dynamics, fluid flow, hydrodynamic stability, magnetohydrodynamics, heat and mass transfer, acoustics, microwave theory, neutron transport, etc. Problems governed by ordinary differential equations, partial differential equations and integral equations have been investigated via Galerkin formulations. Steady, unsteady and eigenvalue problems have proved to be equally amenable to the Galerkin treatment. Essentially, any problem for which governing equations can be written down is a candidate for a Galerkin method [49].

For the first time, we introduce a novel numerical method, namely the Mittag-Leffler–Galerkin (MLG) method, for solving FDREs.

Let $x\left(t\right)$ be the exact solution of Equations (1) and (2). We will use the proposed new MLG method to approximate the exact solution $x\left(t\right)$, as follows:

$$x\left(t\right)\simeq x_{{}_N}\left(t\right)=\sum _{i=0}^N{c}_i{E}_i^{\xi ,\eta }\left(t\right)=\mathrm{C}\mathbb{E}\left(t\right).$$

(10)

Here, $\mathrm{C}$ is an unknown constant matrix of size $1\times (N+1)$ that must be determined, and $\mathbb{E}\left(t\right)$ is a matrix of size $(N+1)\times 1$ that consists of Mittag-Leffler-basis polynomial elements, defined as

$$\begin{array}{c}\hfill \mathrm{C}={\left[{\mathrm{c}}_0,\dots ,{\mathrm{c}}_N\right]}_{1\times (N+1)},\mathbb{E}\left(t\right)={\left[E_0^{\xi ,\eta }\left(t\right),E_1^{\xi ,\eta }\left(t\right),\dots ,E_N^{\xi ,\eta }\left(t\right)\right]}_{(N+1)\times 1}.\end{array}$$

Directly from the Mittag-Leffler-basis polynomial elements, the fractional derivative $D^{\nu }x\left(t\right)$ can be stated as

$$\begin{array}{c}\hfill D^{\nu }x\left(t\right)=\mathrm{C}{D}^{\nu }\mathbb{E}\left(t\right).\end{array}$$

(11)

To use MLG to solve Equation (1), subject to the initial condition (2), first substitute Equations (10) and (11) into Equation (1), to obtain the residual function $\mathbb{Re}\left(t\right)$:

$$\begin{array}{c}\hfill \Re \left(t\right)=\mathrm{C}{D}^{\nu }\mathbb{E}\left(t\right)+a\left(t\right){\left(\mathrm{C}\mathbb{E}\left(t\right)\right)}^2+b\left(t\right)\mathrm{C}\mathbb{E}\left(t\right)-g\left(t\right).\end{array}$$

(12)

We can obtain N nonlinear algebraic equations sets, by using the Galerkin method [49]:

$$\begin{array}{c}\hfill \underset{0}{\overset{1}{\int }}{\Re }_N\left(t\right)E_i^{\xi ,\eta }\left(t\right)\mathrm{d}t=0,i=0,\dots ,N-1.\end{array}$$

(13)

In addition, by inserting the initial condition (2) into Equation (10), we get

$$\begin{array}{c}\hfill x_{{}_N}\left(0\right)=\mathrm{C}\mathbb{E}\left(0\right)=x_0.\end{array}$$

(14)

The resulting combination of (13) and (14) yields $N+1$ of nonlinear equations with unknown coefficients ${\mathrm{c}}_0,{\mathrm{c}}_1,\dots ,{\mathrm{c}}_N$ that can be solved by the Newton method. As a result, the MLG solution to the $x_N\left(t\right)$ problem is obtained.

5. Error Analysis

This section investigates the error $e_{{}_N}\left(t\right)=x\left(t\right)-x_{{}_N}\left(t\right)$, where $x\left(t\right)$ is the exact solution and $x{}_N\left(t\right)$ is the MLG solution. We begin by defining the error by two theorems. The approach is then configured with the residual procedure, which produces estimates of the error, and new approximate solutions.

Theorem 1.

Let us approximate $x\left(t\right)\in C^{\infty }[0,1]$ as $x_{{}_N}\left(t\right)$, given in (10). Then, for every $t\in [0,1]$ there exists $\varpi \in [0,1]$, such that

$${∥x\left(t\right)-x_{{}_N}\left(t\right)∥}_{\infty }\le \frac{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}{(N+1)!}|E_{N+1}^{\xi ,\eta }\left(t\right)|\underset{\varpi \in [0,1]}{max}|x^{(N+1)}\left(\varpi \right)|.$$

(15)

Proof. 

Let $x\left(t\right)\in C^{\infty }[0,1]$ be approximated by $x_N\left(t\right)$, given in (10). Let us define the function:

$$L\left(t\right)=x\left(t\right)-x_{{}_N}\left(t\right)-\theta E_{N+1}^{\xi ,\eta }\left(t\right).$$

Let us select the parameter $\theta$, such that the equation $L\left(t\right)=0$ has a solution $t_0$, but $t_0$ is not a root of $E_{N+1}^{\xi ,\eta }\left(t\right)$, i.e., $E_{N+1}^{\xi ,\eta }\left(t_0\right)\ne 0$. Then, if we solve the equation $L\left(t_0\right)$, with respect to $\theta$, we get $x\left(t_0\right)-x_{{}_N}\left(t_0\right)-\theta E_{N+1}^{\xi ,\eta }\left(t_0\right)=0$, so that

$$\theta =\frac{x\left(t_0\right)-x_{{}_N}\left(t_0\right)}{E_{N+1}^{\xi ,\eta }\left(t_0\right)}.$$

(16)

Given $x\left(t\right)\in C^{\infty }[0,1],E_N^{\xi ,\eta }\left(t_0\right)\in C^N[0,\infty ]$ and $E_{N+1}^{\xi ,\eta }\left(t_0\right)\in C^{N+1}[0,\infty ]$, then $L\left(t\right)\in C^{N+1}[0,1]$ and $L^{(N+1)}\left(t\right)$ has at least one root in the interval; that is:

$$L^{(N+1)}\left(\varpi \right)=x^{(N+1)}\left(\varpi \right)-\theta {\left[E_{N+1}^{\xi ,\eta }\left(\varpi \right)\right]}^{(N+1)}-{\left[E_N^{\xi ,\eta }\left(\varpi \right)\right]}^{\left(N\right)}=0.$$

(17)

By using (9), the last term of (17), given ${\left[E_N^{\xi ,\eta }\left(\varpi \right)\right]}^{\left(N\right)}=0$ and

$$E_{N+1}^{\epsilon ,\eta }\left(\varpi \right)=\frac{{\varpi }^{N+1}}{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}+\mathrm{lower}-\mathrm{degree}\mathrm{polynomials},$$

then

$$\begin{array}{ccc}\hfill {\left[E_{N+1}^{\xi ,\eta }\left(\varpi \right)\right]}^{(N+1)}& =& \frac{(N+1)N(N-1)\times ...\times 3\times 2\times 1}{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}\hfill \\ & =& \frac{(N+1)!}{\mathsf{\Gamma }\left(\right(\xi (N+1))+\eta )}.\hfill \end{array}$$

Substituting ${\left[E_{N+1}^{\xi ,\eta }\left(\varpi \right)\right]}^{(N+1)}$ into (17) yields as follows:

$$\theta =\frac{\mathsf{\Gamma }\left(\right(\xi (N+1))+\eta )}{(N+1)!}{x}^{(N+1)}\left(\varpi \right).$$

(18)

Using Equations (16)–(18), we can write the following equation:

$$x\left(t_0\right)-x_{{}_N}\left(t_0\right)=\frac{\mathsf{\Gamma }\left(\right(\xi (N+1))+\eta )}{(N+1)!}{E}_{N+1}^{\xi ,\eta }\left(t_0\right)x^{(N+1)}\left(\varpi \right),$$

(19)

and so

$$|x\left(t_0\right)-x_{{}_N}\left(t_0\right)|=\frac{\mathsf{\Gamma }\left(\right(\xi (N+1))+\eta )}{(N+1)!}|E_{N+1}^{\epsilon ,\eta }\left(t_0\right)||x^{(N+1)}\left(\varpi \right)|.$$

Finally, taking the maximum of $|x^{(N+1)}\left(\varpi \right)|$ completes the proof. □

Theorem 2.

Let $x\left(t\right)\in C^{\infty }[0,1]$ be the solution, and $x_{{}_N}\left(t\right)$ be the MLG solution of (1)–(2), respectively, and let $p_N\left(t\right)={\sum }_{k=0}^N{c}_k^p{E}_k^{\xi ,\eta }\left(t\right)$. Then, for every $t\in [0,1]$, there exists $\varpi \in [0,1]$, such that the following inequality holds:

$$e_N\left(t\right)\le \frac{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}{(N+1)!}|E_{N+1}^{\xi ,\eta }\left(t\right)|\underset{\varpi \in [0,1]}{max}|p^{(N+1)}\left(\varpi \right)|+\sum _{k=0}^N|c_k^p-c_k^x||E_k^{\xi ,\eta }\left(t\right)|.$$

Proof. 

As $x_{{}_N}\left(t\right)$ is the MLG solution of (1)–(2), we can write the error as

$$\begin{array}{ccc}\hfill |x\left(t\right)-x_{{}_N}\left(t\right)|& =& |x\left(t\right)-p_N\left(t\right)+p_N\left(t\right)-x_{{}_N}\left(t\right)|\hfill \\ & \le & |x\left(t\right)-p_N\left(t\right)|+|p_N\left(t\right)-x_{{}_N}\left(t\right)|\hfill \\ & \le & \frac{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}{(N+1)!}|E_{N+1}^{\xi ,\eta }\left(t\right)|\underset{\varpi \in [0,1]}{max}|p^{(N+1)}\left(\varpi \right)|+|p_N\left(t\right)-x_{{}_N}\left(t\right)|\hfill \\ & \le & \frac{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}{(N+1)!}|E_{N+1}^{\xi ,\eta }\left(t\right)|\underset{\varpi \in [0,1]}{max}|p^{(N+1)}\left(\varpi \right)|+|\sum _{k=0}^N(c_k^p-c_k^x)E_k^{\xi ,\eta }\left(t\right)|\hfill \\ & \le & \frac{\mathsf{\Gamma }\left(\xi \right(N+1)+\eta )}{(N+1)!}|E_{N+1}^{\xi ,\eta }\left(t\right)|\underset{\varpi \in [0,1]}{max}|p^{(N+1)}\left(\varpi \right)|+\sum _{k=0}^N|c_k^p-c_k^x||E_k^{\xi ,\eta }\left(t\right)|.\hfill \end{array}$$

□

Now, we constitute the residual correction procedure.

Theorem 3.

Let $x\left(t\right)$ be the exact solution of (1)–(2), and let $x_N\left(t\right)$ be the MLG solution, respectively. The error $e_N\left(t\right)$ satisfies the following equation:

$$\left\{\begin{array}{cc}{D}^{\nu }{e}_N\left(t\right)=-a\left(t\right)e_N^2\left(t\right)-2a\left(t\right)e_N\left(t\right)x_N\left(t\right)-b\left(t\right)e_N\left(t\right)-{\Re }_N\left(t\right),\hfill & \\ e_N\left(0\right)=x_0-x_N\left(0\right).\hfill \end{array}\right.$$

(20)

Proof. 

We have

$${\Re }_N\left(t\right)=D^{\nu }{x}_N\left(t\right)+a\left(t\right)x_N^2\left(t\right)+b\left(t\right)x_N\left(t\right)-g\left(t\right),$$

and $D^{\nu }x\left(t\right)+a\left(t\right)x^2\left(t\right)+b\left(t\right)x\left(t\right)=g\left(t\right)$. Then:

$$\begin{array}{ccc}\hfill D^{\nu }{e}_N\left(t\right)& =& D^{\mu }x\left(t\right)-D^{\nu }{x}_N\left(t\right)\hfill \\ & =& -a\left(t\right)x^2\left(t\right)-b\left(t\right)x\left(t\right)+g\left(t\right)-{\Re }_N\left(t\right)+a\left(t\right)x_N^2\left(t\right)+b\left(t\right)x_N\left(t\right)-g\left(t\right)\hfill \\ & =& -a\left(t\right)x^2\left(t\right)-b\left(t\right)x\left(t\right)-{\Re }_N\left(t\right)+a\left(t\right)x_N^2\left(t\right)+b\left(t\right)x_N\left(t\right)\hfill \\ & =& -a\left(t\right)x^2\left(t\right)+a\left(t\right)x_N^2\left(t\right)-b\left(t\right)x\left(t\right)+b\left(t\right)x_N\left(t\right)-{\Re }_N\left(t\right)\hfill \\ & =& -a\left(t\right)(x^2\left(t\right)-x_N^2\left(t\right))-b\left(t\right)(x\left(t\right)-x_N\left(t\right))-{\Re }_N\left(t\right)\hfill \\ & =& -a\left(t\right)e_N\left(t\right)(x\left(t\right)+x_N\left(t\right))-b\left(t\right)e_N\left(t\right)-{\Re }_N\left(t\right)\hfill \\ & =& -a\left(t\right)e_N\left(t\right)(e_N\left(t\right)+2x_N\left(t\right))-b\left(t\right)e_{{}_N}\left(t\right)-{\Re }_{{}_N}\left(t\right)\hfill \\ & =& -a\left(t\right)e_{{}_N}^2\left(t\right)-2a\left(t\right)e_{{}_N}\left(t\right)x_{{}_N}\left(t\right)-b\left(t\right)e_{{}_N}\left(t\right)-{\Re }_{{}_N}\left(t\right).\hfill \end{array}$$

□

By solving the error problem (20) by the MLG method, we get the approximation for the error, i.e.,

$$e_{N,M}\left(t\right)=\sum _{i=0}^M{c}_i^e{E}_i^{\xi ,\eta }\left(t\right),$$

where the coefficients $c_i^e,i=0,1,2,\dots ,M$ are unknown constants. Hence, the error can be estimated by using $e_{N,M}\left(t\right)$ in the case of $∥e_{{}_N}\left(t\right)-e_{N,M}\left(t\right)∥<\epsilon .$ On the other hand, $x_{{}_N}\left(t\right)+e_{N,M}\left(t\right)$ is also an approximate solution of (1)–(2). We call the solution $x_{{}_N}\left(t\right)+e_{N,M}\left(t\right)$, as the corrected MLG solution. The corrected MLG solution is a better approximation than the MLG solution in any given norm, whenever

$$∥e_{{}_N}\left(t\right)-e_{N,M}\left(t\right)∥<∥x\left(t\right)-x_{{}_N}\left(t\right)∥.$$

6. Numerical Experiments

In this section, we present some numerical experiments to show how the method is applied, to show the efficiency of the method—by giving the precision of the results obtained using the method—and to compare the method to other methods. All the experiments were performed using Maple on a laptop with an Intel Core i3 processor and 4GB of RAM.

Example 1.

As a first example, let us apply the method to the following FDRE, whose exact solution is $x\left(t\right)=t^2$ [35,37]:

$$\begin{array}{c}{D}^{\frac{1}{2}}x\left(t\right)+x\left(t\right)+x^2\left(t\right)=\frac{8}{3\sqrt{\pi }}{t}^{\frac{3}{2}}+t^2+t^4,0\le t\le 1,\end{array}$$

(21)

with the condition

$$x\left(0\right)=0.$$

(22)

For $N=2$, the MLG solution can be written as

$$\begin{array}{ccc}\hfill x_2\left(t\right)& =& \sum _{k=0}^2{c}_k{E}_k^{\xi ,\eta }\left(t\right)\hfill \\ & =& c_0{E}_0^{\xi ,\eta }\left(t\right)+c_1{E}_1^{\xi ,\eta }\left(t\right)+c_2{E}_2^{\xi ,\eta }\left(t\right).\hfill \end{array}$$

From $\xi =\frac{1}{2},\eta =1$, Equation (12) gives Equation (A1) (see Appendix A). In addition, we have, from Equation (22),

$$c_0+c_1+c_2=0.$$

(23)

Finally, by solving Equations (23) and (A1) (see Appendix A), we obtain

$$c_0=-5.07096009275\times {10}^{-13},c_1=-1.00000000000\mathit{and}{c}_2=1.00000000000.$$

Thus,

$$\begin{array}{ccc}\hfill x_2\left(t\right)& =& c_0{E}_0^{\frac{1}{2},1}\left(t\right)+c_1{E}_1^{\frac{1}{2},1}\left(t\right)+c_2{E}_2^{\frac{1}{2},1}\left(t\right)\hfill \\ & =& c_0+c_1(1-\frac{2t}{\sqrt{\pi }})+c_2(1-\frac{2t}{\sqrt{\pi }}+t^2)\hfill \\ & =& -5.07096009275\times {10}^{-13}+1.00000000000t^2\hfill \\ & \simeq & t^2.\hfill \end{array}$$

From Equation (23), we construct the error equation

$$\left\{\begin{array}{cc}{D}^{\frac{1}{2}}(e_2\left(t\right)+x_2\left(t\right))+e_2\left(t\right)+x\left(t\right)+{(e_2\left(t\right)+x\left(t\right))}^2=\frac{8}{3\sqrt{\pi }}{t}^{\frac{3}{2}}+t^2+t^4,\hfill & \\ e_2\left(0\right)=x_2\left(0\right).\hfill \end{array}\right.$$

(24)

Let us solve error Equation (24) by using the MLG method for $M=3$. The MLG solution of the error equation is obtained as

$$\begin{array}{ccc}\hfill e_{2,3}\left(t\right)& =& 4\times {10}^{-198}-\frac{9.0414\times {10}^{-99}}{\sqrt{\pi }}t+6.6595\times {10}^{-99}{t}^2-\frac{3.9160\times {10}^{-100}}{\sqrt{\pi }}{t}^3.\hfill \end{array}$$

Adding this solution to the MLG solution yields the corrected MLG solution. We present the error, with its estimations and the corrected error, in Figure 1 for $(N,M)=(2,3)$. We can say, from the figures, that the error estimation obtained by the procedure well fits the error.

The development of the basic definition, Equation (8), was realized this way, and affected the convergence of solutions; we added a comparison between the development and the original definition. Now, if we use

$$\begin{array}{c}\hfill E_i^{\xi ,\eta }\left(t\right)=\sum _{k=0}^i\frac{t^k}{\mathsf{\Gamma }(\xi k+\eta )},\end{array}$$

we have

$$c_0=-5.10444760178,c_1=16.9410134078,c_2=-11.8365658060,$$

and then

$$\begin{array}{ccc}\hfill x_2\left(t\right)& =& \frac{10.2088952036t}{\sqrt{\pi }}-11.8365658060t^2\hfill \\ & \ne & t^2.\hfill \end{array}$$

It appeared to us that the original definition had a solution that was not convergent; however, in the developed definition, the solution was convergent.

Example 2.

Let us consider the following FDRE:

$$\begin{array}{c}{D}^{\frac{1}{2}}x\left(t\right)-t{x}^2\left(t\right)=\frac{16}{5}\frac{t^{\frac{5}{2}}}{\sqrt{\pi }}-t^7,0\le t\le 1,\end{array}$$

(25)

$$x\left(0\right)=0.$$

(26)

The exact solution to the problem is

$$x\left(t\right)=t^3.$$

Let us find the MLG solution to the problem for $N=3$, which can be written as follows:

$$\begin{array}{ccc}\hfill x_3\left(t\right)& =& \sum _{k=0}^3{c}_k{E}_k^{\xi ,\eta }\left(t\right)\hfill \\ & =& c_0{E}_0^{\xi ,\eta }\left(t\right)+c_1{E}_1^{\xi ,\eta }\left(t\right)+c_2{E}_2^{\xi ,\eta }\left(t\right)+c_3{E}_3^{\xi ,\eta }\left(t\right).\hfill \end{array}$$

From $\xi =\frac{1}{2},\eta =1$, (12) gives Equation (A2) (see Appendix A). In addition, we have, from Equation (26):

$$c_0+c_1+c_2+c_3=0.$$

(27)

Finally, by solving Equations (27)–(A2), we obtain

$$c_0=1.1901312585\times {10}^{-9},c_1=3.0458221798\times {10}^{-9},c_2=1.3293403878,$$

and $c_3=-1.329340392$. Thus, we obtain the MLG solution for the problem and $N=3$ as

$$\begin{array}{ccc}\hfill x_3\left(t\right)& =& -\frac{2.39\times {10}^{-9}}{\sqrt{\pi }}t+4.24\times {10}^{-9}{t}^2+\frac{1.77245385607}{\sqrt{\pi }}{t}^3\hfill \\ & \simeq & t^3.\hfill \end{array}$$

From Equation (20), we construct the error equation

$$\left\{\begin{array}{cc}{D}^{\frac{1}{2}}{e}_3\left(t\right)=t{e}_3^2\left(t\right)+2t{e}_3\left(t\right)x_3\left(t\right)-{\Re }_3\left(t\right),\hfill & \\ e_3\left(0\right)=-x_3\left(0\right),\hfill \end{array}\right.$$

(28)

where

$$\begin{array}{ccc}\hfill {\Re }_3\left(t\right)& =& D^{\frac{1}{2}}{x}_3\left(t\right)-t{x}_3^2\left(t\right)-(\frac{16}{5}\frac{t^{\frac{5}{2}}}{\sqrt{\pi }}-t^7)\hfill \\ & =& 2.39449496167\times {10}^{-13}\sqrt{t}(7.53982239058\times {10}^{12}{t}^2-26640.7057024t\hfill \\ & & +6354.24705550)\hfill \\ & & -t(\frac{2.39000000000\times {10}^{-9}}{\sqrt{\pi }}t-4.2400000000\times {10}^{-9}{t}^2\hfill \\ & & +\frac{1.77245385607}{\sqrt{\pi }}{t}^3{)}^2-\frac{16}{5\sqrt{\pi }}{t}^{\frac{5}{2}}+t^7.\hfill \end{array}$$

By solving the error problem (28) for $M=4$, we obtain the estimation of the error as

$$\begin{array}{ccc}\hfill e_{3,4}(t)& =& -3\times {10}^{-21}-\frac{2.40732926828\times {10}^{-9}}{\sqrt{\pi }}t+4.30060806282\times {10}^{-9}{t}^2\hfill \\ & & -\frac{5.35088552742\times {10}^{-9}}{\sqrt{\pi }}{t}^3+5.4884232558\times {10}^{-11}{t}^4.\hfill \end{array}$$

For $(N,M)=(3,4)$, the error and the error estimation of Example 2 and the corrected error results are drawn in Figure 2.

In

Table 1. The effects of the Mittag-Leffler parameters $\xi$ and $\eta$ on the error for Example 2 and $N=3$.

$\mathit{\xi }$	$\mathit{\eta }$	Absolute Error
$0.0$	$0.1$	$5.782602642\times {10}^{-10}$
	$0.2$	$1.840492393\times {10}^{-9}$
	$0.3$	$1.900000000\times {10}^{-9}$
	$0.4$	$4.413182132\times {10}^{-9}$
	$0.5$	$1.560493121\times {10}^{-9}$
$0.5$	$0.9$	$1.853612308\times {10}^{-11}$
	1	$1.934584037\times {10}^{-11}$
	2	$1.150677936\times {10}^{-10}$

, we present the effects of the Mittag-Leffler parameters ξ and η on the error, which shows that our results are more accurate:

Example 3.

Consider the following FDRE [43,45]:

$$\begin{array}{c}{D}^{\nu }x\left(t\right)+x^2\left(t\right)-1=0,0\le t\le 1,\end{array}$$

(29)

$$x\left(0\right)=0.$$

(30)

The exact solution of the problem is $x\left(t\right)=tanh\left(t\right)$ for $\nu =1$.

We apply the MLG method to the problem for $\nu =1$, $\xi =\frac{1}{2},\eta =1$. The MLG solutions for $N=4$ and $N=8$ are obtained as follows:

$$\begin{array}{ccc}\hfill x_4\left(t\right)& =& -1\times {10}^{-10}+0.9998413015t+0.0081494057t^2\hfill \\ & & -0.3915419735t^3+0.1451458340t^4,\hfill \\ \\ \hfill x_8\left(x\right)& =& 0.999794070t+0.5516765e-2t^2+0.1917058000t^4+0.4684949362t^6\hfill \\ & & +0.0897t^8-0.381095428t^3-0.2724512057t^5-0.3400392225t^7+{10}^{-9}.\hfill \end{array}$$

By applying the method to the error equations obtained by $x_4$ and $x_8$, the estimations of $e_4$ and $e_8$ are found for $M=7$ and $M=9$, respectively, as follows:

$$\begin{array}{ccc}\hfill e_{4,7}\left(t\right)& =& c_0^e{E}_0^{\frac{1}{2},1}\left(t\right)+c_1^e{E}_1^{\frac{1}{2},1}\left(t\right)+c_2^e{E}_2^{\frac{1}{2},1}\left(t\right)+c_3^e{E}_3^{\frac{1}{2},1}\left(t\right)+c_4^e{E}_4^{\frac{1}{2},1}\left(t\right)+c_5^e{E}_5^{\frac{1}{2},1}\left(t\right)\hfill \\ & & +c_6^e{E}_6^{\frac{1}{2},1}\left(t\right)+c_7^e{E}_7^{\frac{1}{2},1}\left(t\right)\hfill \\ & =& 1.8629\times {10}^{-4}t-0.1208t^6-0.0090t^2+0.0663t^3-0.1822t^4+0.2240t^5+0.0214t^7,\hfill \\ \\ \hfill e_{8,9}\left(t\right)& =& -7.3650\times {10}^{-6}t+2.6532e-04t^2+0.0171t^4+0.0892t^6+0.0480t^8-0.0031t^3\hfill \\ & & -0.0514t^5-0.0894t^7-0.0107t^9-4\times {10}^{-100},\hfill \end{array}$$

The results are given in Figure 3. As a result, we can say that the absolute errors $e_4$ and $e_8$ are estimated by $e_{4,7}$ and $e_{8,9}$, respectively. On the other hand, for each case, the corrected MLG solutions are better than the MLG solutions.

Table 2. Comparison of the approximate solutions of $x\left(t\right)$ with the MLG method to the methods in [43,45,50] schemes, with $\nu =1$ for Example 3.

t	Exact Solution	Method [43]	Method [45]	Method [50]	MLG Method
$0.0$	$0.00000000000$	$0.000000000000$	$0.0000000011$	$0.000000000000$	$0.00000000000$
$0.2$	$0.19737532022$	$0.197375320493$	$0.1973918880$	$0.19737532017$	$0.19737532019$
$0.4$	$0.37994896226$	$0.379948962506$	$0.3799632287$	$0.379948962207$	$0.37994896222$
$0.6$	$0.53704956700$	$0.537049567214$	$0.5370622335$	$0.53704956701$	$0.53704956700$
$0.8$	$0.66403677027$	$0.664036770562$	$0.6640456511$	$0.66403677030$	$0.66403677029$
$1.0$	$0.76159415596$	$0.761594224400$	$0.7616019763$	$0.76159415595$	$0.76159415595$

presents the comparison of the approximate solution of $x\left(t\right)$ for $\nu =1,N=12$ to the Bernoulli wavelet operational matrix [50], the computational intelligence approach [45] and the Bernstein collocation method [43]. For $\nu =1$, the results were compared to the operational matrices method [51], the differential squared method [52] method, the method dependent on shifted Chebyshev polynomials [53] and the hybrid functions approach [54] in

Table 3. Comparison of the errors obtained by the MLG method to the methods in [51,52,53,54] schemes, with $\nu =1$ for Example $3$.

t	Method [51]	Method [54]	Method [52]	Method [53]	MLG Method
	$N=12$	$N=20$	$N=12$	$N=12$	$N=12$
$0.1$	$1.11\times {10}^{-10}$	$7.2701\times {10}^{-6}$	$8.3141\times {10}^{-11}$	$3.0881\times {10}^{-11}$	$2.8244\times {10}^{-11}$
$0.2$	$2.04\times {10}^{-10}$	$1.0922\times {10}^{-5}$	$9.1576\times {10}^{-11}$	$4.8024\times {10}^{-11}$	$3.4678\times {10}^{-11}$
$0.3$	$2.10\times {10}^{-12}$	$1.3476\times {10}^{-5}$	$7.5812\times {10}^{-11}$	$4.6345\times {10}^{-11}$	$3.5020\times {10}^{-11}$
$0.4$	$2.23\times {10}^{-10}$	$1.4755\times {10}^{-5}$	$1.1151\times {10}^{-10}$	$4.8119\times {10}^{-11}$	$2.7373\times {10}^{-11}$
$0.5$	$4.03\times {10}^{-10}$	$1.4778\times {10}^{-5}$	$5.5890\times {10}^{-11}$	$9.7303\times {10}^{-12}$	$9.2649\times {10}^{-12}$
$0.6$	$1.79\times {10}^{-10}$	$1.3730\times {10}^{-5}$	$7.8642\times {10}^{-11}$	$1.7617\times {10}^{-11}$	$1.0959\times {10}^{-11}$
$0.7$	$8.59\times {10}^{-11}$	$1.1891\times {10}^{-5}$	$6.2746\times {10}^{-11}$	$4.1180\times {10}^{-11}$	$2.2709\times {10}^{-11}$
$0.8$	$2.70\times {10}^{-10}$	$9.5751\times {10}^{-6}$	$5.3920\times {10}^{-11}$	$3.2346\times {10}^{-11}$	$2.3907\times {10}^{-11}$
$0.9$	$1.89\times {10}^{-10}$	$7.0732\times {10}^{-6}$	$4.8389\times {10}^{-11}$	$2.6919\times {10}^{-11}$	$1.5571\times {10}^{-11}$

. Figure 4 shows the results of the MLG method for $\nu =1$ and different values of N. The approximate values of $x\left(t\right)$ for $N=10$ and $\nu =\{0.75,0.9,0.95,1\}$ are given in Figure 5. Although the results of the MLG method are better than the results of the hybrid functions approach [54], we can say from Table 2 and Table 3 that it yields approximation results similar to the other methods for this problem. We conclude from Figure 4 that increasing $=N$ yields better approximation results for the problem.

Figure 6 for $N=12$, the absolute error for Example 3.

Example 4.

Let us consider the following FDRE [55]:

$$\begin{array}{c}{D}^{\nu }x\left(t\right)-x^2\left(t\right)-1=0,0\le t\le 1,\end{array}$$

(31)

$$x\left(0\right)=0.$$

(32)

which has the exact solution for $\nu =1$ is $x\left(t\right)=tan\left(t\right)$.

The results obtained when the method is applied for $\nu =1$ and $N=\{4,8,12\}$ are given below:

$$\begin{array}{ccc}\hfill x_4\left(t\right)& =& 0.93859t+0.53891t^2-0.98546t^3+1.06529t^4;\hfill \\ \\ \hfill x_8\left(t\right)& =& 0.99943t+0.0192t^2+1.04658t^4+4.07934t^6+1.03275t^8+0.12481t^3\hfill \\ & & -2.65049t^5-3.09424t^7-{10}^{-13};\hfill \\ \\ \hfill x_{12}\left(t\right)& =& 0.99998t+8.9376\times {10}^{-04}{t}^2+0.0986t^4+0.79925t^6+0.28047t^8-0.0774t^{10}\hfill \\ & & +0.12974t^{12}+0.31956t^3-0.24372t^5-0.81591t^7+0.27357t^9-0.20760t^{11}.\hfill \end{array}$$

The solution for N = 12 and $\nu =1$ is given in

Table 4. Comparison of the approximate solution of Example 4, using the MLG method, to the solutions obtained by the methods presented in [44,53,55] for $\nu =1$:

t	Exact Solution	Method [44]	Method [55]	Method [53]	Our Method
				$N=14$	$N=12$
$0.0$	$0.00000000000$		$0.0000000000$	$0.0000000000$	$0.00000000000$
$0.1$	$0.10033467208$	$0.100342$	$0.1003346713$	$0.1003346714$	$0.10033465034$
$0.2$	$0.20271003550$	$0.202726$	$0.2027099297$	$0.2027100349$	$0.20271006703$
$0.3$	$0.30933624961$	$0.309372$	$0.3093343442$	$0.3093362509$	$0.30933621959$
$0.4$	$0.42279321873$	$0.422832$	$0.4227777155$	$0.4227932186$	$0.42279323380$
$0.5$	$0.54630248984$	$0.546363$	$0.5462212762$	$0.5463024891$	$0.68413677385$
$0.6$	$0.68413680834$	$0.684251$	$0.6838056920$	$0.6841368110$	$0.68413677385$
$0.7$	$0.84228838046$	$0.842411$	$0.8411449022$	$0.8422883779$	$0.84228842425$
$0.8$	$1.02963855705$	$1.029849$	$1.0261001110$	$1.0296385599$	$1.02963851386$
$0.9$	$1.26015821755$	$1.260573$	$1.2499664940$	$1.2601582184$	$1.26015825498$
$1.0$	$1.55740772465$	$1.557938$	$1.5293009690$	$1.5574077258$	$1.55740772465$

and

Table 5. Comparison of the absolute errors obtained using the MLG method for $\nu =1$, in Example 4, to the absolute errors obtained using the methods in [53,55]:

t	Method [55]	Method [53]	Our Method
		$N=12$	$N=12$
$0.0$	$0.00000000$	$0.00000000$	$0.00000000$
$0.1$	$8.162\times {10}^{-10}$	$2.8897\times {10}^{-8}$	$2.1736\times {10}^{-8}$
$0.2$	$1.0580\times {10}^{-7}$	$5.1478\times {10}^{-8}$	$3.1526\times {10}^{-8}$
$0.3$	$1.9050\times {10}^{-6}$	$4.6086\times {10}^{-8}$	$3.0013\times {10}^{-8}$
$0.4$	$1.5500\times {10}^{-5}$	$3.4828\times {10}^{-8}$	$1.5069\times {10}^{-8}$
$0.5$	$8.1210\times {10}^{-5}$	$2.3389\times {10}^{-8}$	$1.1149\times {10}^{-8}$
$0.6$	$3.3110\times {10}^{-4}$	$5.0755\times {10}^{-8}$	$3.4490\times {10}^{-8}$
$0.7$	$1.1430\times {10}^{-3}$	$7.3355\times {10}^{-8}$	$4.3796\times {10}^{-8}$
$0.8$	$3.5380\times {10}^{-3}$	$4.4578\times {10}^{-8}$	$4.3187\times {10}^{-8}$
$0.9$	$1.1090\times {10}^{-3}$	$5.8748\times {10}^{-8}$	$3.7439\times {10}^{-8}$
$1.0$	$2.8110\times {10}^{-3}$	$2.2418\times {10}^{-8}$	$1.7293\times {10}^{-15}$

, by comparing the results obtained to some other numerical methods, such as the wavelet operational matrix method [44], the method dependent on shifted Chebyshev polynomials [53] and the decomposition algorithm [55]. We can conclude, from the tables, that the MLG method produces results similar to the other methods cited, except for the decomposition algorithm [55]. Figure 7 depicts the function Log${}_{10}$ of the absolute error for $\nu =1$ and various values of N, whereas Figure 8 shows the approximate values of $x\left(t\right)$ for $N=10$ and some values of ν. In conclusion, we can say that increasing N gives better approximations for this problem. The errors are given in Figure 9 and Figure 10 for $N=12$ and $N=10$, respectively. In Figure 10, the estimation of absolute error and the corrected absolute error for Example 4 are given. We can say from Figure 9 and Figure 10 that the MLG method gives more accurate results, and that the residual correction procedure estimates the error well.

Figure 9 for $N=12$, the absolute error for Example 4.

Example 5.

As a final example, let us consider the following FDRE [43,53]:

$$\begin{array}{c}{D}^{\nu }x\left(t\right)-2x\left(t\right)+x^2\left(t\right)-1=0,0\le t\le 1,\end{array}$$

(33)

$$x\left(0\right)=0.$$

(34)

which has the exact solution for $\nu =1$ is $x\left(t\right)=1-\sqrt{2}\frac{\sqrt{2}tanh\left(\sqrt{2}t\right)-1}{tanh\left(\sqrt{2}t\right)-\sqrt{2}}$.

Applying the MLG method to the problem for $\nu =1$ and $N=\{4,8,12\}$ yields the following MLG solutions:

$$\begin{array}{ccc}\hfill x_4\left(t\right)& =& 0.9655t+1.26574t^2-0.14446t^3-0.39736t^4;\hfill \\ \\ \hfill x_8\left(t\right)& =& 0.99991t+1.00259t^2-0.26776t^4-0.58544t^6-0.35438t^8+0.31134t^3\hfill \\ & & -0.46069t^5+1.04393t^7+4\times {10}^{-91};\hfill \\ \\ \hfill x_{12}\left(t\right)& =& t+0.99999t^2-0.33748t^4-0.36952t^6-1.60334t^8-2.81446t^{10}-0.234263t^{12}\hfill \\ & & +0.33358t^3-0.42847t^5+0.93948t^7+2.90224t^9+1.30175t^{11}.\hfill \end{array}$$

The results for $\nu =1,N=12$ are given in

Table 6. Comparison of the MLG solution to the methods in [43,53] for $\nu =1$ and Example $5$.

t	Exact Solution	Method [43]	Method [53]	Our Method
			$N=14$	$N=12$
$0.0$	$0.0000000000$	$0.000000000$	$0.0000000000$	$0.0000000000$
$0.2$	$0.2419767996$	$0.241977035$	$0.2419767995$	$0.24197679964$
$0.4$	$0.5678121663$	$0.567812472$	$0.5678121662$	$0.56781216535$
$0.6$	$0.9535662164$	$0.953566555$	$0.9535662164$	$0.95356621379$
$0.8$	$1.3463636554$	$1.346363997$	$1.3463636552$	$1.34636365269$
$1.0$	$1.6894983916$	$1.689510190$	$1.6894983916$	$1.68949839159$

and Figure 11 with a comparison to the approximate solutions obtained by the methods in [43,53]. In addition, Figure 12 displays Log${}_{10}$ of absolute errors for $\nu =1$ and different values of N. We conclude that, as in the other examples, more accurate results are obtained by increasing N. The approximate values of $x\left(t\right)$ for $N=10$, and some values of ν, are given in Figure 13. The absolute error, with its estimation, and the corrected absolute error for Example 5 are drawn in Figure 14 for $N=4$. The estimation absolute error for Example 5 are drawn in Figure 15 for $N=5$ and $N=8$.

The procedure estimates the error well, and leads to a better approximate solution. We also compared the absolute errors estimated by the MLG method to the Bernstein collocation method [43] for $\nu =0.9$, and we give the results in

Table 7. Comparison of the estimated absolute errors for various values of N and M for $\nu =0.9$ of Equation (33) to $\xi =0.5,\eta =1$.

t	$(\mathit{N},\mathit{M})=(5,7)$		$(\mathit{N},\mathit{M})=(8,9)$
	Method [43]	Our Method	Method [43]	Our Method
0	0	0	0	0
$0.2$	$2.4585\times {10}^{-3}$	$1.8653\times {10}^{-3}$	$5.5413\times {10}^{-5}$	$9.2910\times {10}^{-6}$
$0.4$	$2.3010\times {10}^{-3}$	$1.4361\times {10}^{-3}$	$6.6222\times {10}^{-5}$	$1.1642\times {10}^{-5}$
$0.6$	$2.5766\times {10}^{-3}$	$2.3933\times {10}^{-4}$	$6.5671\times {10}^{-5}$	$6.2586\times {10}^{-6}$
$0.8$	$1.6519\times {10}^{-3}$	$4.8353\times {10}^{-4}$	$5.6089\times {10}^{-5}$	$9.4982\times {10}^{-6}$
$1.0$	$6.3526\times {10}^{-3}$	$1.1605\times {10}^{-4}$	$6.4831\times {10}^{-4}$	$4.2162\times {10}^{-6}$

.

7. Conclusions

In this study, we introduced a new numerical method for solving fractional differential Riccati equations. The method is based on the Mittag-Leffler function and the Galerkin method. Some theorems related to the error analysis of the method were presented. The error can be bounded by these theorems. The residual correction method allows for estimating the absolute error and obtaining a new approximate solution, which was presented for the method. We applied the method to some test examples, to illustrate its effectiveness and how the method is applied. The results obtained by the method were given by comparing the results of some other previously known methods for the solutions of fractional Riccati equations. The numerical test results showed that the method gave good approximation results for the examples. In addition, while only similar results were obtained by some of the known methods, more accurate results were obtained for the problems under consideration than were obtained by others of the known methods.