Next Article in Journal
Degenerate Daehee Numbers of the Third Kind
Next Article in Special Issue
On a Length Problem for Univalent Functions
Previous Article in Journal
Oriented Algebras and the Hochschild Cohomology Group
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Operational Matrix of Fractional Derivatives to Solve Systems of Fractional Differential Equations via Legendre Wavelets

1
Department of Mathematical Engineering, Yildiz Technical University, Istanbul 34200, Turkey
2
Yildiz Technical University, Istanbul 34200, Turkey
*
Author to whom correspondence should be addressed.
Submission received: 2 October 2018 / Revised: 29 October 2018 / Accepted: 31 October 2018 / Published: 5 November 2018

Abstract

:
This paper introduces a new numerical approach to solving a system of fractional differential equations (FDEs) using the Legendre wavelet operational matrix method (LWOMM). We first formulated the operational matrix of fractional derivatives in some special conditions using some notable characteristics of Legendre wavelets and shifted Legendre polynomials. Then, the system of fractional differential equations was transformed into a system of algebraic equations by using these operational matrices. At the end of this paper, several examples are presented to illustrate the effectivity and correctness of the proposed approach. Comparing the methodology with several recognized methods demonstrates that the advantages of the Legendre wavelet operational matrix method are its accuracy and the understandability of the calculations.

1. Introduction

Differential and integral operators are the basis of mathematical models, and they are also used as a means of understanding the working principles of natural and artificial systems. Therefore, differential and integral equations are of great importance both theoretically and practically. Such equations have a wide range of applications, including in the physical sciences (such as in physics and engineering) as well as in social science. Systems of differential equations, as differential equations, are often used in issues such as theories of elasticity, dynamics, fluid mechanics, oscillation, and quantum dynamics.
Interest in differential and integral operators has led to the exploration of fractional differential and integral operators by examining these issues further in depth. Owing to a question, the origin of fractional calculus arose in a message from Leibniz to L’Hôpital in 1695. Fractional calculus has received attention in recent years due to its ability to simplify numerous physical, engineering, and economics phenomena such as the fluid dynamic traffic model, damping laws, continuum and statistical mechanics, diffusion processes, solid mechanics, control theory, colored noise, viscoelasticity, electrochemistry, and electromagnetism, among others.
Because a variety of solutions of fractional differential equations (FDEs) cannot be found analytically, numerical and approximate methods are needed. There are a lot of techniques that have been studied by many researchers in solving FDEs and the system of such equations numerically. Some of these techniques are the Adomian decomposition method presented by Song and Wang [1], the collocation method, the improved operational matrix method [2,3,4], the perturbation iteration method introduced by Şenol and Dolapçı [5], the computational matrix method illustrated by Khader et al. [6], the differential transform method demonstrated by Ertürk and Momani [7], the variational iteration method, the Laplace transform method given by Gupta et al. [8], and the fractional complex transform method studied by Ghazanfari and Ghazanfari [9], among others. Kilbas et al. [10] inclusively examined fractional differential and fractional integro-differential equations. In addition, numerical solutions of FDEs and the system of such equations have been presented using the Legendre polynomial operational matrix method [11], Bernstein operational matrix method [12], Genocchi operational matrix method [13], Jacobi operational matrix method [14], Chebyshev wavelet operational matrix method [15], polynomial least squares method (PLSM) [16], Legendre wavelet-like operational matrix method (LWPT) [17], and the Genocchi wavelet-like operational matrix method [18].
This paper focuses on the numerical analysis of a system of fractional order differential equations using the Legendre wavelet operational matrix method. The most important advantage of the proposed method is that it presents an understandable procedure to reduce FDEs and the system of such equations to a system of algebraic equations. First, we begin by presenting some basic definitions and fundamental relations in Section 2 and Section 3, respectively. Then, in Section 4, the operational matrix of the fractional derivate is natively formulated to linear and nonlinear systems of fractional differential equations. Section 5 presents five illustrative examples that were tested with the introduced method. Finally, the last section includes the conclusions.

2. Basic Definitions

