#### Abstract

The Riccati differential equation is a well-known nonlinear differential equation and has different applications in engineering and science domains, such as robust stabilization, stochastic realization theory, network synthesis, and optimal control, and in financial mathematics. In this study, we aim to approximate the solution of a fractional Riccati equation of order with Atangana–Baleanu derivative (ABC). Our numerical scheme is based on Laplace transform (LT) and quadrature rule. We apply LT to the given fractional differential equation, which reduces it to an algebraic equation. The reduced equation is solved for the unknown in LT space. The solution of the original problem is retrieved by representing it as a Bromwich integral in the complex plane along a smooth curve. The Bromwich integral is approximated using the trapezoidal rule. Some numerical experiments are performed to validate our numerical scheme.

#### 1. Introduction

In applied mathematics, the fractional differential equation (FDE) is an equation which contains derivatives of arbitrary order. Fractional derivatives first appeared in 1695 [1]. Recently, the research community took more interest in fractional calculus because of its applications in engineering and other sciences [2–4]. Many physical systems show fractional order behavior that may vary with time or space. Fractional derivatives have many kinds. The three important and most commonly used fractional derivatives are Grnwald–Letnikov derivative, Riemann–Liouville fractional derivative, and Caputo derivative [5]. However, these classical fractional derivatives have a singular kernel, and hence, they may face difficulties in describing the nonlocality of real-world dynamics. In order to handle the nonlocal systems in a better way, recently, new fractional derivatives with nonsingular kernels are defined such as the Caputo–Fabrizio (CF) derivative and Atangana–Baleanu (ABC) derivative [6, 7]. Fractional derivatives with nonsingular kernels become more valuable due to the fact that numerous phenomena cannot properly be modeled by fractional derivatives with singular kernels [8].

In this work, we aim to approximate a Riccati differential equation (RDE) with ABC derivative. RDEs have many applications such as random processes, optimal control, and diffusion process [9]. The RDE of fractional order has been studied by many authors; for example, in [10], the authors developed the Adomain decomposition method for the solution of RDE of fractional order. In [11], some analytic techniques are presented for the solution of RDE. The authors [12] obtained the solution of RDE using the differential transform method. In [13], the authors have developed a Laplace-transform Adomain decomposition method for the solution of RDE. Other works on the analytic solution of RDE can be found in [14, 15] and the references therein.

Most of the time, the exact/analytical solution of FDEs cannot be found, so numerical approximations must be utilized [16–19]. Numerous numerical methods have been developed for numerical approximation of FDEs such as the Chebyshev collocation method [20], variation iteration method [16], reproducing kernel Hilbert space method [18, 19], and homotopy perturbation method ([21] and references therein). The numerical solutions of fractional-order RDE have been studied by a large number of researchers. For example, the authors [20] studied fractional-order RDE with ABC derivatives, and they established the existence and uniqueness results using Banach fixed-point theorem. A reproducing kernel Hilbert space method [19] for approximating RDEs and Bernoulli differential equations of fractional order with ABC derivative has been presented. In [18], the authors developed an iterative reproducing kernel Hilbert space method for numerical approximation of fractional RDE. The authors [22] studied the numerical solution of fractional-order RDE using the modified homotopy perturbation method. Numerical solution of fractional-order RDEs using Bernstein polynomials is considered [23]. A fractional Chebyshev finite difference method [24] for numerical investigation of RDE of fractional order is proposed. Jafari and Tajadodi [25] proposed a variational iteration method for solving the RDE of fractional order. The authors in [26] proposed a hybrid method based on the quasilinearization technique and reproducing kernel method for RDE of fractional order. In [27], the solution of fractional-order RDE is approximated using fractional polynomial approximations. The authors in [28] developed a method based on finite difference and Padé-variational iteration method for solving the RDE of noninteger order. A numerical method based on the path following methods and the Tau Legendre is presented [29] for the solution of fractional-order RDE. The authors of [30] proposed a modified variational iteration method based on Adomain polynomials for the solution of RDE. Khashan et al. [31] have utilized the Haar wavelets for the approximation of fractional-order RDE. The authors of [32] utilized the finite difference scheme for the approximation of RDE. Other works on the numerical approximation of the solution of fractional-order RDE can be found in [33–42] and references therein. In this work, we approximate the solution of fractional-order RDE with ABC derivative of the following form:here, denotes Atangana–Baleanu fractional derivative of order .

#### 2. Preliminaries

In this section, we present some basic results and definitions. For details about fractional calculus, we refer [43–49], and for details about modeling, we refer [50, 51] and references therein.

*Definition 1. *The Mittag–Leffler (ML) function with one parameter is defined as [52, 53]

*Definition 2. *The two-parameter Mittag–Leffler (ML) function is defined as [53, 54]

*Definition 3. *The ABC fractional derivative of order of with base point at is defined as [7]where is a first-order Sobolev space equipped with -norm over the region , which is defined asand the term is given as

*Definition 4. *The (ABC) fractional integral of order of with base point at is defined as [7]

*Definition 5. *The Laplace transform of a piecewise continuous function is defined as

*Definition 6. *The LT of function with one parameter is defined as [4]

*Definition 7. *The LT of two-parameter function is defined as [4]

*Definition 8. *If , then the LT of the ABC derivative is defined by [7]

#### 3. Laplace Transform Method for Fractional Riccati Equation with ABC Derivative

