Advances in Mathematical Physics

Volume 2013 (2013), Article ID 954015, 9 pages

http://dx.doi.org/10.1155/2013/954015

## Approximate Analytical Solution for Nonlinear System of Fractional Differential Equations by BPs Operational Matrices

^{1}Department of Mathematics, Imam Khomeini International University, P.O. Box 34149-16818, Qazvin, Iran^{2}Department of Mathematics and Computer Sciences, Cankaya University, 06530 Ankara, Turkey^{3}Institute of Space Sciences, P.O. Box MG-23, 077125 Magurele-Bucharest, Romania^{4}Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, P.O. Box 80204, Jeddah 21589, Saudi Arabia

Received 20 February 2013; Revised 6 March 2013; Accepted 7 March 2013

Academic Editor: José Tenreiro Machado

Copyright © 2013 Mohsen Alipour and Dumitru Baleanu. 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 present two methods for solving a nonlinear system of fractional differential equations within Caputo derivative. Firstly, we derive operational matrices for Caputo fractional derivative and for Riemann-Liouville fractional integral by using the Bernstein polynomials (BPs). In the first method, we use the operational matrix of Caputo fractional derivative (OMCFD), and in the second one, we apply the operational matrix of Riemann-Liouville fractional integral (OMRLFI). The obtained results are in good agreement with each other as well as with the analytical solutions. We show that the solutions approach to classical solutions as the order of the fractional derivatives approaches 1.

#### 1. Introduction

Differential equations of fractional order have been subjected to many studies due to their frequent appearance in various applications in fluid mechanics, viscoelasticity, biology, physics, engineering, and so on. Recently, a large amount of literature was developed regarding the application of fractional differential equations in nonlinear dynamics (see, e.g., [1–11] and the references therein). Thus, a huge attention has been given to the solution of fractional ordinary differential equations, integral equations, and fractional partial differential equations of physical interest. As it is known, there exists no method that yields an exact solution for fractional differential equations. Various methods have been proposed in order to solve the fractional differential equations. These methods include the homotopy perturbation method [12–15], Adomian’s decomposition method [16–20], variation iteration method [12–14, 21–23], homotopy analysis method [24], differential transform method [25], operational matrices [26–28], and nonstandard finite difference scheme [29].

In this paper, we investigate the nonlinear system of fractional differential equations as and the initial condition where and . Also, are multivariable polynomial functions.

The structure of the paper is given later. In Section 2, we present some preliminaries and properties in fractional calculus and Bernstein polynomials. In Section 3, we make operational matrices for product, power, Caputo fractional derivative, and Riemann-Liouville fractional integral by BPs. In Section 4, we apply two methods for solving nonlinear system of fractional differential equations by BPs. In Section 5, numerical examples are simulated to demonstrate the high performance of the proposed method. Conclusions are presented in Section 6.

#### 2. Basic Tools

In this section, we recall some basic definitions and properties of the fractional calculus and Bernstein polynomials.

*Definition 1 (see [2, 7, 10]). *The Riemann-Liouville fractional integral operator of order , of a function , is defined as
and for , , , , the fractional derivative of in the Caputo sense is defined as
where for we have
Also, if , and , then

*Definition 2 (see [30]). *The Bernstein polynomials (BPs) of th degree are defined on the interval as follows:

Lemma 3. *One can write , where is a matrix upper triangular, , and .*

*Proof. *(see [26]).

*Definition 4. *We can define the dual matrix on the basis of Bernstein polynomials of th degree as follows:
where

Lemma 5. *Let be a Hilbert space with the inner product and . Then, we can find the unique vector such that is the best approximation of from space . Moreover, one can get such that .*

*Proof . *(see [31]).

Lemma 6. * Suppose that the function is times continuously differentiable . If is the best approximation out of , then
**
where . Also, if , then the error bound vanishes.*

*Proof . *(see [32]).

#### 3. Operational Matrices of Bernstein Polynomials

In Section 3, we recall the operational matrices for product, power, Caputo fractional derivative and Riemann-Liouville fractional integral by BPs.

Lemma 7. *Suppose that is an arbitrary vector. The operational matrix of product using BPs can be given as follows:
*

*Proof . *(see [27]).

Corollary 8. *Suppose that , , and is the operational matrix of product using BPs for vector . One can get the approximate function for using BPs as follows:
*

*Proof. *By using Lemma 7, it is clear.

Corollary 9. *Suppose that and is the operational matrix of product using BPs for vector . One can get the approximate function for , using BPs as follows:
**
where .*

*Proof . *(see [26]).

Theorem 10. *One can get BPs operational matrix from order for the Caputo fractional derivative as follows:
*

*Proof. *See [26] for details.

Theorem 11. *One can obtain the operational matrix from order for the Riemann-Liouville fractional integral on the basis of BPs from order as
*

*Proof. *See [28] for details.

#### 4. Solving System of Fractional Differential Equations

