• Views 787
• Citations 16
• ePub 29
• PDF 582
`Abstract and Applied AnalysisVolume 2014, Article ID 626275, 8 pageshttp://dx.doi.org/10.1155/2014/626275`
Research Article

## New Ultraspherical Wavelets Spectral Solutions for Fractional Riccati Differential Equations

1Department of Mathematics, Faculty of Science, King Abdulaziz University, Saudi Arabia
2Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt

Received 19 March 2014; Accepted 29 April 2014; Published 22 May 2014

Copyright © 2014 W. M. Abd-Elhameed and Y. H. Youssri. 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.

#### Abstract

We introduce two new spectral wavelets algorithms for solving linear and nonlinear fractional-order Riccati differential equation. The suggested algorithms are basically based on employing the ultraspherical wavelets together with the tau and collocation spectral methods. The main idea for obtaining spectral numerical solutions depends on converting the differential equation with its initial condition into a system of linear or nonlinear algebraic equations in the unknown expansion coefficients. For the sake of illustrating the efficiency and the applicability of our algorithms, some numerical examples including comparisons with some algorithms in the literature are presented.

#### 1. Introduction

The fractional calculus deals with derivatives and integrals to an arbitrary order (real or complex order). Fractional ordinary differential equations are crucial in many fields. For example, many physical phenomena in acoustics, damping laws and electroanalytical chemistry, neuron modeling, diffusion processing, and material sciences (see, e.g., [13]) are described by differential equations of fractional order.

The use of fractional differential and integral operators in mathematical models has gained great importance in recent years. Several forms of fractional differential equations have been proposed in standard models, and there has been significant interest in developing numerical schemes for their solutions. There are numerous methods that deal with these types of equations; some of these methods are Laplace transform [4], Adomian decomposition method (ADM) [5], variational iteration method (VIM) [6], fractional difference method (FDM) [7], a quadrature tau method [8], and a shifted Jacobi spectral method [9].

Spectral methods are a class of techniques used extensively in various fields such as applied mathematics and scientific computing (see, e.g., [1012]). These methods aim to find numerical solutions to different kinds of differential equations. They are used in solving ordinary, partial, and fractional differential equations. When spectral methods are used, the numerical solution is assumed as an expansion in terms of certain special functions which are called “basis functions.” A great advantage of spectral methods is that they take on a global approach unlike finite element methods which use a local approach. There are three well-known versions of spectral methods, namely, tau, collocation, and Galerkin methods (see, e.g., [13]). The choice of the suitable spectral method suggested for solving the given equation depends certainly on the type of the differential equation and on the type of the boundary conditions governed by it.

The ultraspherical polynomials have received considerable attention in recent decades, from both theoretical and practical points of view (see, e.g., [14]). Some authors are interested in employing these polynomials for solving various kinds of differential equations. In this respect, Elgindy and Smith-Miles in [15] treated boundary value problems and integral and integrodifferential equations using ultraspherical integration matrices. Moreover, Doha and Abd-Elhameed employed ultraspherical polynomials for solving one- and two-dimensional second-order differential equations [16]. In addition, the same authors in [17] developed some accurate spectral solutions for the parabolic and elliptic partial differential equations based on the ultraspherical tau method.

The use of wavelets in solving ordinary and fractional differential equations has attracted a large number of authors. For example, in a huge number of articles, Legendre and Chebyshev wavelets are employed for treating ordinary differential equations as well as fractional differential equations (see, e.g., [1820]). To the best of our knowledge, no articles deal with the application of ultraspherical wavelets in handling fractional-order differential equations. This motivates our interest in such kind of wavelets hoping to employ them in various practical applications. Another motivation is that the Chebyshev and Legendre wavelets can be deduced as special cases of the ultraspherical wavelets.

The fractional-order Riccati differential equation is one of the important fractional-order differential equations. This equation has applications in various disciplines, such as random processes, optimal control, and diffusion problems. We refer the interested readers in the fundamental theories of this equation and its applications to the book of Reid [21].

The main objective in this paper is to introduce two efficient spectral algorithms based on using ultraspherical wavelets for solving fractional-order Riccati differential equation. The two algorithms, namely, ultraspherical wavelets tau method (UWTM) and ultraspherical wavelets collocation method (UWCM), are implemented and presented for handling both of linear and nonlinear fractional-order Riccati differential equations.

The paper is organized as follows. Section 2 is concerned with presenting some necessary definitions and mathematical preliminaries of the fractional calculus theory which are required for establishing our results. Also, in this section, some properties of ultraspherical polynomials and their shifted ones are presented. Section 3 is devoted to constructing the ultraspherical wavelets. In Section 4, some new spectral wavelets solutions for fractional-order Riccati differential equations are developed based on employing the ultraspherical wavelets which were constructed in Section 3. In Section 5, the convergence of the suggested ultraspherical wavelets expansion is carefully investigated. Some numerical examples and discussions are given in Section 6 aiming to illustrate the efficiency, simplicity, and applicability of the suggested algorithms. Some conclusions are given in Section 7.

#### 2. Preliminaries

##### 2.1. Basic Definitions of Fractional Calculus

This section is concerned with presenting some notations, definitions, and preliminary facts of the fractional calculus theory which will be useful throughout this paper.

Definition 1. The Riemann-Liouville fractional integral operator of order on the usual Lebesgue space is defined as The operator has the following properties: where , , and .

Definition 2. The Riemann-Liouville fractional derivative of order is defined by The following well-known properties for the operator are useful: where the ceiling notation denotes the smallest integer greater than or equal to and . Recall that, for , the fractional-differential operator coincides with the usual differential operator of integer-order. Similar to integer-order differentiation, also is a linear operator; that is, where and are constants.

Definition 3. The Caputo definition of fractional differential operator is given by where , .
The operator satisfies the following two basic properties for For more details on mathematical properties of fractional derivatives and integrals, see, for example, [4, 22].

##### 2.2. Some Properties of Ultraspherical Polynomials and Their Shifted Ones

The ultraspherical polynomials (a special type of Jacobi polynomials) associated with the real parameter are a sequence of orthogonal polynomials defined on (−1, 1), with respect to the weight function . The orthogonality relation is given by where

It should be noted here that the ultraspherical polynomials are normalized such that . This normalization is characterized by an advantage that the polynomials are identical with the Chebyshev polynomials of the first kind , are the Legendre polynomials , and is equal to , where are the Chebyshev polynomials of the second kind.

The polynomials may be generated by using the recurrence relation with the initial values and .

For more properties and relations of ultraspherical polynomials, see, for instance, [23].

The shifted ultraspherical polynomials are a sequence of orthogonal polynomials defined on , with respect to the weight function ; that is, They also may be generated by using the recurrence relation with the initial values and .

All relations and properties of ultraspherical polynomials can be easily transformed to give the corresponding relations and properties of the shifted ultraspherical polynomials.

Now, the following integral formula (see [23]) is needed: Also, the following theorem is essential in investigating the convergence analysis for the suggested ultraspherical wavelets expansion.

Theorem 4 (Bernstein-type inequality; see [24]). The following inequality holds for ultraspherical polynomials:

#### 3. Construction of Ultraspherical Wavelets

Wavelets constitute a family of functions constructed from dilation and translation of single function called the mother wavelet. When the dilation parameter and the translation parameter vary continuously, we have the following family of continuous wavelets: We define the ultraspherical wavelets as they have five arguments: can be assumed to be any positive integer, is the order for the ultraspherical polynomial, is the known ultraspherical parameter, and is the normalized time. Explicitly, they are defined on the interval as where , , and

Remark 5. It is worth noting here that is identical to Legendre wavelets [25, 26], is identical to first kind Chebyshev wavelets [27, 28], and is identical to second kind Chebyshev wavelets [29, 30].

Now, consider a function defined on and suppose that may be expanded in terms of the ultraspherical wavelets as where Also, we can assume that the function can be approximated in terms of ultraspherical wavelets as

#### 4. Spectral Solutions of Fractional-Order Riccati Differential Equation

In this section, with the aid of the ultraspherical wavelets expansion together with the two well-known spectral methods, namely, tau and collocation methods, we numerically solve the following fractional-order Riccati differential equation: subject to the initial condition The following lemma is needed for developing the wavelets solutions of fractional-order Riccati differential equation.

Lemma 6. The fractional derivative of order of the ultraspherical wavelets (16) is given by

