Legendre Wavelet Operational Matrix Method for Solution of Riccati Differential Equation

ALegendrewavelet operationalmatrixmethod (LWM) is presented for the solution of nonlinear fractional-orderRiccati differential equations, having variety of applications in quantum chemistry and quantum mechanics. The fractional-order Riccati differential equations converted into a system of algebraic equations using Legendre wavelet operational matrix. Solutions given by the proposed scheme are more accurate and reliable and they are compared with recently developed numerical, analytical, and stochastic approaches. Comparison shows that the proposed LWM approach has a greater performance and less computational effort for getting accurate solutions. Further existence and uniqueness of the proposed problem are given and moreover the condition of convergence is verified.


Introduction
In recent years, use of fractional-order derivative goes very strongly in engineering and life sciences and also in other areas of science.One of the best advantages of use of fractional differential equation is modeling and control of many dynamic systems.Fractional-order derivatives are used in fruitful way to model many remarkable developments in those areas of science such as quantum chemistry, quantum mechanics, damping laws, rheology, and diffusion processes [1][2][3][4][5] described through the models of fractional differential equations (FDEs).Modeling of a physical phenomenon depends on two parameters such as the time instant and the prior time history; because of this reason, reasonable modeling through fractional calculus was successfully achieved.The abovementioned advantages and applications of FDEs attracted researchers to develop efficient methods to solve FDEs in order to get accurate solutions to such problems and more active research is still going on in those areas.Most of the FDEs are complicated in their structure; hence finding exact solutions for them cannot be simple.Therefore, one can approach the best accurate solution of FDEs through analytical and numerical methods.Designing accurate or best solution to FDEs, many methods are developed in recent years; each method has its own advantages and limitations.This paper aims to solve a FDE called fractional-order Riccati differential equation, one of the important equations in the family of FDEs.The Riccati equations play an important role in engineering and applied science [6], especially in quantum mechanics [7] and quantum chemistry [8,9].Therefore, solutions to the Riccati differential equations are important to scientists and engineers.Solving fractional-order Riccati differential equation, the most significant methods are Adomian decomposition method [10], homotopy perturbation method [11][12][13][14], homotopy analysis method [15,16], Taylor matrix method [17], Haar wavelet method [18], and combination of Laplace, Adomian decomposition, and Padé approximation [19] methods.
Several numerical methods for approximating the solution of nonlinear fractional-order Riccati differential equations are known.Raja et al. [20] developed a stochastic technique based on particle swarm optimization and simulated annealing.They were used as a tool for rapid global search method and simulated annealing for efficient local search method.A fractional variational iteration method described in the Riemann-Liouville derivative has been applied in [21], to give an analytical approximate solution to nonlinear fractional Riccati differential equation.A combination of finite 2 International Journal of Mathematics and Mathematical Sciences difference method and Padé-variational iteration numerical scheme was proposed by Sweilam et al. [22].Moreover, an analytical scheme comprising the Laplace transform, the Adomian decomposition method (ADM), and the Padé approximation is given in [19].
However, the abovementioned methods have some restrictions and disadvantages in their performance.For example, very complicated and toughest Adomian polynomials are constructed in the Adomian decomposition method.In the variational iteration method identification Lagrange multiplier yields an underlying accuracy.The homotopy perturbation method needs a linear functional equation in each iteration to solve nonlinear equations; forming these functional equations is very difficult.The performance of the homotopy analysis method very much depends on the choice of the auxiliary parameter ℎ of the zero-order deformation equation.Moreover, the convergence region and implementation of these results are very small.
In recent years, wavelets theory is one of the growing and predominantly new methods in the area of mathematical and engineering research.It has been applied in vast range of engineering sciences; particularly, they are used very successfully for waveform representation and segmentations in signal analysis and time-frequency analysis and in the mathematical sciences it is used in thriving manner for solving variety of linear and nonlinear differential and partial differential equations and fast algorithms for easy implementation [23].Moreover, wavelets build a connection with fast numerical algorithms [24]; this is due to the fact that wavelets admit the exact representation of a variety of functions and operators.The application of Legendre wavelet and its operational matrix for solving differential, integral, and fractionalorder differential equations is thoroughly considered in [25,26].
In this work, the nonlinear Riccati differential equations of fractional-order are approached analytically by using Legendre wavelets method.The operational matrix of Legendre wavelet is generalized for fractional calculus in order to solve fractional and classical Riccati differential equations.The Legendre wavelet method (LWM) is illustrated by application, and obtained results are compared with recently proposed method for the fractional-order Riccati differential equation.We have adopted Legendre wavelet method to solve Riccati differential equations not only due to its emerging application but also due to its greater convergence region.
The rest of the paper is as follows.In Section 2, definitions and mathematical preliminaries of fractional calculus are presented.In Section 3, Legendre wavelet, its properties, function approximations, and generalized Legendre wavelet operational matrix fractional calculus are discussed.Section 4 establishes application of proposed method in the solution of Riccati differential equations, existence and uniqueness solution of the proposed problem, and convergence analyses of the proposed approach.Section 5 deals with the illustrative examples and their solutions by the proposed approach.Section 6 ends with our conclusion.