In this section, we use two methods for solving system of fractional differential equations. In the first method, we use the operational matrix for Caputo fractional derivative (OMCFD), and in the second method, we apply the operational matrix for Riemann-Liouville fractional integral (OMRLFI).

##### 4.1. Solving the Problem by OMCFD

Using Lemma 5, we can approximate the functions as follows: where .

From (17) and (15) we can write Therefore, problem (1) and (2) reduces to the following problem: and the initial condition Now, using Lemma 5 we can approximate all of the known functions in the system (19). Then, by using Lemma 7 and Corollaries 8 and 9, since functions are polynomial, we obtain the following approximations: where .

Also, for each , by using tau method [33] we can generate algebraic equations from (19) and (21) as follows and from (23) we set .

Finally, problem (1) and (2) has been reduced to the system of algebraic equations The aforementioned system can be solved for by Newton’s iterative method. Then, we get the approximate value of the functions from (17).

##### 4.2. Solving the Problem by OMRLFI

This method consists of two steps.

*Step 1. * Initial conditions are used to reduce a given initial-value problem to a problem with zero initial conditions. Therefore we have a modified system, incorporating the initial values.

*Step 2. * The BPs operational matrix of Riemann-Liouville fractional integral is used to transform the problem into a system of algebraic equations.

Now, from (2) we define where , are the new unknown functions.

Substituting (24) in (1) and (2), we have the following system: and the initial condition where and are multivariable polynomial functions. We use the following approximation: where are unknown vectors. From (7), (27), and Theorem 11, we can write So, by (27) and (28), problem (25) and (26) reduces to the following problem: As we saw in the previous section, we can obtain the following approximations: where . So, from (29) and (30) we have Therefore, we have reduced problem (1) and (2) to the system of algebraic equations as follows: where this system can be solved for by Newton’s iterative method. Finally we obtain the approximate of the functions by

#### 5. Examples

To demonstrate the applicability and to validate the numerical scheme, we apply the present method for the following examples.

*Example 12. *Consider the following linear system of fractional differential equations [24, 25]:
with initial condition
For this problem we have the exact solution in the case of as
We solved this problem by OMCFD and OMRLFI. Figures 1 and 2 show the approximate solutions of and , respectively, as a function of time for , for different values of , . The results show that numerical solutions are in good agreement with each other, in both methods. Also, these figures show that as , approach close to 1, the numerical solutions approach to the solutions for as expected. In Figures 3 and 4, we see the absolute error of both methods, for , . In these figures, we can see that obtained results using the presented methods agree well with the analytical solutions for .

*Example 13. *Let us consider the following nonlinear fractional system [24] as follows:
such that

The exact solution of this system, when , is Figures 5 and 6 show the approximate solutions of and , respectively, for different values of , by OMCFD and OMRLFI. We conclude that as , approach close to 1, the numerical solutions approach solutions for as expected. Furthermore, in both methods, the results agree well with each other. Figures 7 and 8 show that, the absolute error of obtained results for and using OMCFD and OMRLFI is in good agreement with the exact solution.

*Example 14. * Consider the nonlinear system of fractional differential equations [24]:
with the initial conditions given by
The exact solution of this system, when , becomes
We can see the approximate solutions of and , by OMCFD and OMRLFI for and different values of , and , in Figures 9, 10, and 11. These figures show that, when , , and approach close to 1, the numerical solutions approach the solutions for as expected. In Figures 9–11, we observe that results of OMCFD and OMRLFI overlap. In Figures 12, 13, and 14, we see the absolute error of the obtained results for and in both methods.

#### 6. Conclusion

In this paper, we get operational matrices of the product, Caputo fractional derivative, and Riemann-Liouville fractional integral by Bernstein polynomials. Then by using these matrices, we proposed two methods that reduced the nonlinear systems of fractional differential equations to the two system of algebraic equations that can be solved easily. Finally, numerical examples are simulated to demonstrate the high performance of the proposed method. We saw that the results of both methods were in good agreement with each other, and the classical solutions were recovered when the order of the fractional derivative goes to 1.

#### References

- J.-H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 167, no. 1-2, pp. 57–68, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Podlubny,
*Fractional Differential Equations*, Academic Press, New York, NY, USA, 1999. View at MathSciNet - R. Hilfer,
*Applications of Fractional Calculus in Physics*, World Scientific, Singapore, 2000. - X. Gao and J. Yu, “Synchronization of two coupled fractional-order chaotic oscillators,”
*Chaos, Solitons and Fractals*, vol. 26, no. 1, pp. 141–145, 2005. View at Publisher · View at Google Scholar · View at Scopus - J. G. Lu, “Chaotic dynamics and synchronization of fractional-order Arneodo's systems,”
*Chaos, Solitons & Fractals*, vol. 26, no. 4, pp. 1125–1133, 2005. View at Publisher · View at Google Scholar - J. G. Lu and G. Chen, “A note on the fractional-order Chen system,”
*Chaos, Solitons & Fractals*, vol. 27, no. 3, pp. 685–688, 2006. View at Publisher · View at Google Scholar - A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and Applications of Fractional Differential Equations*, Elsevier, San Diego, Calif, USA, 2006. View at MathSciNet - D. Baleanu, O. G. Mustafa, and R. P. Agarwal, “An existence result for a superlinear fractional differential equation,”
*Applied Mathematics Letters*, vol. 23, no. 9, pp. 1129–1132, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Baleanu, O. G. Mustafa, and R. P. Agarwal, “On the solution set for a class of sequential fractional differential equations,”
*Journal of Physics A*, vol. 43, no. 38, Article ID 385209, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus Models and Numerical Methods*, Series on Complexity, Nonlinearity and Chaos, World Scientific, Hackensack, NJ, USA, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Bhalekar, V. Daftardar-Gejji, D. Baleanu, and R. L. Magin, “Transient chaos in fractional Bloch equations,”
*Computers & Mathematics with Applications*, vol. 64, no. 10, pp. 3367–3376, 2012. View at Publisher · View at Google Scholar - S. Momani and Z. Odibat, “Numerical approach to differential equations of fractional order,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 96–110, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Momani and Z. Odibat, “Homotopy perturbation method for nonlinear partial differential equations of fractional order,”
*Physics Letters A*, vol. 365, no. 5-6, pp. 345–350, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Momani and Z. Odibat, “Numerical comparison of methods for solving linear differential equations of fractional order,”
*Chaos, Solitons & Fractals*, vol. 31, no. 5, pp. 1248–1255, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Momani and K. Al-Khaled, “Numerical solutions for systems of fractional differential equations by the decomposition method,”
*Applied Mathematics and Computation*, vol. 162, no. 3, pp. 1351–1365, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Jafari and V. Daftardar-Gejji, “Solving a system of nonlinear fractional differential equations using Adomian decomposition,”
*Journal of Computational and Applied Mathematics*, vol. 196, no. 2, pp. 644–651, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Lesnic, “The decomposition method for Cauchy advection-diffusion problems,”
*Computers & Mathematics with Applications*, vol. 49, no. 4, pp. 525–537, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Lesnic, “The decomposition method for initial value problems,”
*Applied Mathematics and Computation*, vol. 181, no. 1, pp. 206–213, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 7, no. 1, pp. 27–34, 2006. View at Google Scholar · View at Scopus - S. Momani and Z. Odibat, “Analytical approach to linear fractional partial differential equations arising in fluid mechanics,”
*Physics Letters A*, vol. 355, no. 4-5, pp. 271–279, 2006. View at Publisher · View at Google Scholar · View at Scopus - V. Daftardar-Gejji and H. Jafari, “An iterative method for solving nonlinear functional equations,”
*Journal of Mathematical Analysis and Applications*, vol. 316, no. 2, pp. 753–763, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Zurigat, S. Momani, Z. Odibat, and A. Alawneh, “The homotopy analysis method for handling systems of fractional differential equations,”
*Applied Mathematical Modelling*, vol. 34, no. 1, pp. 24–35, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. S. Ertürk and S. Momani, “Solving systems of fractional differential equations using differential transform method,”
*Journal of Computational and Applied Mathematics*, vol. 215, no. 1, pp. 142–151, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Alipour, D. Rostamy, and D. Baleanu, “Solving multi-dimensional FOCPs with inequality constraint by BPs operational matrices,”
*Journal of Vibration and Control*, 2012. View at Publisher · View at Google Scholar - D. Rostamy and K. Karimi, “Bernstein polynomials for solving fractional heat- and wave-like equations,”
*Fractional Calculus and Applied Analysis*, vol. 15, no. 4, pp. 556–571, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - D. Rostamy, M. Alipour, H. Jafari, and D. Baleanu, “Solving multi-term orders fractional differential equations by operational matrices of BPs with convergence analysis,”
*Romanian Reports in Physics*, vol. 65, no. 2, 2013. View at Google Scholar - S. Momani, A. Abu Rqayiq, and D. Baleanu, “A nonstandard finite difference scheme for two-sided space-fractional partial differential equations,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 22, no. 4, Article ID 1250079, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - E. W. Cheney,
*Introduction to Approximation Theory*, AMS Chelsea Publishing, Providence, RI, USA, 2nd edition, 1982. View at MathSciNet - E. Kreyszig,
*Introduction Functional Analysis with Applications*, John Wiley & Sons, New York, NY, USA, 1978. View at MathSciNet - M. Alipour and D. Rostamy, “Bernstein polynomials for solving Abel's integral equation,”
*The Journal of Mathematics and Computer Science*, vol. 3, no. 4, pp. 403–412, 2011. View at Google Scholar - C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang,
*Spectral Methods in Fluid Dynamic*, Prentice-Hall, Englewood Cliffs, NJ, USA, 1988. View at MathSciNet