Proof. The result can be immediately obtained by applying the operator defined in (3) on the ultraspherical wavelets defined in (16).

##### 4.1. Case I: Linear Fractional-Order Riccati Differential Equation

In this subsection, we consider (21) in case of ; that is, our aim is to discuss the following fractional-order Riccati linear differential equation: subject to the initial condition (22).

Now, assume that the approximate solution of (24) can be expanded in terms of the ultraspherical wavelets as and if we substitute (25) into (24), then the residual of this equation is given by The application of the typical tau method (see, e.g., [13]) implies that where Moreover, the use of the initial condition (22) yields Equations (27) and (29) generate a set consisting of equations in the unknown expansion coefficients, . This system of equations can be solved with the aid of any suitable solver such as Gauss-Seidel method (see, [31]). Hence, the unknown components of the vector can be obtained, and accordingly the required spectral wavelets solution given by (25) can be obtained.

##### 4.2. Case II: Nonlinear Fractional-Order Riccati Differential Equation

In this subsection, we aim to handle the following nonlinear fractional-order Riccati differential equation: subject to the initial condition (22).

Assume that the approximate solution of (30) can be expanded in terms of the ultraspherical wavelets as in (25). Substitution of (25) into (30) enables one to write the residual of (30) as The application of the collocation method (see, e.g., [13]) implies that where the points , are selected to be the first roots of . Equations (29) and (32) generate a set of equations in the unknown expansion coefficients . This system of equations can be solved with the aid of Newton’s iterative method (see [32]), to obtain the unknown components of the vector , and hence the required spectral wavelets solution given by (25) can be obtained.

#### 5. Convergence Analysis

In this section, we state and prove the following important theorem which ascertains that the ultraspherical wavelets expansion of a function with bounded second derivative converges uniformly to .

Theorem 7. A function , , , can be expanded as an infinite series of ultraspherical wavelets, which converges uniformly to , provided . Explicitly, the expansion coefficients in (20) satisfy the inequality

Proof. From relations (16) and (19), one can write and the integration of the right hand side of (34) by parts with the aid of relation (13) yields If we integrate (35) by parts again, and make use of the substitution: , then the coefficient takes the form Now, assuming that , taking into account the assumption , and with the aid of Theorem 4, we obtain Knowing that , , and with the aid of relation (17), we get This completes the proof of the theorem.

Note. It is to be noted here that, for large values of and , and making use of the well-known Stirling’s formula (see, [33]), it can be easily shown that is of .

#### 6. Numerical Results and Comparisons

In this section, we apply the presented algorithms in Section 4, for obtaining numerical solutions for linear and nonlinear fractional-order Riccati differential equations. In addition, we compare our spectral ultraspherical algorithms with some other numerical algorithms that exist in the literature aiming to demonstrate the efficiency and rapid convergence of the proposed algorithms.

Example 8. Consider the nonlinear Riccati differential equation (see, [3436]) The exact solution for (39), when , is given by In Table 1, we compare the numerical solution resulting from the application of UWCM, for the case that corresponds to , , with the generalized Euler method (GEM) and Laplace-Adomian-Padé method (LAPM) [34], while in Table 2, we compare the solution resulting from the application of UWCM, for the case that corresponds to , , , with the two solutions resulting from the application of homotopy perturbation method (HPM) developed in [35] and the variational iteration method (VIM) developed in [36]. In addition, Figure 1 illustrates different solutions for Example 8 in case of for various values of .

Table 1: Comparison between UWCM, GEM, and LAPM for Example 8 for .
Table 2: Comparison between UWCM, HPM, and VIM for Example 8 for .
Figure 1: Different solutions of Example 8.

Example 9. Consider the nonlinear Riccati differential equation (see, [36, 37]) The exact solution of (41), when , is . In Table 3, the maximum absolute error is listed for the case that corresponds to , and for various values of , while in Table 4, we give a comparison between the numerical solution resulting from the application of UWCM, for the case that corresponds to , , , with the two numerical solutions resulting from the application of the two methods, namely, variational iteration method (VIM) and second kind Chebyshev wavelets method (CUWM) developed, respectively, in [36, 37].

Table 3: Maximum absolute error for Example 9.
Table 4: Comparison between UWCM, VIM, and CUWM for Example 9 for .

