A Collocation Method for Solving Fractional Riccati Differential Equation

We have introduced a Taylor collocation method, which is based on collocation method for solving fractional Riccati differential equation with delay term. This method is based on first taking the truncated Taylor expansions of the solution function in the fractional Riccati differential equation and then substituting their matrix forms into the equation. Using collocation points, we have the system of nonlinear algebraic equation. Then, we solve the system of nonlinear algebraic equation using Maple 13, and we have the coefficients of the truncated Taylor sum. In addition, illustrative examples are presented to demonstrate the effectiveness of the proposed method. Comparing the methodology with some known techniques shows that the present approach is relatively easy and highly accurate.


Introduction
The concept of fractional or noninteger order derivation and integration can be traced back to the genesis of integer order calculus itself [1,2]. The recent investigations in science and engineering have demonstrated that the dynamics of many systems may be described more accurately by using differential equations of noninteger order. The fractional differential equations FDEs have shown to be adequate models for various physical phenomena in areas like damping laws, diffusion processes, and so forth. For example, the nonlinear oscillation of earthquake can be modeled with fractional derivatives [3] and the fluid-dynamic traffic model with fractional derivatives [4], psychology [5] and so forth, [6][7][8][9].
In this paper, we present numerical and analytical solutions for the fractional Riccati differential equation with delay term * ( ) = ( ) + ( ) ( + ) + ( ) 2 ( ) , > 0, 0 < ≤ 1 (1) subject to the initial conditions where ( ), ( ), and ( ) are given functions, is a parameter describing the order of the fractional derivative and , are appropriate constants, and + > 0 for all ∈ [0, 1]. The general response expression contains a parameter describing the order of the fractional derivative that can be varied to obtain various responses. In the case of = 1, the fractional equation reduces to the classical Riccati differential equation. The importance of this equation usually arises in the optimal control problems [10]. The existing literature on fractional differential equations tends to focus on particular values for the order . In modern applications (see e.g., [11]) much more general values of the order appear in the equations, and therefore one needs to consider numerical and analytical methods to solve differential equations of arbitrary order. This equation is solved the numerically in [12][13][14].

Journal of Applied Mathematics
We seek the approximate solution of (1) under the conditions (2) with the fractional Taylor series as ( ) ∈ ( , ], where 0 < ≤ 1. In recently, collocation method has become a very useful technique for solving equations [15][16][17][18][19][20][21][22]. This method transforms each part of the equation into matrix form and using the collocation points as we get the nonlinear algebraic equation. Then this equation is solved, we obtained the coefficients, the approximate solutions for various . All computations are performed on the computer algebraic system Maple 13 in this paper.

Basic Definitions
In this section, we first give some basic definitions and then present properties of fractional calculus [1,2,23].

Fundamental Relations
In this section, we consider the fractional Ricatti differential equations (1). We use the Taylor matrix method [15][16][17][18][19][20][21][22] to find the truncated Taylor series expansions of each term in expression at = and their matrix representations for solving th order linear fractional part and nonlinear part. We first consider the solution ( ) of (1) defined by a truncated Taylor series (3). Then, we have the matrix form of the solution ( ) where ] .
Journal of Applied Mathematics 3 Then, the matrix representation of the function * ( ) becomes * where we compute the * X( ), where ] .
Then, so the matrix representation of fractional differential part as * Additionally, using (11) we can write where Moreover, since [21,22] where and using collocation points in (11) where then we construct the following relation Hence, the fundamental matrix relation of (1) is Finally, we obtained matrix representation of the condition in (2)

Method of Solution
Using collocation points in (4), we can write (24) or briefly the fundamental matrix equation where 4 Journal of Applied Mathematics where To obtain the solution of (1) under conditions (2) ] .
So, we obtained a system of ( + 1) nonlinear algebraic equations with unknown Taylor coefficients. We can easily check the accuracy of the method. Since the truncated fractional Taylor series (3)  (32)

Examples
In order to illustrate the effectiveness of the method proposed in this paper, several numerical examples are carried out in this section. In the following computations, for convenience, absolute errors between th-order approximate values and the corresponding exact values ex as = | − ex | are determined and all computations performed computer algebraic system with mathematical programing in Maple 13. Then, ( ) = Γ(3) 3/2 /Γ(5/2) − 2 3 − 2 , ( ) = − 2 , ( ) = 1. We assume that = 1/2, 0 ≤ ≤ 1, and we seek the approximate solutions by Taylor series, for = 0, = 4 with collocation points being Fundamental matrix relation of this problem is where  ] , Comparison of numerical results with the exact solution is shown in Table 1 for various .
Example 2. Let us consider the following fractional Riccati equation [13] * ( ) = 2 ( ) + 1, 0 < ≤ 1 subject to the initial condition Then, ( ) = 1, ( ) = 0, ( ) = 1. Fundamental matrix relation of this problem is Also, we have the matrix representation of conditions, The exact solution, when = 1, is We approximately solve the fractional Ricatti equation for = 12 and obtained the approximate solution for = 1,   Table 2, present method is in high agreement with the exact solution than homotopy perturbation method [13]. Moreover, using the numerical result in Table 2, Figures 1 and  2 are plotted. Figure 2 shows that numerical solution PM is so closed to the exact solution.

Conclusion
In real world systems, delays can be recognised everywhere and there has been widespread interest in the study of delay differential equations for many years. Although it seems natural to model certain processes and systems in engineering and other sciences with this kind of equation, only in the last few years has the attention of the scientific community been devoted to them.
In this study, we present a Taylor collocation method for the numerical solutions of fractional Riccati differential equation with delay term. This method transforms fractional Riccati differential equation with delay term into matrix equations. The desired approximate solutions can be determined by solving the resulting system, which can be effectively computed using symbolic computing codes on Maple 13. Examples show that Taylor collocation method has been successfully applied to find the approximate solutions of the fractional Riccati differential equation. Graphics and tables show that this method is extremely effective and practical for this sort of approximate solutions.