In this section, we give a detailed description of our proposed numerical scheme for approximating the inverse Laplace transform. First, we will apply the Laplace transform to the given fractional problem, which will transform it to an algebraic equation. After solving the reduced equation, the solution of the original problem can be obtained by representing it as a contour integral in the left half of the complex plane. The trapezoidal rule is then utilized to approximate the contour integral. Applying the Laplace transform to equation (1), we obtainwhich can be written in simplified form aswhere

In our method, first we represent the solution of original problem (1) as a contour integral:where, for , is appropriately large and is an initially appropriately chosen line perpendicular to the real axis in the complex plane, with . Integral (15) is just the inverse transform of , with the condition that it must be analytic to the right of . To make sure the contour of integration remains in the domain of analyticity of , we select as a deformed contour in the set , which behaves as a pair of asymptotes in the left half plane, with when , which force to decay towards both ends of . In our work, we choose aswhere

By writing , we notice that (16) is the left branch of the following hyperbola:The asymptotes for (18) are and *x*-intercept at . Condition (17) confirms that lies in the sector and grows into the left half plane. From (16) and (15), we obtain

The trapezoidal rule is used for the approximation of equation (19) with step as follows:where .

#### 4. Error Analysis

In the process of obtaining the solution of problem equations (1)–(13), the fractional integrodifferential equation is first transformed to an algebraic equation using Laplace transform, and this causes no error. The transformed equation is then solved for the unknown in the Laplace space. Finally, the solution is obtained using inverse Laplace transform via integral representation (19). The integral is then approximated using quadrature rule. In the process of approximating integral (19), convergence is achieved at different rates depending on the path . In approximating integral (19), the convergence order relies on the step of the quadrature rule and the time domain . The proof for the order of quadrature error is given in the next theorem.

Theorem 1 ([54] Theorem 2.1). *Let be the solution of (1) with being analytic in . Let , and define by , where , , and and let . Then, for equation (20), with , we have , for , , , , , and . Hence, the error estimate for the proposed scheme is*

#### 5. Results and Discussion

Most of the time analytical methods cannot be applied to handle a real-world problem. So we need numerical methods to approximate the solutions of the problems. In this section, we consider fractional-order Riccati equations to validate our method.

*Problem 1. *Here, we consider the fractional Riccati equation as follows:with .

The problem has exact solution . In this experiment, the optimal parameters utilized are . The quadrature nodes are generated using the MATLAB commands . The values of approximate solution for different fractional orders and are depicted in Table 1. The absolute errors for various quadrature nodes and fractional order are depicted in Table 2. The method produced almost exact values for different values of fractional order . The plots of absolute error for fractional orders are displayed in Figure 1(a), where Figure 1(b) shows the comparison between the absolute error and error estimate for the fractional order . Figures 2(a) and 2(b) show the error functions for fractional orders and , respectively. It can be seen that the method can solve fractional Riccati equation with ABC derivative efficiently.

**(a)**

**(b)**

**(a)**

**(b)**

*Problem 2. *Here, we consider the fractional Riccati equation as follows:with .

The problem has exact solution . In this experiment, the optimal parameters utilized are . The quadrature nodes are generated using the MATLAB commands . The values of approximate solution for different fractional orders and are depicted in Table 3. The method produced exact values for . The plots of numerical numerical solutions for different fractional orders are displayed in Figure 3, where Figure 4(a) shows the absolute error for the fractional order and Figure 4(b) shows the error function for .

**(a)**

**(b)**

*Problem 3. *Here, we consider the fractional Riccati equation as follows:withwhere,The problem has exact solution . In this experiment, the same set of optimal parameters is utilized. The values of approximate solution for different fractional orders and are depicted in Table 4. The absolute errors for various values of and fractional order are shown in Table 5. Figure 5 shows the comparison between the absolute error and error_est, and a good agreement between them is observed. The plots of numerical and exact solutions for different fractional orders are displayed in Figures 6–9. It can be seen that the proposed method has produced good results and the exact and numerical solutions are in good agreement. This shows that this method can solve fractional Riccati equation with ABC derivative efficiently.

*Problem 4. *Here, we consider the fractional Riccati equation as follows:with

The problem has exact solution . In this experiment, the same set of optimal parameters is utilized. The values of approximate solution for different fractional orders and are depicted in Table 6. Figure 10(a) displays the comparison between the absolute error and error estimate. Figure 10(b) shows the error function for .

**(a)**

**(b)**

*Problem 5. *Here, we consider the fractional Riccati equation as follows:with

The problem has exact solution . In this experiment, the same set of optimal parameters is utilized. The values of approximate solution for different fractional orders and are displayed in Table 7. The comparison between absolute error and error estimate is shown in Figure 11(a), and Figure 11(b) shows the error function for fractional order .

**(a)**

**(b)**

#### 6. Conclusion

In this work, we developed a numerical scheme based on LT and inverse LT for approximation of the solution of fractional Riccati equations with ABC derivative. The inverse LT is approximated using the quadrature rule. The proposed method approximated the fractional Riccati equation with ABC derivative accurately and efficiently. From the results, it can be seen that this method is an excellent alternative for approximation of such type of equations.

#### Data Availability

Data are included within this paper.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

Xin Liu improved literature review and gave numerical applications of results and present examples. Kamran wrote the paper, and Yukun Yao gave error analysis and plotted error functions.

#### Acknowledgments

This research was partitively supported by the funds of HEC Pakistan.