Research Article  Open Access
Xin Liu, Kamran, Yukun Yao, "Numerical Approximation of Riccati Fractional Differential Equation in the Sense of CaputoType Fractional Derivative", Journal of Mathematics, vol. 2020, Article ID 1274251, 12 pages, 2020. https://doi.org/10.1155/2020/1274251
Numerical Approximation of Riccati Fractional Differential Equation in the Sense of CaputoType Fractional Derivative
Abstract
The Riccati differential equation is a wellknown 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 realworld 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 Laplacetransform 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 fractionalorder RDE have been studied by a large number of researchers. For example, the authors [20] studied fractionalorder RDE with ABC derivatives, and they established the existence and uniqueness results using Banach fixedpoint 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 fractionalorder RDE using the modified homotopy perturbation method. Numerical solution of fractionalorder 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 fractionalorder 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 fractionalorder 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 fractionalorder 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 fractionalorder RDE can be found in [33–42] and references therein. In this work, we approximate the solution of fractionalorder 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 twoparameter 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 firstorder 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 twoparameter 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 xintercept 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 realworld problem. So we need numerical methods to approximate the solutions of the problems. In this section, we consider fractionalorder 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.
References
 U. N. Katugampola, “A new approach to generalized fractional derivatives,” Bulletin of Mathematical Analysis and Applications, vol. 6, no. 4, pp. 1–15, 2014. View at: Google Scholar
 D. Baleanu, A. Jajarmi, S. S. Sajjadi, and D. Mozyrska, “A new fractional model and optimal control of a tumorimmune surveillance with nonsingular derivative operator,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 29, no. 8, Article ID 083127, 2019. View at: Publisher Site  Google Scholar
 A. Atangana and J. F. GómezAguilar, “Decolonisation of fractional calculus rules: breaking commutativity and associativity to capture more natural phenomena,” The European Physical Journal Plus, vol. 133, no. 4, p. 166, 2018. View at: Publisher Site  Google Scholar
 I. Podlubny, Fractional Differential Equations: an Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, vol. 198, Elsevier, Amsterdam, Netherlands, 1998.
 A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204, Elsevier Science Limited, Amsterdam, Netherlands, 2006.
 M. Caputo and M. Fabrizio, “A new definition of fractional derivative without singular kernel,” Progress in Fractional Differentiation and Applications, vol. 1, no. 2, pp. 1–13, 2015. View at: Google Scholar
 A. Atangana and D. Baleanu, “New fractional derivatives with nonlocal and nonsingular kernel: theory and application to heat transfer model,” Thermal Science, vol. 20, no. 2, pp. 763–769, 2016. View at: Publisher Site  Google Scholar
 J. Hristov, “On the AtanganaBaleanu derivative and its relation to the fading memory concept: the diffusion equation formulation,” in Fractional Derivatives with MittagLeffler Kernel. Studies in Systems, Decision and Control, Springer, Cham, Switzerland, 2019. View at: Publisher Site  Google Scholar
 W. T. Reid, Riccati Differential Equations, Elsevier, Amsterdam, Netherlands, 1972.
 S. Momani and N. Shawagfeh, “Decomposition method for solving fractional Riccati differential equations,” Applied Mathematics and Computation, vol. 182, no. 2, pp. 1083–1092, 2006. View at: Publisher Site  Google Scholar
 Y. Pala and M. O. Ertas, “An analytical method for solving general riccati equation,” International Journal of Mathematical and Computational Sciences, vol. 11, no. 3, pp. 125–130, 2017. View at: Google Scholar
 J. Biazar and M. Eslami, “Differential transform method for quadratic riccati differential equation,” International Journal of Nonlinear Science, vol. 9, no. 4, pp. 444–447, 2010. View at: Google Scholar
 P.Y. Tsai and C. o.K. Chen, “An approximate analytic solution of the nonlinear riccati differential equation,” Journal of the Franklin Institute, vol. 347, no. 10, pp. 1850–1862, 2010. View at: Publisher Site  Google Scholar
 M. A. Z. Raja, I. M. Qureshi, and J. A. Khan, “A new stochastic approach for solution of riccati differential equation of fractional order,” Annals of Mathematics and Artificial Intelligence, vol. 60, no. 34, pp. 229–250, 2010. View at: Publisher Site  Google Scholar
 H. Aminikhah, A. H. R. Sheikhani, and H. Rezazadeh, “Approximate analytical solutions of distributed order fractional riccati differential equation,” Ain Shams Engineering Journal, vol. 9, no. 4, pp. 581–588, 2018. View at: Publisher Site  Google Scholar
 N. H. Sweilam, M. M. Khader, and R. F. AlBar, “Numerical studies for a multiorder fractional differential equation,” Physics Letters A, vol. 371, no. 12, pp. 26–33, 2007. View at: Publisher Site  Google Scholar
 M. A. Snyder, Chebyshev Methods in Numerical Approximation, vol. 2, PrenticeHall, Upper Saddle River, NJ, USA, 1966.
 M. G. Sakar, A. Akgül, and D. Baleanu, “On solutions of fractional Riccati differential equations,” Advances in Difference Equations, vol. 2017, p. 39, 2017. View at: Publisher Site  Google Scholar
 O. A. Arqub and B. Maayah, “Modulation of reproducing kernel Hilbert space method for numerical solutions of Riccati and Bernoulli equations in the AtanganaBaleanu fractional sense,” Chaos, Solitons & Fractals, vol. 125, pp. 163–170, 2019. View at: Google Scholar
 M. I. Syam and M. AlRefai, “Fractional differential equations with AtanganaBaleanu fractional derivative: analysis and applications,” Chaos, Solitons & Fractals: X, vol. 2, Article ID 100013, 2019. View at: Publisher Site  Google Scholar
 M. M. Khader, “Introducing an efficient modification of the homotopy perturbation method by using Chebyshev polynomials,” Arab Journal of Mathematical Sciences, vol. 18, no. 1, pp. 61–71, 2012. View at: Publisher Site  Google Scholar
 Z. Odibat and S. Momani, “Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order,” Chaos, Solitons & Fractals, vol. 36, no. 1, pp. 167–174, 2008. View at: Publisher Site  Google Scholar
 Ş. Yüzbaşı, “Numerical solutions of fractional Riccati type differential equations by means of the Bernstein polynomials,” Applied Mathematics and Computation, vol. 219, no. 11, pp. 6328–6343, 2013. View at: Google Scholar
 M. M. Khader, “Numerical treatment for solving fractional Riccati differential equation,” Journal of the Egyptian Mathematical Society, vol. 21, no. 1, pp. 32–37, 2013. View at: Publisher Site  Google Scholar
 H. Jafari and H. Tajadodi, “He’s variational iteration method for solving fractional Riccati differential equation,” International Journal of Differential Equations, vol. 2010, Article ID 764738, 8 pages, 2010. View at: Publisher Site  Google Scholar
 X. Y. Li, B. Y. Wu, and R. T. Wang, “Reproducing kernel method for fractional Riccati differential equations,” Abstract and Applied Analysis, vol. 2014, Article ID 970967, 6 pages, 2014. View at: Publisher Site  Google Scholar
 M. Izadi, “Fractional polynomial approximations to the solution of fractional Riccati equation,” Punjab University Journal of Mathematics, vol. 51, no. 11, pp. 123–141, 2019. View at: Google Scholar
 N. H. Sweilam, M. M. Khader, and A. M. S. Mahdy, “Numerical studies for solving fractional Riccati differential equation,” Applications and Applied Mathematics, vol. 7, no. 2, pp. 595–608, 2012. View at: Google Scholar
 M. I. Syam, H. I. Siyyam, and I. AlSubaihi, “TauPath following method for solving the Riccati equation with fractional order,” Journal of Computational Methods in Physics, vol. 2014, Article ID 207916, 7 pages, 2014. View at: Publisher Site  Google Scholar
 H. Jafari, H. Tajadodi, and D. Baleanu, “A modified variational iteration method for solving fractional Riccati differential equation by Adomian polynomials,” Fractional Calculus and Applied Analysis, vol. 16, no. 1, pp. 109–122, 2013. View at: Publisher Site  Google Scholar
 M. M. Khashan, R. Amin, and M. I. Syam, “A new algorithm for fractional riccati type differential equations by using haar wavelet,” Mathematics, vol. 7, no. 6, p. 545, 2019. View at: Publisher Site  Google Scholar
 B. S. Kashkari and M. I. Syam, “A numerical approach for investigating a special class of fractional riccati equation,” Results in Physics, vol. 17, Article ID 103080, 2020. View at: Publisher Site  Google Scholar
 W. M. AbdElhameed and Y. H. Youssri, “New ultraspherical wavelets spectral solutions for fractional Riccati differential equations,” Abstract and Applied Analysis, vol. 2014, Article ID 626275, 8 pages, 2014. View at: Publisher Site  Google Scholar
 M. Merdan, “On the solutions fractional Riccati differential equation with modified RiemannLiouville derivative,” International Journal of Differential Equations, vol. 2012, Article ID 346089, 17 pages, 2012. View at: Publisher Site  Google Scholar
 H. Yaslan, “Numerical solution of fractional Riccati differential equation via shifted Chebyshev polynomials of the third kind,” Journal of Engineering Technology and Applied Sciences, vol. 2, no. 1, pp. 1–11, 2017. View at: Publisher Site  Google Scholar
 B. Lu, “Bäcklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations,” Physics Letters A, vol. 376, no. 2829, pp. 2045–2048, 2012. View at: Google Scholar
 S. W. Jeng and A. Kilicman, “Fractional riccati equation and its applications to rough heston model using numerical methods,” Symmetry, vol. 12, no. 6, p. 959, 2020. View at: Publisher Site  Google Scholar
 H. Z. Mjthap and S. N. AlAzzawi, “Mixing sumudu transform and adomain decomposition method for solving riccati equation of variable fractional order,” Journal of Interdisciplinary Mathematics, vol. 22, no. 8, pp. 1559–1563, 2019. View at: Publisher Site  Google Scholar
 A. Turan Dincel, “Solution to fractionalorder riccati differential equations using euler wavelet method,” Scientia Iranica, vol. 26, no. 3, pp. 1608–1616, 2019. View at: Google Scholar
 S. N. TuralPolat, “Thirdkind Chebyshev wavelet method for the solution of fractional order riccati differential equations,” Journal of Circuits, Systems and Computers, vol. 28, no. 14, Article ID 1950247, 2019. View at: Publisher Site  Google Scholar
 W. M. AbdElhameed and Y. H. Youssri, “Explicit shifted secondkind Chebyshev spectral treatment for fractional riccati differential equation,” Computer Modeling in Engineering & Sciences, vol. 121, no. 3, pp. 1029–1049, 2019. View at: Publisher Site  Google Scholar
 J. Biazar, L. B. Hasani, and Z. Ayati, “Asymptotic decomposition method for fractional order riccati differential equation,” Computational Methods for Differential Equations, 2020. View at: Google Scholar
 S. Mehmood, G. Farid, and G. Farid, “Fractional integrals inequalities for exponentially \(m\)convex functions,” Open Journal of Mathematical Sciences, vol. 4, no. 1, pp. 78–85, 2020. View at: Publisher Site  Google Scholar
 S. Mehmood, G. Farid, K. A. Khan, and M. Yussouf, “New hadamard and fejér–hadamard fractional inequalities for exponentially mconvex function,” Engineering and Applied Science Letter, vol. 3, pp. 45–55, 2020. View at: Publisher Site  Google Scholar
 S. Mehmood, G. Farid, K. A. Khan, and M. Yussouf, “New fractional hadamard and fejér–hadamard inequalities associated with exponentially (h, m)convex functions,” Engineering and Applied Science Letter, vol. 3, pp. 9–18, 2020. View at: Google Scholar
 S. I. Butt, M. Nadeem, and G. Farid, “On caputo fractional derivatives via exponential (s, m)convex functions,” Engineering and Applied Science Letter, vol. 3, no. 2, pp. 32–39, 2020. View at: Google Scholar
 Y. C. Kwun, G. Farid, W. Nazeer, S. Ullah, and S. M. Kang, “Generalized riemannliouville kfractional integrals associated with ostrowski type inequalities and error bounds of hadamard inequalities,” IEEE Access, vol. 6, pp. 64946–64953, 2018. View at: Publisher Site  Google Scholar
 Y. C. Kwun, M. S. Saleem, M. Ghafoor, W. Nazeer, and S. M. Kang, “Hermite–hadamardtype inequalities for functions whose derivatives are ηconvex via fractional integrals,” Journal of Inequalities and Applications, vol. 2019, p. 44, 2019. View at: Publisher Site  Google Scholar
 S. Kang, G. Abbas, G. Farid, and W. Nazeer, “A generalized fejér–hadamard inequality for harmonically convex functions via generalized fractional integral operator and related results,” Mathematics, vol. 6, no. 7, p. 122, 2018. View at: Publisher Site  Google Scholar
 A. Singh and A. Prakash, “Risk evaluation in information systems using continuous and discrete distribution laws,” Engineering and Applied Science Letters, vol. 3, no. 1, pp. 35–44, 2020. View at: Publisher Site  Google Scholar
 C. M. D. Simarmata, N. Susyanto, I. J. Hammadi, and C. Rahmaditya, “A mathematical model of smoking behaviour in Indonesia with densitydependent death rate,” Open Journal of Mathematical Sciences, vol. 4, no. 1, pp. 118–125, 2020. View at: Google Scholar
 G. M. MittagLeffler, “Sur la nouvelle fonction E(x),” Comptes Rendus Mathématique, vol. 137, no. 2, pp. 554–558, 1903. View at: Google Scholar
 P. Humbert and R. P. Agarwal, “Sur la fonction de MittagLeffler et quelquesunes de ses généralisations,” Bulletin des Sciences Mathématiques, vol. 77, no. 2, pp. 180–185, 1953. View at: Google Scholar
 W. McLean and V. Thomée, “Numerical solution via Laplace transforms of a fractional order evolution equation,” Journal of Integral Equations and Applications, vol. 22, no. 1, pp. 57–94, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Xin Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.