Preliminaries and Notations
The notations, definitions, and preliminary facts present in this section will be used in forthcoming sections of this work.Several definitions of fractional integrals and derivatives have been proposed after the logical definition given by Liouville.Important and few of these definitions include the Riemann-Liouville, the Caputo, the Weyl, the Hadamard, the Marchaud, the Riesz, the Grunwald-Letnikov, and the Erdelyi-Kober.As stated in [25], the Caputo fractional derivative uses initial and boundary conditions of integer order derivatives having some physical interpretations.Because of this specific reason, in this work, we will use the Caputo fractional derivative   proposed by Caputo [27] in the theory of viscoelasticity.
The Caputo fractional derivative of order  > 0( ∈ , −1 <  ≤ ,  ∈ ) and ℎ : (0, ∞) →  is continuous and is defined by where is the Riemann-Liouville fractional integral operator of order  > 0 and Γ is the gamma function.
The fractional integral of   ,  > −1 is given as Properties of fractional integrals and derivatives are as follows [28], for ,  > 0.
The fractional-order integral satisfies the semigroup property The integer order derivative   and fractional-order derivative   commute with each other: The fractional integral operator and fractional derivative operator do not satisfy the commutative property.In general, But in the reverse way we have International Journal of Mathematics and Mathematical Sciences 3

Generalized Legendre Wavelet Operational Matrix to Fractional Integration
3.1.Legendre Wavelets.A family of functions were constituted by wavelets and constructed from dilation and translation of a single function called mother wavelet.When the parameters  of dilation and  of translation vary continuously, following are the family of continuous wavelets [29]: If the parameters  and  are restricted to discrete values as  =  − 0 ,  =  0  − 0 ,  0 > 1,  0 > 0, and  and  are positive integers, following are the family of discrete wavelets: where  , () form a wavelet basis for  2 ().In particular, when  0 = 2 and  0 = 1,  , () form an orthonormal basis [29].

Operational Matrix of the Fractional Integration.
The integration of the Ψ() defined in ( 14) can be approximated by Legendre wavelet series with Legendre wavelet coefficient matrix : where the m-square matrix  is called Legendre wavelet operational matrix of integration.The -set of block-pulse functions is defined on [0, l) as follows: where  = 0, 1, 2, 3, . . ..

International Journal of Mathematics and Mathematical Sciences
The functions   are disjoint and orthogonal.That is, The orthogonality property of block-pulse function is obtained from the disjointness property.
An arbitrary function  ∈  2 [0, 1) can be expanded into block-pulse functions as where   are the coefficients of the block-pulse function, given by The Legendre wavelets can be expanded into -set of blockpulse functions as where The fractional integral of block-pulse function vector can be written as where   × is given in [32].Now, we introduce the derivation process of the Legendre wavelet operational matrix of the fractional integration: where the m-square matrix   m× m is called Legendre wavelet operational matrix of the fractional integration.
Using ( 23) and ( 25), we have From ( 25) and ( 26), we get and by ( 23), (27) becomes Then, the Legendre wavelet operational matrix   m× m of fractional integration is given by Following is the Legendre wavelet operational matrix   m× m of fractional-order integration, for the particular values of  = 2,  = 3, and  = 0.5:

Application to Fractional Riccati Differential Equation
In this section, we use the generalized Legendre wavelet operational matrix to solve nonlinear Riccati differential equation and we discuss the existence and uniqueness of solutions with initial conditions and convergence criteria of the proposed LWM approach.Consider the fractional-order Riccati differential equation of the form subject to the initial condition Let us suppose that the functions   (), (), (), and () are approximated using Legendre wavelet as follows: where , , , , and Ψ() are given in (14).Using (6), we can write By ( 25) and ( 32), (34) leads to where (0) =  ≈   0 Ψ(),  = (    m× m +   0 )  .Substituting (33) and ( 35) into (31), we have (36) Substituting ( 23) into (36), we have where , , , and  m× m are known.Equation (37) represents a system of nonlinear equations with unknown vector .This system of nonlinear equations can be solved by Newton method for the unknown vector  and we can get the approximation solution by including  into (35).