Remark 10. It is worth noting here that the results displayed in Table 3 show that the best absolute error resulting from the application of UWCM does not exceed (), while the best error resulting from the application of the two methods, namely, VIM and CUWM developed in [36, 37], does not exceed (). Moreover, the results presented in Table 4 appear to agree closely where the UWCM deviates somewhat, because of the expected high accuracy of the UWCM if compared with the results obtained by VIM and CUWM [36, 37].

Remark 11. The results of Table 3 ascertain that the results that correspond to the first kind of Chebyshev wavelets expansion (in case of ) are not always better than the other ultraspherical wavelets expansions (see, [17, 38]).

Example 12. Consider the composite fractional oscillation equation (see, [36]) The exact solution for (42) is . We apply the UWTM for problem (42), for the case that corresponds to , . In detail, the approximate solution is given by where, , is as given in (17), .
In virtue of (27), the solution of the resulting linear system is and consequently which is the exact solution.

#### 7. Conclusions

This paper is concerned with presenting and implementing new algorithms for obtaining numerical spectral solutions for fractional-order Riccati differential equations. These methods are essentially based on employing the ultraspherical wavelets along with the application of the two spectral methods, namely, tau and collocation methods. The developed algorithms are applicable and easy in implementation and of high efficiency. The obtained numerical results indicate that the first kind of Chebyshev wavelets expansion is not always better than the other ultraspherical wavelets expansions.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. (130-269-D1435). The authors, therefore, acknowledge with thanks DSR technical and financial support.

#### References

