• Views 1,650
• Citations 10
• ePub 56
• PDF 929
`Advances in Mathematical PhysicsVolume 2013, Article ID 954015, 9 pageshttp://dx.doi.org/10.1155/2013/954015`
Research Article

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

1Department of Mathematics, Imam Khomeini International University, P.O. Box 34149-16818, Qazvin, Iran
2Department of Mathematics and Computer Sciences, Cankaya University, 06530 Ankara, Turkey
3Institute of Space Sciences, P.O. Box MG-23, 077125 Magurele-Bucharest, Romania
4Department 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

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., [111] 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 [1215], Adomian’s decomposition method [1620], variation iteration method [1214, 2123], homotopy analysis method [24], differential transform method [25], operational matrices [2628], 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 .

Figure 1: Approximate solutions of for and different values of , by OMCFD and OMRLFI in Example 12.
Figure 2: Approximate solutions of for and different values of , by OMCFD and OMRLFI in Example 12.
Figure 3: Plot of absolute error function for and by OMCFD and OMRLFI in Example 12.
Figure 4: Plot of absolute error function for and by OMCFD and OMRLFI in Example 12.

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.

Figure 5: Approximate solutions of for and different values of , by OMCFD and OMRLFI in Example 13.
Figure 6: Approximate solutions of for and different values of , by OMCFD and OMRLFI in Example 13.
Figure 7: Plot of absolute error function for and by OMCFD and OMRLFI in Example 13.
Figure 8: Plot of absolute error function for and by OMCFD and OMRLFI in Example 13.

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 911, 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.

Figure 9: Approximate solutions of for and different values of , , by OMCFD and OMRLFI in Example 14.
Figure 10: Approximate solutions of for and different values of , , by OMCFD and OMRLFI in Example 14.
Figure 11: Approximate solutions of for and different values of , , by OMCFD and OMRLFI in Example 14.
Figure 12: Plot of absolute error function for and by OMCFD and OMRLFI in Example 14.
Figure 13: Plot of absolute error function for and by OMCFD and OMRLFI in Example 14.
Figure 14: Plot of absolute error function for and by OMCFD and OMRLFI in Example 14.

#### 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

1. 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.
2. I. Podlubny, Fractional Differential Equations, Academic Press, New York, NY, USA, 1999.
3. R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.
4. 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.
5. J. G. Lu, “Chaotic dynamics and synchronization of fractional-order Arneodo's systems,” Chaos, Solitons & Fractals, vol. 26, no. 4, pp. 1125–1133, 2005.
6. 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.
7. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, San Diego, Calif, USA, 2006.
8. 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.
9. 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.
10. 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.
11. 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.
12. 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.
13. 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.
14. 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.
15. 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.
16. 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.
17. 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.
18. D. Lesnic, “The decomposition method for Cauchy advection-diffusion problems,” Computers & Mathematics with Applications, vol. 49, no. 4, pp. 525–537, 2005.
19. D. Lesnic, “The decomposition method for initial value problems,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 206–213, 2006.
20. 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.
21. 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.
22. 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.
23. 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.
24. 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.
25. 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.
26. M. Alipour, D. Rostamy, and D. Baleanu, “Solving multi-dimensional FOCPs with inequality constraint by BPs operational matrices,” Journal of Vibration and Control, 2012.
27. 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.
28. 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.
29. 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.
30. E. W. Cheney, Introduction to Approximation Theory, AMS Chelsea Publishing, Providence, RI, USA, 2nd edition, 1982.
31. E. Kreyszig, Introduction Functional Analysis with Applications, John Wiley & Sons, New York, NY, USA, 1978.
32. 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.
33. C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamic, Prentice-Hall, Englewood Cliffs, NJ, USA, 1988.