International Journal of Mathematics and Mathematical Sciences
Now, again from (42), (43), and (46), we get which implies that the integral equation ( 46) is equivalent to the initial value problem (32) and the theorem is proved.

Numerical Examples
In order to show the effectiveness of the Legendre wavelets method (LWM), we implement LWM to the nonlinear fractional Riccati differential equations.All the numerical experiments were carried out on a personal computer with some MATLAB codes.The specifications of PC are Intel core i5 processor and with Turbo boost up to 3.1 GHz and 4 GB of DDR3 memory.The following problems of nonlinear Riccati differential equations are solved with real coefficients.
Example 1.Consider the following nonlinear fractional Riccati differential equation: with initial condition Exact solution for  = 1 was found to be The integral representation of ( 54) and ( 55) is given by Let and then By substituting (59) and ( 60) into (58), we get the following system of algebraic equations: By solving the above system of linear equations, we can find the vector .Numerical results are obtained for different values of , , and .Solution obtained by the proposed LWM approach for  = 1,  = 1, and  = 3 is given in Figure 1 and for different values of  = 0.6, 0.7, 0.8, and 0.9 and for  = 2 and  = 5 is graphically given in Figure 2.
It can be seen from Figure 1 that the solution obtained by the proposed LWM approach is more close to the exact solution.Table 1 describes the efficiency of the proposed method by comparing with the methods in [20,22] through their absolute error.The following is used for the errors of the approximation ŷ() of (); that is, ‖− ŷ‖ = max |()− ŷ()|.
Table 1 shows that very high accuracies are obtained for = 3 and = 5 by the present method.Example 2. Consider another fractional-order Riccati differential equation: with initial condition Exact solution for the above equation was found to be The integral representation of (62) and ( 63) is given by Let and then By substituting (66) and (67) in (62), we get the following system of algebraic equations: By solving the above system of linear equations, we can find the vector .Numerical results are obtained for different values of , , and .Results obtained by LWM for  = 1,  = 2, and  = 3 are shown in Figure 3 and it can be seen from the figure that solution given by the LWM merely coincides with the exact solution.Figure 4 shows the obtained results of ( 62) and ( 63) by LWM for different values of  and for  = 2 and International Journal of Mathematics and Mathematical Sciences  = 5.Table 2 describes the efficiency of the proposed method by comparing with the methods in [20,22] through their absolute error.Table 1 shows that very high accuracies are obtained for  = 3 and  = 5 by the present method and from these results we can identify that guarantee of convergence of the proposed LWM approach is very high.

Example 3. Let us consider another problem of nonlinear
Riccati differential equation: with initial condition When  = 1, its exact solution is given by where   () is the Bessel function of first kind.The integral representation of (69) and ( 70) is given by Let and then By substituting (73) and ( 74) into (69), we get the following system of algebraic equations: Ψ () = 1 + 2 Γ ( + 1) + 2Γ () (  Ψ ())
(75) By solving the above system of linear equations, we can find the vector .Numerical results are obtained for different values of , , and .Obtained results for (69) and (70) are shown in Figures 5 and 6. Figure 5 shows the solutions obtained by LWM for different values of  and for  = 2 and  = 4. Figure 6 compares the solution obtained by LWM with the exact solution of (69) and (70) when  = 1,  = 1, and  = 2.So far, there are no published results of absolute error for this problem and hence we are unable to compare absolute error of our method with the existing methods.From these results, we can see that the proposed LWM approach gives the solution which is very close to the exact solution and outperformed recently developed approaches for the nonlinear fractional-order Riccati differential equations in terms of solution quality and convergence criteria.

Conclusions
Nonlinear fractional-order Riccati differential equations play an important role in the modeling of many biological, physical, chemical, and real life problems.Therefore, it is necessary to develop a method which would give more accurate solutions to such type of problems with greater convergence criteria.In this work, a Legendre's wavelet operational matrix method called LWM was proposed for solving nonlinear fractional-order Riccati differential equations.Comparison was made for the solutions obtained by the proposed method and with the other recent approaches developed for the same problem through their error analysis; obtained results show that the proposed LWM yields more accurate and reliable solutions even for small values of  and , which assures the best approximate solution in less computational effort.Further, we have discussed the convergence criteria of proposed scheme, which indeed provides the guarantee of consistency and stability of the proposed LWM scheme for the solutions of nonlinear fractional Riccati differential equations.

Figure 2 :
Figure 2: Numerical results of Example 1 by LWM for different values of .

)Figure 4 :
Figure 4: Numerical results of Example 2 by LWM for different values of .

Figure 5 :Figure 6 :
Figure 5: Numerical results of Example 3 by LWM for different values of .