1. Q. M. Al-Mdallal, M. I. Syam, and M. N. Anwar, “A collocation-shooting method for solving fractional boundary value problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 12, pp. 3814–3822, 2010.
2. Y. Çenesiz, Y. Keskin, and A. Kurnaz, “The solution of the Bagley-Torvik equation with the generalized Taylor collocation method,” Journal of the Franklin Institute: Engineering and Applied Mathematics, vol. 347, no. 2, pp. 452–466, 2010.
3. V. Daftardar-Gejji and H. Jafari, “Adomian decomposition: a tool for solving a system of fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 301, no. 2, pp. 508–518, 2005.
4. 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 of Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1999.
5. H. Jafari and S. Seifi, “Solving a system of nonlinear fractional partial differential equations using homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 1962–1969, 2009.
6. S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009.
7. M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for two-sided space-fractional partial differential equations,” Applied Numerical Mathematics: An IMACS Journal, vol. 56, no. 1, pp. 80–90, 2006.
8. A. H. Bhrawy, A. S. Alofi, and S. S. Ezz-Eldien, “A quadrature tau method for fractional differential equations with variable coefficients,” Applied Mathematics Letters, vol. 24, no. 12, pp. 2146–2152, 2011.
9. E. H. Doha, A. H. Bhrawy, D. Baleanu, and S. S. Ezz-Eldien, “On shifted Jacobi spectral approximations for solving fractional differential equations,” Applied Mathematics and Computation, vol. 219, no. 15, pp. 8042–8056, 2013.
10. W. M. Abd-Elhameed, E. H. Doha, and Y. H. Youssri, “Efficient spectral-Petrov-Galerkin methods for third- and fifth-order differential equations using general parameters generalized Jacobi polynomials,” Quaestiones Mathematicae, vol. 36, no. 1, pp. 15–38, 2013.
11. E. H. Doha, W. M. Abd-Elhameed, and Y. H. Youssri, “Second kind Chebyshev operational matrix algorithm for solving differential equations of Lane-Emden type,” New Astronomy, vol. 23-24, pp. 113–117, 2013.
12. E. H. Doha, W. M. Abd-Elhameed, and M. A. Bassuony, “New algorithms for solving high even-order differential equations using third and fourth Chebyshev-Galerkin methods,” Journal of Computational Physics, vol. 236, pp. 563–579, 2013.
13. C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Series in Computational Physics, Springer, New York, NY, USA, 1988.
14. K. T. Elgindy and K. A. Smith-Miles, “On the optimization of Gegenbauer operational matrix of integration,” Advances in Computational Mathematics, vol. 39, no. 3-4, pp. 511–524, 2013.
15. K. T. Elgindy and K. A. Smith-Miles, “Solving boundary value problems, integral, and integro-differential equations using Gegenbauer integration matrices,” Journal of Computational and Applied Mathematics, vol. 237, no. 1, pp. 307–325, 2013.
16. E. H. Doha and W. M. Abd-Elhameed, “Efficient spectral-Galerkin algorithms for direct solution of second-order equations using ultraspherical polynomials,” SIAM Journal on Scientific Computing, vol. 24, no. 2, pp. 548–571, 2002.
17. E. H. Doha and W. M. Abd-Elhameed, “Accurate spectral solutions for the parabolic and elliptic partial differential equations by the ultraspherical tau method,” Journal of Computational and Applied Mathematics, vol. 181, no. 1, pp. 24–45, 2005.
18. W. M. Abd-Elhameed, E. H. Doha, and Y. H. Youssri, “New wavelets collocation method for solving second-order multipoint boundary value problems using Chebyshev polynomials of third and fourth kinds,” Abstract and Applied Analysis, vol. 2013, Article ID 542839, 9 pages, 2013.
19. L. Zhu and Q. Fan, “Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 6, pp. 2333–2341, 2012.
20. I. Sadek, T. Abualrub, and M. Abukhaled, “A computational method for solving optimal control of a system of parallel beams using Legendre wavelets,” Mathematical and Computer Modelling, vol. 45, no. 9-10, pp. 1253–1264, 2007.
21. W. T. Reid, Riccati Differential Equations, vol. 86 of Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1972.
22. K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
23. G. E. Andrews, R. Askey, and R. Roy, Special Functions, vol. 71 of Encyclopedia of Mathematics and its Applications, Cambridge University Press, Cambridge, UK, 1999.
24. C. Giordano and A. Laforgia, “On the Bernstein-type inequalities for ultraspherical polynomials,” Journal of Computational and Applied Mathematics, vol. 153, pp. 243–248, 2003.
25. M. Razzaghi and S. Yousefi, “Legendre wavelets direct method for variational problems,” Mathematics and Computers in Simulation, vol. 53, no. 3, pp. 185–192, 2000.
26. S. A. Yousefi, “Legendre wavelets method for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1417–1422, 2006.
27. E. Babolian and F. Fattahzadeh, “Numerical solution of differential equations by using Chebyshev wavelet operational matrix of integration,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 417–426, 2007.
28. Y. Li, “Solving a nonlinear fractional differential equation using Chebyshev wavelets,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2284–2292, 2010.
29. K. Maleknejad, S. Sohrabi, and Y. Rostami, “Numerical solution of nonlinear Volterra integral equations of the second kind by using Chebyshev polynomials,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 123–128, 2007.
30. W. M. Abd-Elhameed, E. H. Doha, and Y. H. Youssri, “New spectral second kind Chebyshev wavelets algorithm for solving linear and nonlinear second-order differential equations involving singular and Bratu type equations,” Abstract and Applied Analysis, vol. 2013, Article ID 715756, 9 pages, 2013.
31. R. Barrett, M. Berry, T. F. Chan et al., Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods, Society for Industrial and Applied Mathematics (SIAM), 1994.
32. R. S. Dembo, S. C. Eisenstat, and T. Steihaug, “Inexact Newton methods,” SIAM Journal on Numerical Analysis, vol. 19, no. 2, pp. 400–408, 1982.
33. Y.-C. Li, “A note on an identity of the gamma function and Stirling's formula,” Real Analysis Exchange, vol. 32, no. 1, pp. 267–272, 2007.
34. N. A. Khan, A. Ara, and N. A. Khan, “Fractional-order Riccati differential equation: analytical approximation and numerical results,” Advances in Difference Equations, vol. 2013, article 185, 2013.
35. 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.
36. F. Yin, J. Song, J. H. Leng, and F. Lu, “Couple of the variational iteration method and fractional-order legendre functions method for fractional differential equations,” The Scientific World Journal, vol. 2014, Article ID 928765, 9 pages, 2014.
37. Y. Wang and Q. Fan, “The second kind Chebyshev wavelet method for solving fractional differential equations,” Applied Mathematics and Computation, vol. 218, no. 17, pp. 8592–8601, 2012.
38. W. A. Light, “Are Chebysheve expansions really better,” Bulletin of the Institute of Mathematics and Its Applications, vol. 22, pp. 180–181, 1986.