The Liouville_Caputo fractional_order derivative, shifted Legendre polynomials, and Legendre wavelets are defined below [19,20].
Definition 1.
The Liouville_Caputo fractional derivative of u is defined as [19]
D α u ( t ) = 1 Γ ( n α ) 0 t u ( n ) ( ξ ) ( t ξ ) α + 1 n d ξ , n 1 < α n ,   n N .
Some characteristics of the Liouville_Caputo fractional derivative are as follows:
D α C = 0 ,
where C is a constant. In addition, there is
D α t β = { 0 , β N 0   a n d   β < α Γ ( β + 1 ) Γ ( β + 1 α ) x β α , β N 0   a n d   β α   o r   β N   a n d   β > α ,
in which α and α respectively imply that the largest integer is less than or equal to α , and the smallest integer is greater than or equal to α .
The Liouville_Caputo fractional order derivative is a linear operation of the integer order derivative
D α ( η u ( t ) + ζ v ( t ) ) = η D α u ( t ) + ζ D α v ( t ) ,
where η and ζ are constant.
Definition 2.
Let a and b respectively be the parameters of dilation and translation of a single function called the mother wavelet. If a and b change continuously, then we obtain the following family of continuous wavelets [21,22]:
ψ a b ( t ) = | a | 1 / 2 ψ ( t b a ) , a , b R , a 0 .
Definition 3.
Let P m ( t ) imply the shifted Legendre polynomials of order m . Then P m ( t ) can be formulated as [21]
P m ( t ) = k = 0 m ( 1 ) m + k ( m + k ) ! ( m k ) ! t k ( k ! ) 2 ,
and the orthogonality condition is
0 1 P m ( t ) P n ( t ) d t = { 1 2 m + 1 , f o r   m = n 0 , f o r   m n .
Definition 4.
Let n and k be any positive integer, m be the order of shifted Legendre polynomials, and t be the normalized time. Then the Legendre wavelets ψ n m ( t ) = ψ ( k , n , m , t ) are defined on the interval [ 0 , 1 ] by [21,22].
ψ n m ( t ) = { 2 k + 1 2 m + 1 2 P m ( 2 k t n ) , n 2 k t n + 1 2 k 0 , o t h e r w i s e ,
where m = 0 , 1 , , M ;   n = 0 , 1 , , ( 2 k 1 ) . The coefficient m + 1 2 is for orthonormality.
Definition 5.
Let u ( t ) and v ( t ) be functions defined over [ 0 , 1 ] and then expanded in the terms of the Legendre wavelet as [21,22]
u ( t ) = n = 0 m = 0 c n m ψ n m ( t ) ,
where c n m = ( u ( t ) , ψ n m ( t ) ) , in which ( . , . ) implies the inner product. If the infinite series in Equation (9) is truncated, then it can be expressed as
u ( t ) n = 0 2 k 1 m = 0 M c n m ψ n m ( t ) = C T ψ ( t ) ,
where C and ψ ( t ) are matrices, as presented by
C = [ c 0 , 0 , c 0 , 1 , , c 0 , M , , c 2 , M , , c 2 k 1 , 0 , c 2 k 1 , 1 , , c 2 k 1 , M ] T ψ = [ ψ 0 , 0 , ψ 0 , 1 , , ψ 0 , M , , ψ 2 , M , , ψ 2 k 1 , 0 , ψ 2 k 1 , 1 , , ψ 2 k 1 , M ] T .

3. Fundamental Relations

Saadatmandi and Dehghan [11] derived the operational matrix of a fractional derivative by using shifted Legendre polynomials. In this section, we show how we derived the Legendre wavelet operational matrix of fractional derivatives in some special conditions by drawing from Saadatmandi and Dehghan [11]. Additionally, the theorem and corollary related to the Legendre wavelet operational matrix of derivatives illustrated by Mohammadi [21] are cited here as follows.
Theorem 1.
Let ψ ( t ) be the Legendre wavelet vector introduced in Equation (8). Then ψ ( t ) is expressed as [21,22]
d ψ ( t ) d t = D ψ ( t ) ,
where D is the 2 k ( M + 1 ) operational matrix of the derivative, which can be stated as
D = [ U O O O U O O O U ] ,
where U is an ( M + 1 ) ( M + 1 ) matrix and its ( r , s ) th element is written as
U r , s = { 2 k + 1 ( 2 r 1 ) ( 2 s 1 ) , r = 2 , , ( M + 1 ) ,   s = 1 , , r 1   and   ( r + s )   odd 0 , otherwise .
Corollary 1.
Using Equation (12), the operational matrix for the n th derivative can be stated as [21]
d n ψ ( t ) d t n = D n ψ ( t ) ,
where D n is the n th power of matrix D .
Lemma 1.
Let ψ ( t ) be the Legendre wavelets vector introduced in Equation (8). Assuming that k = 0 , then
D α ψ r ( t ) = 0 ,   r = 0 , 1 , , α 1 ,   α > 0 .
Proof. 
The desired result can be obtained by using Equations (2) and (4) in Equation (8). □
Theorem 2.
Let ψ ( t ) be the Legendre wavelets vector introduced in Equation (8). Supposing that k = 0 and α > 0 , then
D α ψ ( t ) D ( α ) ψ ( t ) ,
where D ( α ) is the ( M + 1 ) x ( M + 1 ) operational matrix of the fractional derivative of the order α > 0 ,   N 1 < α N in the Liouville_Caputo sense and can be stated as
D ( α ) = ( 0 0 0 0 0 0 h = α α ξ α , 0 , h h = α α ξ α , 1 , h h = α α ξ α , m , h h = α r ξ r , 0 , h h = α r ξ r , 1 , h h = α r ξ r , m , h h = α m ξ m , 0 , h h = α m ξ m , 1 , h h = α m ξ m , m , h ) ,
where ξ r , s , h is written as
ξ r , s , h = 2 r + 1 2 s + 1 l = 0 s ( 1 ) r + s + h + l ( r + h ) ! ( s + l ) ! ( r h ) ! h ! Γ ( h α + 1 ) ( s l ) ! ( l ! ) 2 ( h + l α + 1 ) .
Consider in D ( α ) that the first α rows are all zero.
Proof. 
Presume that ψ r ( t ) is the r th element of the vector ψ ( t ) introduced in Equation (11), where r = n M + ( m + 1 ) ,   m = 0 , 1 , , M ,   n = 0 , 1 , , ( 2 k 1 ) . Then ψ r ( t ) can be stated as
ψ r ( t ) = 2 k + 1 2 r + 1 2 P r ( 2 k t n ) χ [ n 2 k , n + 1 2 k ] .
Accepting that k = 0 , and by using the shifted Legendre polynomial, we obtain
ψ r ( t ) = 2 r + 1 2 h = 0 r ( 1 ) r + h ( r + h ) ! ( r h ) ! ( h ! ) 2 t h χ [ 0 , 1 ] .
If we use Equations (3), (4), and (21), then we have
D ( α ) ψ r ( t ) = 2 r + 1 2 h = 0 r ( 1 ) r + h ( r + h ) ! ( r h ) ! ( h ! ) 2 D α ( t h ) χ [ 0 , 1 ] = 2 r + 1 h = α r ( 1 ) r + h ( r + h ) ! ( r h ) ! ( h ! ) Γ ( h α + 1 ) t h α χ [ 0 , 1 ] , r = α , , m .
Approximating t h α by ( m + 1 ) terms of the Legendre wavelets, then we obtain
t h α s m b h , s ψ s ( t ) ,
where
b h , s = 0 1 t h α ψ s ( t ) d t = 2 s + 1 2 l = 0 s ( 1 ) s + l ( s + l ) ! ( s l ) ! ( l ! ) 2 0 1 t h + l α d t = 2 s + 1 l = 0 s ( 1 ) s + l ( s + l ) ! ( s l ) ! ( l ! ) 2 ( h + l α + 1 ) .
Utilizing Equations (22) and (24), we get
D α ψ r ( t ) 2 r + 1 h = α r s = 0 m ( 1 ) r + h ( r + h ) ! ( r h ) ! ( h ! ) Γ ( h α + 1 ) b h , s ψ s ( t ) χ [ 0 , 1 ] = s = 0 m ( h = α r ξ r , s , h ) ψ s ( t ) χ [ 0 , 1 ] , r = α , , m ,
in which ξ r , s , h is presented in Equation (19). In addition, if we use Lemma 1, then we can write
D α ψ r ( t ) = 0 , r = 0 , 1 , , α 1 , α > 0 .
Combining Equations (25) and (27), the result can be obtained. □

4. Solving Systems of Fractional Order Differential Equations

In this section, the Legendre wavelet operational matrix method was implemented to obtain the numerical solution of the system of fractional order differential equations. Consider a system of fractional differential equations as follows:
D η 1 u 1 ( t ) = U 1 ( t , u 1 , u 2 , , u m ) , D η 2 u 2 ( t ) = U 2 ( t , u 1 , u 2 , , u m ) , D η n u m ( t ) = U m ( t , u 1 , u 2 , , u m ) ,
where U i is a linear/nonlinear function of t , u 1 , u 2 , , u m , D η i is the derivative of u i with the order of η i in the Liouville–Caputo sense and N 1 η i < N , and they are subjected to the initial conditions
u 1 ( t 0 ) = u 10 , d u 1 d t ( t 0 ) = u 11 , d 2 u 1 d t 2 ( t 0 ) = u 12 , , d n 1 u 1 d t n 1 ( t 0 ) = u 1 ( n 1 ) u 2 ( t 0 ) = u 20 , d u 2 d t ( t 0 ) = u 21 , d 2 u 2 d t 2 ( t 0 ) = u 22 , , d n 1 u 2 d t n 1 ( t 0 ) = u 2 ( n 1 ) u m ( t 0 ) = u m 0 , d u m d t ( t 0 ) = u m 1 , d 2 u m d t 2 ( t 0 ) = u m 2 , , d n 1 u m d t n 1 ( t 0 ) = u m ( n 1 ) .
First of all, approximating u 1 ( t ) , u 2 ( t ) , , u m ( t ) and D η 1 u 1 ( t ) , D η 2 u 2 ( t ) , , D η n u m ( t ) , we obtain
u 1 ( t ) n = 0 2 k 1 m = 0 M c 1 n , m ψ n , m = C 1 T ψ ( t ) u 2 ( t ) n = 0 2 k 1 m = 0 M c 2 n , m ψ n , m = C 2 T ψ ( t ) u m ( t ) n = 0 2 k 1 m = 0 M c n n , m ψ n , m = C m T ψ ( t ) ,
where C i ,   i = 1 , 2 , , m is an unknown vector and ψ ( t ) is the vector introduced in Equation (8). If we utilize Equation (17), then we have
D η 1 u 1 ( t ) C 1 T D ( η 1 ) ψ ( t ) D η 2 u 2 ( t ) C 2 T D ( η 2 ) ψ ( t ) D η n u m ( t ) C m T D ( η n ) ψ ( t ) .
Substituting Equations (29) and (30) into Equation (27), we obtain
R 1 ( t ) = C 1 T D ( η 1 ) ψ ( t ) U 1 ( t , C 1 T ψ ( t ) , C 2 T ψ ( t ) , , C m T ψ ( t ) ) R 2 ( t ) = C 2 T D ( η 2 ) ψ ( t ) U 2 ( t , C 1 T ψ ( t ) , C 2 T ψ ( t ) , , C m T ψ ( t ) ) R m ( t ) = C m T D ( η n ) ψ ( t ) U m ( t , C 1 T ψ ( t ) , C 2 T ψ ( t ) , , C m T ψ ( t ) )
If U i is a linear function of t , u 1 , u 2 , , u m , then we produce 2 k ( M + 1 ) m n linear equations by implementing
0 1 ψ j ( t ) R i ( t ) d t = 0 ,   j = 1 , , 2 k ( M + 1 ) m n ,   i = 1 , 2 , , m .
Also, by substituting the initial conditions in Equation (28) into Equation (30), then we obtain
u 1 ( t 0 ) C 1 T ψ ( t 0 ) = u 10 , d u 1 d t ( t 0 ) C 1 T D ψ ( t 0 ) = u 11 , , d n 1 u 1 d t n 1 ( t 0 ) C 1 T D n 1 ψ ( t 0 ) = u 1 ( n 1 ) u 2 ( t 0 ) C 2 T ψ ( t 0 ) = u 20 , d u 2 d t ( t 0 ) C 2 T D ψ ( t 0 ) = u 21 , , d n 1 u 2 d t n 1 ( t 0 ) C 2 T D n 1 ψ ( t 0 ) = u 2 ( n 1 ) u m ( t 0 ) C m T ψ ( t 0 ) = u m 0 , d u m d t ( t 0 ) C m T D ψ ( t 0 ) = u m 1 , , d n 1 u m d t n 1 ( t 0 ) C m T D n 1 ψ ( t 0 ) = u m ( n 1 ) .
A 2 k ( M + 1 ) set of linear equations is generated by combining Equations (32) and (33). The solution of these linear equations can be obtained for unknown coefficients of the vector C . Consequently, u 1 ( t ) , u 2 ( t ) , , u m ( t ) , introduced in Equation (27), can be computed.
If U i is a nonlinear function of t , u 1 , u 2 , , u m , then we first compute R 1 ( t ) , R 2 ( t ) , , R m ( t ) at 2 k ( M + 1 ) m n points and, for a better result, use the first 2 k ( M + 1 ) m n roots of shifted Legendre P 2 k ( M + 1 ) ( t ) . Then these equations, collectively with Equation (33), produce 2 k ( M + 1 ) nonlinear equations. The solution of these nonlinear equations can be obtained by employing Newton’s iterative method. Consequently, u 1 ( t ) , u 2 ( t ) , , u m ( t ) , introduced in Equation (27), can be computed.

5. Illustrative Examples

In this section, to show the applicability and powerfulness of the introduced method, we present the solutions to five linear and nonlinear systems of fractional order differential equations.
Example 1.
We first considered the following linear system of fractional differential equations [7,8]:
D α u ( t ) = u ( t ) + v ( t ) D α v ( t ) = u ( t ) + v ( t ) ,
subject to
u ( 0 ) = 0 , v ( 0 ) = 1 .
The exact solution of this system when α = 1 is known to be
u ( t ) = e t sin t , v ( t ) = e t cos t .
This example was examined for M = 2 , k = 0 , and α = 0.9 , 0.7 , 0.5 . When the obtained results were matched against the exact solution when α = 1 , as demonstrated in Figure 1, we can clearly observe that when α approached 1, our results approached the exact solution. We also solved this problem by using Legendre polynomial operational matrix method (LPOMM), and we compared the results with the LWOMM. The numerical computations for u(t) and v(t) when α = 0.9 are revealed in Table 1 and Table 2.
Example 2.
We considered the following nonlinear system of fractional differential equations [13]:
D 3 2 u ( t ) = 8 u ( t ) + v 2 ( t ) 4 t 6 + 4 t 3 + 8 t 3 2 π 1 D 1 2 v ( t ) = t 2 D u ( t ) + v ( t ) 3 t 4 2 t 3 + 32 t 5 2 5 π 1
u ( 0 ) = 0 , v ( 0 ) = 1 , u ( 1 ) = 1 , v ( 1 ) = 3 , u ( 0 ) = 0 , u ( 1 ) = 3 .
The exact solution of this system is known to be
u ( t ) = t 3 , v ( t ) = 2 t 3 + 1
Using the parameters M = 3 and k = 0 , we applied both the proposed method and the LPOMM to solve this problem and show that our approach is more efficient and useful. Our numerical results supported the idea that our solution approaches the exact solution more than the approximate solution LPOMM. Comparisons of the approximate and exact solutions are presented in Table 3 and Table 4.
Example 3.
We considered the following nonlinear system of fractional differential equations with the initial conditions [8]
D α u ( t ) = u ( t ) 2 D α v ( t ) = u 2 ( t ) + v ( t ) ,
u ( 0 ) = 1 , v ( 0 ) = 0 .
The exact solution of this system when α = 1 is known to be
u ( t ) = e ( t 2 ) , v ( t ) = t e t .
The parameters M = 2 ,   k = 0 , and α = 0.5 , 0.7 , 0.9 were utilized. A comparison of our results and the exact solution when α = 1 is displayed in Figure 2. The figures support that when α approximated 1, our results approximated the exact solution. We also solved this problem by using the LPOMM, and compared the results to the LWOMM. Finally, we present the numerical computations for u(t) and v(t) when α = 0.9 in Table 5 and Table 6.
Example 4.
We considered the following nonlinear system of FDEs with initial conditions [13]
D α u ( t ) = 1002 u ( t ) + 1000 v 2 ( t ) D α v ( t ) = u ( t ) v ( t ) v 2 ( t )
u ( 0 ) = 1 , v ( 0 ) = 1
The exact solution of this system when α = 1 is known to be
u ( t ) = e 2 t , v ( t ) = e t
This example was analyzed for M = 3 , k = 0 , and α = 0.9 , 0.7 , 0.5 . When the obtained results were matched against the exact solution when α = 1 , as demonstrated in Figure 3, we can clearly observe that when α approached 1, our results approached the exact solution. We also solved this problem by using the LPOMM, and compared the results with the LWOMM. The numerical computations for u(t) and v(t) when α = 0.99 are also revealed in Table 7 and Table 8.
Example 5.
We considered the following fractional order Brusselator system [16,17]:
D α u ( t ) = 2 u ( t ) + u 2 ( t ) v ( t ) D α v ( t ) = u ( t ) u 2 ( t ) v ( t )
u ( 0 ) = 1 , v ( 0 ) = 1 .
The approximate solutions of this system when α = 1 and α = 0.98 were presented by Chang and Isah using the LWPT [17] and by Bota and Caruntu using the PLSM [16]. These solutions when α = 98 are given by
u L W P T ( t ) = 1 1.0791 t + 0.2711 t 2 0.0638 t 3 , v L W P T ( t ) = 1 + 0.0151 t + 0.4185 t 2 0.2624 t 3 u P L S M ( t ) = 1 1.08655 t + 0.311138 t 2 + 0.0243682 t 3 , v P L S M ( t ) = 1 + 0.0349127 t + 0.333424 t 2 0.184414 t 3 .
The parameters M = 2 , k = 0 , and α = 0.98 were used. A comparison of our results to the approximate solutions introduced by Bota and Caruntu [16] and Chang and Isah [17] when α = 0.98 is displayed in Figure 4. Finally, we also present the numerical computations for u(t) and v(t) when α = 0.98 in Table 9 and Table 10.

6. Conclusions

In this paper, a system of fractional_order differential equations was examined by drawing from a new operational matrix of the fractional derivative in some special conditions. We also systematized a very operational algorithm in order to attain the solution of the linear and nonlinear systems of fractional differential equations in Maple. All numerical results and graphical presentations generated by Maple affirmed that the Legendre wavelet operational matrix method is very effective and applicable. As the next step, the method introduced in this paper can be applied to fractional partial differential equations and the system of such equations, fractional integral equations and the system of such equations, and fractional integro-differential equations. These equations are at least as important as fractional differential equations and they are very significant in science, engineering, and technology.

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Song, L.; Wang, W. A new improved Adomian decomposition method and its application to fractional differential equations. Appl. Math. Model. 2013, 37, 1590–1598. [Google Scholar] [CrossRef]
  2. Balaji, S. Legendre wavelet operational matrix method for solution of fractional order Riccati differential equation. J. Eqypt. Math. Soc. 2015, 23, 263–270. [Google Scholar] [CrossRef]
  3. Rehman, M.; Khan, R.A. The Legendre wavelet method for solving fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4163–4173. [Google Scholar] [CrossRef]
  4. Mohammadi, F.; Hosseini, M.M.; Mohyud-Din, S.T. A new operational matrix for Legendre wavelets and its applications for solving fractional order boundary value problems. Int. J. Syst. Sci. 2011, 32, 7371–7378. [Google Scholar] [CrossRef]
  5. Şenol, M.; Dolapçı, İ.T. On the perturbation-iteration algorithm for fractional differential equations. J. King Saud Univ.-Sci. 2016, 28, 69–74. [Google Scholar] [CrossRef]
  6. Khader, M.M.; Danaf, T.S.; Hendy, A.S. A computational matrix method for solving systems of high order fractional differential equations. Appl. Math. Model. 2013, 37, 4035–4050. [Google Scholar] [CrossRef]
  7. Ertürk, V.S.; Momani, S. Solving systems of fractional differential equations using differential transform method. J. Comput. Appl. Math. 2008, 215, 142–151. [Google Scholar] [CrossRef]
  8. Gupta, S.; Kumar, D.; Singh, J. Numerical study for systems of fractional differential equations via Laplace transform. J. Egypt. Math. Soc. 2015, 23, 256–262. [Google Scholar] [CrossRef]
  9. Ghazanfari, B.; Ghazanfari, A.G. Solving system of fractional differential eqautions by fractional complex transform method. Asian J. Appl. Sci. 2012, 5, 438–444. [Google Scholar] [CrossRef]
  10. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; Elsevier Science: Amsterdam, The Netherlands, 2006; p. 540. [Google Scholar]
  11. Saadatmandi, A.; Dehghan, M. A new operational matrix for solving fractional-order differential equations. Comput. Math. Appl. 2010, 59, 1326–1336. [Google Scholar] [CrossRef]
  12. Alshbool, M.H.T.; Bataineh, A.S.; Hashim, I.; Işık, O.R. Solution of fractional-order differential equations based on the operational matrices of new fractional Bernstein functions. J. King Saud Univ.-Sci. 2017, 29, 1–18. [Google Scholar] [CrossRef]
  13. Isah, A.; Phang, C. New operational matrix of derivative for solving non-linear fractional differential equations via Genocchi polynomials. J. King Saud Univ.-Sci. 2017. [Google Scholar] [CrossRef]
  14. Doha, E.H.; Bhrawy, A.H.; Ezz-Eldien, S.S. A new Jacobi operational matrix: An application for solving fractional differential equations. Appl. Math. Model. 2012, 36, 4931–4943. [Google Scholar] [CrossRef]
  15. Mohammadi, F. Numerical solution of Bagley-Torvik equation using Chebyshev wavelet operational matrix of fractional derivative. Int. J. Adv. Appl. Math. Mech. 2014, 2, 83–91. [Google Scholar]
  16. Bota, C.; Caruntu, B. Approximate analytical solutions of the fractional-order brusselator system using the polynomial least squares method. Corp. Adv. Differ. Equ. 2015. [Google Scholar] [CrossRef]
  17. Chang, P.; Isah, A. Legendre wavelet operational matrix of fractional derivative through wavelet-polynomial transformation and its applications in solving fractional order brusselator system. J. Phys. Conf. Ser. 2016, 693, 012001. [Google Scholar] [CrossRef]
  18. Isah, A.; Phang, C. Genocchi wavelet-like operational matrix and its application for solving non-linear fractional differential equations. Open Phys. 2016, 14, 463–472. [Google Scholar] [CrossRef]
  19. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations, 1st ed.; Wiley: Hoboken, NJ, USA, 1993; p. 384. [Google Scholar]
  20. Oldham, K.B.; Spanier, J. The Fractional Calculus, 1st ed.; Academic Press: New York, NY, USA, 1974; p. 234. [Google Scholar]
  21. Mohammadi, F.; Hosseini, M.M. A new Legendre wavelet operational matrix of derivative and its applications in solving the singular ordinary differential equations. J. Frankl. Inst. 2011, 348, 1787–1796. [Google Scholar] [CrossRef]
  22. Mohammadi, F.; Hosseini, M.M.; Mohyud-Din, S.T. Legendre wavelet Galerkin method for solving ordinary differential equations with non-analytic solution. Int. J. Syst. Sci. 2011, 42, 579–585. [Google Scholar] [CrossRef]
Figure 1. Comparison of our solutions and the exact solution when α = 0.9 , 0.7 , 0.5 in Example 1: (a) Our solution u(t); and (b) Our solution v(t).
Figure 1. Comparison of our solutions and the exact solution when α = 0.9 , 0.7 , 0.5 in Example 1: (a) Our solution u(t); and (b) Our solution v(t).
Mathematics 06 00238 g001aMathematics 06 00238 g001b
Figure 2. Comparison of our solutions to the exact solution when α = 0.9 , 0.7 , 0.5 for Example 3: (a) Our solution u(t); and (b) Our solution v(t).
Figure 2. Comparison of our solutions to the exact solution when α = 0.9 , 0.7 , 0.5 for Example 3: (a) Our solution u(t); and (b) Our solution v(t).
Mathematics 06 00238 g002
Figure 3. Comparison of our solutions to the exact solution when α = 0.9 , 0.7 , 0.5 for Example 4: (a) Our solution u(t); and (b) Our solution v(t).
Figure 3. Comparison of our solutions to the exact solution when α = 0.9 , 0.7 , 0.5 for Example 4: (a) Our solution u(t); and (b) Our solution v(t).
Mathematics 06 00238 g003
Figure 4. Comparison of our solutions to the approximate solution LWPT and the approximate solution PLSM when α = 0.98 for Example 5: (a) Our solution u(t); and (b) Our solution v(t).
Figure 4. Comparison of our solutions to the approximate solution LWPT and the approximate solution PLSM when α = 0.98 for Example 5: (a) Our solution u(t); and (b) Our solution v(t).
Mathematics 06 00238 g004
Table 1. Numerical solutions of u(t) when α = 0.9 attained by the introduced method and the LPOMM for Example 1.
Table 1. Numerical solutions of u(t) when α = 0.9 attained by the introduced method and the LPOMM for Example 1.
t u L W O M M u L P O M M Absolute Error
0.00.3 × 10 9 0.00000000000.3 × 10 9
0.10.14837843300.14837843250.12 × 10 9
0.20.32176452830.32176452770.633 × 10 9
0.30.52015828620.52015828550.65 × 10 9
0.40.74355970670.74355970590.83 × 10 9
0.50.99196878980.99196878900.8 × 10 9
0.61.2653855361.2653855350.53 × 10 9
0.71.5638099441.5638099430.95 × 10 9
0.81.8872420141.8872420140.733 × 10 9
0.92.2356817482.2356817480.62 × 10 9
1.02.6091291442.6091291440.3 × 10 9
Table 2. Numerical solutions of v(t) when α = 0.9 attained by the introduced method and the LPOMM for Example 1.
Table 2. Numerical solutions of v(t) when α = 0.9 attained by the introduced method and the LPOMM for Example 1.
t v L W O M M v L P O M M Absolute Error
0.01.0000000001.0000000000.2 × 10 9
0.11.1522708991.152270900−0.5 × 10 9
0.21.2748018581.2748018580.599 × 10 9
0.31.3675928771.367592878−0.67 × 10 9
0.41.4306439561.430643957−0.13 × 10 8
0.51.4639550941.463955094−0.1 × 10 8
0.61.4675262911.467526293−0.13 × 10 8
0.71.4413575491.441357551−0.171 × 10 8
0.81.3854488661.385448868−0.1461 × 10 8
0.91.2998002441.299800245−0.15 × 10 8
1.01.1844116801.184411682−0.18 × 10 8
Table 3. The numerical results attained by using the introduced method in comparison to the approximate solution LPOMM and the exact solution u(t) in Example 2.
Table 3. The numerical results attained by using the introduced method in comparison to the approximate solution LPOMM and the exact solution u(t) in Example 2.
t Exact Solution u L W O M M u L P O M M
0.00.000−0.12 × 10 9 0.000000000000
0.10.0010.010000000050.001000000000
0.20.0080.020000000160.008000000000
0.30.0270.037500000210.027000000000
0.40.0640.070000000200.064000000000
0.50.1250.125000000010.125000000000
0.60.2160.210000000000.216000000000
0.70.3430.332499999980.343000000000
0.80.5120.499999999960.512000000000
0.90.7290.719999999930.729000000000
1.01.0000.999999999891.000000000000
Table 4. The numerical results attained by using the introduced method in comparison with the approximate solution LPOMM and exact solution v(t) in Example 2.
Table 4. The numerical results attained by using the introduced method in comparison with the approximate solution LPOMM and exact solution v(t) in Example 2.
t Exact Solution v L W O M M v L P O M M
0.01.0001.0000000000.9999999998
0.11.0021.0348413671.165803114
0.21.0161.0575876281.283854751
0.31.0541.0865376681.374911456
0.41.1281.1399903701.459729770
0.51.2501.2362446181.559066242
0.61.4321.3935992961.693677414
0.71.6861.6303532901.884319831
0.82.0241.9648054822.151750038
0.92.4582.4152547582.516724579
1.03.0003.0000000003.000000000
Table 5. Our solutions u(t) when α = 0.9 attained by the presented method and the LPOMM for Example 3.
Table 5. Our solutions u(t) when α = 0.9 attained by the presented method and the LPOMM for Example 3.
t u L W O M M u L P O M M Absolute Error
0.01.0000000001.000000000−0.37 × 10 10
0.11.0648163201.064816320−0.132 × 10 9
0.21.1302485881.130248589−0.1375 × 10 9
0.31.1962968071.1962968070.44 × 10 10
0.41.2629609741.262960975−0.808 × 10 9
0.51.3302410901.330241091−0.532 × 10 9
0.61.3981371561.398137156−0.168 × 10 9
0.71.4666491711.466649171−0.756 × 10 9
0.81.5357771341.535777134−0.1375 × 10 9
0.91.6055210461.6055210460.468 × 10 9
1.01.6758809081.6758809080.163 × 10 9
Table 6. Our solutions v(t) when α = 0.9 attained by the presented method and the LPOMM for Example 3.
Table 6. Our solutions v(t) when α = 0.9 attained by the presented method and the LPOMM for Example 3.
t v L W O M M v L P O M M Absolute Error
0.0−0.1 × 10 9 0.0000000000−0.1 × 10 9
0.10.13817626030.1381762609−0.7 × 10 9
0.20.31336033610.3133603363−0.2 × 10 9
0.30.52555222590.5255522263−0.37 × 10 9
0.40.77475193040.7747519309−0.5 × 10 9
0.51.0609594501.060959450−0.5 × 10 9
0.61.3841747831.384174784−0.8 × 10 9
0.71.7443979321.744397932−0.47 × 10 9
0.82.1416288942.141628895−0.7 × 10 9
0.92.5758676722.5758676720.3 × 10 9
1.03.0471142643.047114264−0.1 × 10 9
Table 7. Numerical solutions of u(t) when α = 0.99 for Example 4.
Table 7. Numerical solutions of u(t) when α = 0.99 for Example 4.
t u L W O M M u L P O M M Absolute Error
0.01.0000000001.000000000−0.4 × 10 10
0.10.81443515290.8144351528−0.21 × 10 10
0.20.66394252330.6639425233−0.15 × 10 9
0.30.54299472290.5429947230−0.34 × 10 9
0.40.44606436360.4460643636−0.42 × 10 9
0.50.36762405680.36762405680
0.60.30214641420.3021464147−0.68 × 10 9
0.70.24410404870.2441040481−0.24 × 10 9
0.80.18796956990.18796956900.53 × 10 9
0.90.12821559200.12821559050.461 × 10 9
1.00.05931472480.05931472220.144 × 10 8
Table 8. Numerical solutions of v(t) when α = 0.99 for Example 4.
Table 8. Numerical solutions of v(t) when α = 0.99 for Example 4.
t v L W O M M v L P O M M Absolute Error
0.01.0000000000.99999999990.79 × 10 10
0.10.90256018370.90256018370.1498 × 10 9
0.20.81526464870.8152646488−0.119 × 10 9
0.30.73711165770.7371116578−0.8 × 10 10
0.40.66709947370.6670994739−0.11 × 10 9
0.50.60422635970.6042263600−0.24 × 10 9
0.60.54749057850.5474905788−0.19 × 10 9
0.70.49589039320.4958903935−0.26 × 10 9
0.80.44842406640.4484240670−0.443 × 10 9
0.90.40408986140.4040898620−0.5898 × 10 9
1.00.36188604080.3618860417−0.719 × 10 9
Table 9. Numerical solutions of u(t) when α = 0.98 obtained by the introduced method, the LWPT, and the PLSM for Example 5.
Table 9. Numerical solutions of u(t) when α = 0.98 obtained by the introduced method, the LWPT, and the PLSM for Example 5.
t u L W O M M u L W P T u P L S M
0.01.0000000001.00000001.0000000000
0.10.89420248260.89473720.8944807482
0.20.79506969160.79451360.7953304656
0.30.70260162680.69894640.7026953614
0.40.61679828830.60765280.6167216448
0.50.53765967610.52025000.5375555250
0.60.46518579020.43635520.4653432112
0.70.39937663060.35558560.4002309126
0.80.34023219730.27755840.3423648384
0.90.28775249020.20189080.2918911978
1.00.24193750950.12820000.2489562000
Table 10. Numerical solutions of v(t) when α = 0.98 obtained by the introduced method, the LWPT, and the PLSM for Example 5.
Table 10. Numerical solutions of v(t) when α = 0.98 obtained by the introduced method, the LWPT, and the PLSM for Example 5.
t v L W O M M v L W P T v P L S M
0.01.0000000001.00000001.000000000
0.11.0080693071.00543261.006641096
0.21.0194799611.01766081.018844188
0.31.0342319591.03511021.035502792
0.41.0523253041.05620641.055510424
0.51.0737599951.07937501.077760600
0.61.0985360321.10304161.101146836
0.71.1266534151.12563181.124562648
0.81.1581121431.14557121.146901552
0.91.1929122171.16128541.167057064
1.01.2310536381.17120001.183922700

Share and Cite

MDPI and ACS Style

Secer, A.; Altun, S. A New Operational Matrix of Fractional Derivatives to Solve Systems of Fractional Differential Equations via Legendre Wavelets. Mathematics 2018, 6, 238. https://0-doi-org.brum.beds.ac.uk/10.3390/math6110238

AMA Style

Secer A, Altun S. A New Operational Matrix of Fractional Derivatives to Solve Systems of Fractional Differential Equations via Legendre Wavelets. Mathematics. 2018; 6(11):238. https://0-doi-org.brum.beds.ac.uk/10.3390/math6110238

Chicago/Turabian Style

Secer, Aydin, and Selvi Altun. 2018. "A New Operational Matrix of Fractional Derivatives to Solve Systems of Fractional Differential Equations via Legendre Wavelets" Mathematics 6, no. 11: 238. https://0-doi-org.brum.beds.ac.uk/10.3390/math6110238

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop