Next Article in Journal
AI-Enabled Consensus Algorithm in Human-Centric Collaborative Computing for Internet of Vehicle
Next Article in Special Issue
Dynamical Study of Coupled Riemann Wave Equation Involving Conformable, Beta, and M-Truncated Derivatives via Two Efficient Analytical Methods
Previous Article in Journal
Generalized Support Vector Regression and Symmetry Functional Regression Approaches to Model the High-Dimensional Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations

1
Institute of Mathematical Sciences, Universiti Malaya, Kuala Lumpur 50603, Malaysia
2
Department of Mathematics, AL-Qunfudhah University College, Umm Al-Qura University, Al Qunfudhah 24382, Saudi Arabia
3
Department of Mathematics, Academy of Engineering and Medical Sciences, Khartoum 12045, Sudan
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(6), 1263; https://doi.org/10.3390/sym15061263
Submission received: 24 May 2023 / Revised: 7 June 2023 / Accepted: 12 June 2023 / Published: 15 June 2023

Abstract

:
In this paper, we investigate the solution to a class of symmetric non-homogeneous two-dimensional fractional integro-differential equations using both analytical and numerical methods. We first show the differences between the Caputo derivative and the symmetric sequential fractional derivative and how they help facilitate the implementation of numerical and analytical approaches. Then, we propose a numerical approach based on the operational matrix method, which involves deriving operational matrices for the differential and integral terms of the equation and combining them to generate a single algebraic system. This method allows for the efficient and accurate approximation of the solution without the need for projection. Our findings demonstrate the effectiveness of the operational matrix method for solving non-homogeneous fractional integro-differential equations. We then provide examples to test our numerical method. The results demonstrate the accuracy and efficiency of the approach, with the graph of exact and approximate solutions showing almost complete overlap, and the approximate solution to the fractional problem converges to the solution of the integer problem as the order of the fractional derivative approaches one. We use various methods to measure the error in the approximation, such as absolute and L 2 errors. Additionally, we explore the effect of the derivative order. The results show that the absolute error is on the order of 10 14 , while the L 2 error is on the order of 10 13 . Next, we apply the Laplace transform to find an analytical solution to a class of fractional integro-differential equations and extend the approach to the two-dimensional case. We consider all homogeneous cases. Through our examples, we achieve two purposes. First, we show how the obtained results are implemented, especially the exact solution for some 1D and 2D classes. We then demonstrate that the exact fractional solution converges to the exact solution of the ordinary derivative as the order of the fractional derivative approaches one.

1. Introduction

Fractional calculus is a branch of mathematical analysis that extends the concepts of differentiation and integration to non-integer orders. The origins of fractional calculus can be traced back to the work of Leibniz and Euler in the 18th century, who introduced the concept of fractional differentiation and integration. However, in the late 19th and early 20th centuries, fractional calculus began to be studied systematically, thanks to the work of mathematicians such as Liouville, Riemann, and Grunwald; see [1]. One of the main results of fractional calculus is the fractional derivative, which is a generalization of the classical derivative to non-integer orders. The fractional derivative has numerous applications in physics, engineering, and other sciences, particularly in problems involving non-locality or memory effects [2].
Another important concept in fractional calculus is the fractional integral, which is a generalization of the classical integral to non-integer orders. The fractional integral is useful for solving differential equations involving fractional derivatives and has applications in areas such as signal processing and control theory.
One of the most significant advancements in the field of fractional calculus is the development of the Caputo fractional derivative. Caputo introduced the fractional derivative in 1967, which is defined as the Riemann–Liouville fractional derivative with the initial conditions being in the form of a non-singular integer. The Caputo fractional derivative is widely used in various applications due to its ability to handle initial conditions better than other types of fractional derivatives; see [2] for more details.
The concept of symmetric sequential fractional derivatives was first introduced by Kilbas, Srivastava, and Trujillo in [3]. They defined the sequential fractional derivative of order α and β , α , β > 0 as the composition of two fractional derivatives of orders α and β , respectively. Later, a generalization of this concept was proposed by Zhang and Wei [4], who defined the sequential fractional derivative of arbitrary order as the composition of n fractional derivatives of different orders. In [5], Atanackovic and Pilipovic studied the symmetric sequential fractional derivative in the context of the Caputo sense. They derived the Laplace transform of the symmetric sequential fractional derivative and discussed its basic properties. They also investigated the relationship between the sequential fractional derivative and the fractional integral.
Li et al. [6] introduced a new class of symmetric sequential fractional derivative models and demonstrated their applications in the fields of fluid mechanics, heat transfer, and finance. Zhong et al. [7] discussed the properties of sequential fractional derivatives and their application to fractional partial differential equations. They also showed that the sequential fractional derivative is more general than the Riemann–Liouville fractional derivative.
One of the key properties of symmetric sequential fractional derivatives is their non-commutativity. This means that the order in which the fractional derivatives are applied can affect the final result. This property has been studied in detail by several researchers, such as Baleanu, Gülsuand Mohammadi [8], who have shown that it can have important implications for the behavior of complex systems. Symmetric sequential fractional derivatives have also been used in numerical methods for solving differential equations. Fractional Adams–Bashforth methods that use sequential fractional derivatives to approximate the solution of fractional differential equations have been proposed by Zhang, Li, and Shen [9].
Moreover, different applications in science and engineering are modeled by integro-differential equations (IDEs) and partial differential equations (PDEs); see [10,11]. They arise naturally in many fields, including physics, engineering, economics, and biology. One of the earliest studies of IDEs was performed by Volterra in [12], who developed a theory of integral equations based on his work on functional analysis. Then, the study of IDEs has grown significantly, with many researchers contributing to the development of the field.
On the other hand, in [2], the authors showed one of the key properties of IDEs is their non-locality, which means that the solution at a given point depends not only on the local values of the function but also on its integral over a certain range. This property makes the analysis of IDEs more challenging than that of ordinary differential equations (ODEs).
Several techniques have been developed for solving IDEs, including numerical methods, Laplace transforms, and Fourier transforms. One of the most popular methods for solving IDEs, which has been used by Baleanu et al. [13], is the method of characteristics, which involves finding curves along which the solution to the equation is constant. This method has been applied to a wide range of problems, including those in fluid mechanics, heat transfer, and finance.
The stability analysis of IDEs is another important topic in the study of these equations. Stability analysis involves studying the behavior of the solutions to the IDEs as time passes and determining whether the solutions converge or diverge. Several researchers have developed stability criteria for IDEs, which have been used in the study of problems in biology, finance, and other fields. Wang, Sun, and Li [14] studied the stability of solutions to fractional differential equations with the Caputo derivative. They presented a new sufficient condition for the asymptotic stability of the zero solution to the equation, based on the Lyapunov direct method. In recent years, there has been a growing interest in the use of operational matrices in machine learning and data analysis. Operational matrices can be used to represent data in a low-dimensional space, which can make it easier to analyze and process large datasets. Operational matrices are an important tool in applied mathematics and engineering for solving differential equations and other problems. They are matrices that used to calculate the estimated solution of a differential equation in terms of the values of the function and its derivatives at a finite set of points. Wu and Liao [15] presented an efficient approach to obtain higher-order approximations of fractional derivatives using the generalized Taylor matrix method.
Canuto et al. [16] provided one of the most common methods for solving differential equations using the spectral method. This method involves approximating the solution of a differential equation using a series of orthogonal functions, such as Fourier series or Chebyshev polynomials and then using the corresponding operational matrices to solve for the coefficients of the series. For more detail, we refer the reader to the references [17,18,19,20].
This article is divided into six sections. The first two sections provide a brief overview of the literature on the applications of and solutions to fractional IDEs. We also mention some preliminary definitions and theorems used in this work. More attention is paid to some symmetric sequential derivatives and fractional Laplace transform results. A derivation of the operational matrix method (OMM) can be found in Section 3. A proof of these matrices is given. Section 4 derives exact solutions for 1D and 2D fractional IDEs. All homogeneous cases are considered. Finally, the last two sections provide some examples for two purposes. First, we show how the results obtained are implemented, specifically the exact solution for some IDE classes. Next, we show the efficiency of OMM. In addition, we end the article with some conclusions and concluding remarks.
It is worth mentioning that our investigation focuses on solving a class of symmetric non-homogeneous two-dimensional fractional integro-differential equations using both analytical and numerical methods. Initially, we highlight the distinctions between the Caputo derivative and the symmetric sequential fractional derivative and discuss how these differences aid in the implementation of numerical and analytical approaches. Subsequently, we propose a numerical approach based on the operational matrix method. This approach entails deriving operational matrices for the differential and integral terms of the equation and combining them to form a unified algebraic system. This method enables efficient and precise approximation of the solution without requiring projection.

2. Preliminaries

We start this section by defining two important operators, which are the Caputo fractional derivative and the Riemann–Liouville fractional integral operator.
Definition 1
([3,21]). Let q , s > 0 and M N with M 1 < q < M . Then, the Caputo derivative of z(s) of order q is defined as
c D q z ( s ) = 1 Γ ( M q ) 0 s ( s r ) M q 1 z ( M ) ( r ) d r
and the fractional integral operator is given by
I q z ( s ) = 1 Γ ( q ) 0 s ( s r ) q 1 z ( r ) d r .
The power rule is one of the important rules, and it is given as follows for both operators
c D q s γ = 0 , γ < q , γ 0 , 1 , 2 , . . . Γ ( γ + 1 ) Γ ( γ q + 1 ) z γ q otherwise
and
I q s γ = Γ ( γ + 1 ) Γ ( γ + q + 1 ) z γ + q .
The relations between these operators are given as follows:
c D q I q z ( s ) = z ( s )
and
I q c D q z ( x ) = z ( x ) j = 0 M 1 z j ( 0 ) j ! s j .
Formula (3) implies that if 1 2 < q 1 , then
c D 2 q s q = Γ ( q + 1 ) Γ ( q 2 q + 1 ) s q 2 q = Γ ( q + 1 ) Γ ( 1 q ) s q
and
c D q c D q s q = c D q Γ ( q + 1 ) Γ ( 1 ) = 0 .
Thus,
c D 2 q s q c D q c D q s q .
Hence, the Caputo derivative is not symmetric-sequential. The definition of a symmetric-sequential derivative is given as follows.
Definition 2
([22,23]). If q ( 0 , 1 ) with
c D M q z ( s ) = c D q c D ( M 1 ) q z ( s ) , M N ,
then the derivative is called a symmetric-sequential Caputo derivative of order p, which we denote by  s c D M q z ( s ) . For simplicity, we use  D M q z ( s )  to mean symmetric sequential derivative.
Now, we define a list of interesting functions that we use in this paper.
Definition 3.
Let α 1 > 0 , α 2 > 0 , and λ , γ . Then,
1. 
The one and two Mittag–Leffler functions are
E α 1 ( s ) = j = 0 s j Γ ( j α 1 + 1 ) , E α 1 , α 2 ( s ) = j = 0 s j Γ ( j α 1 + α 2 ) ,
respectively.
2. 
The fractional sine and cosine functions are
s i n α 1 , λ , γ ( s ) = E α 1 ( ( λ + i γ ) s ) E α 1 ( ( λ i γ ) s ) 2 i , c o s α 1 , λ , γ ( s ) = E α 1 ( ( λ + i γ ) s ) + E α 1 ( ( λ i γ ) s ) 2 ,
respectively.
For the sequential fractional derivative, we mention two important results in the following theorem.
Theorem 1.
Let q > 0 and L ( r ( t ) ) = R ( s ) . Then, the Laplace transforms D q r ( t ) and D 2 q r ( t ) on [ 0 , ) are given as
L ( D q r ( t ) ) = s q R ( s ) s q 1 r ( 0 )
and
L ( D 2 q r ( t ) ) = s 2 q R ( s ) s 2 q 1 r ( 0 ) s q 1 D q r ( 0 ) .
In the next theorem, we list the Laplace transforms of some functions that we use in this paper.
Theorem 2.
Let q > 0 and α . Then, the Laplace transform of Mittag–Leffler and fractional trigonometric functions are given as follows:
L ( E q ( ± α t q ) ) = s q 1 s q α , s q > α ,
L ( t q E q , q ( ± α t q ) ) = q s q 1 ( s q α ) 2 , s q > α ,
L ( sin q , λ , γ ( t q ) ) = μ s q 1 ( s q λ ) 2 + μ 2 ,
L ( cos q , λ , γ ( t q ) ) = s q 1 ( s q λ ) ( s q λ ) 2 + μ 2 ,
L ( t 2 q 1 G E q ( α t q ) ) = 1 ( s q α ) 2 ,
where
G E q ( t ) = j = 0 ( j + 1 ) t j Γ ( q ( j + 2 ) ) .

3. Non-Homogeneous Two-Dimensional Fractional Integro-Differential Equations

This section is devoted to the method of solving a class of two-dimensional fractional inegro-differential Equations of the form
w x x ( x , t ) + D t 2 q w ( x , t ) + ν 1 D t q w ( x , t ) + w x ( x , t ) + ν 2 ( t ) w ( x , t ) + ν 3 I t q w ( x , t ) = g ( x , t )
with
w ( 0 , t ) = w ( 1 , t ) = 0 , t 0
and
w ( x , 0 ) = r ( x ) , D t q w ( x , 0 ) = 0 , 0 < x < 1
on the domain [ 0 , 1 ] × [ 0 , η ] , where ν 1 , ν 3 are real constants, r C 3 ( [ 0 , 1 ] ) , ν 2 C [ 0 , η ] and g C ( [ 0 , 1 ] × [ 0 , η ) . Before we start the method of solution, we need the following definition.
Definition 4
([24,25,26]). Let M 1 , M 2 N and Δ x = 1 M 1 and Δ t = η M 2 be two step sizes in the x and t directions, respectively. Let x i = i Δ x and t j = j Δ t for i A x = { 0 , 1 , , M 1 1 } and j A t = { 0 , 1 , , M 2 1 } . The two-dimensional Block Pulse function (BPF) is a function β i , j : [ 0 , 1 ) × [ 0 , η ) defined by
β i , j ( x , t ) = 1 , x [ x i , x i + 1 ) , t [ t j , t j + 1 ) 0 , otherwise
where i A x and j A t .
The two-dimensional BPF can be split into the product of two one-dimensional BPFs as follows.
Theorem 3.
For any i A x and j A t , we have
β i , j ( x , t ) = β x i ( x ) β t j ( t )
where β x i ( x ) and β t j ( t ) are BPFs on [ 0 , η ) and [ 0 , T m a x ) , respectively.
Proof. 
The proof follows directly from the fact that ( x , t ) [ a , b ) × [ c , d ) if and only if x [ a , b ) and t [ c , d ) . □
Since [ 0 , 1 ) × [ t , η ) is divided into disjoint sets { [ x i , x i + 1 ) × [ t j , t j + 1 ) : i A x , j A t } , one can see the following product and orthogonality relations
β i 1 , j 1 ( x , t ) β i 2 , j 2 ( x , t ) = β i 1 , j 1 ( x , t ) , i 1 = i 2 , j 1 = j 2 0 , otherwise ,
0 1 0 η β i 1 , j 1 ( x , t ) β i 2 , j 2 ( x , t ) d t d x = Δ x Δ t , i 1 = i 2 , j 1 = j 2 0 , otherwise .
Using these two properties, we can prove the completeness property.
Theorem 4.
Let w L 2 ( [ 0 , 1 ) × [ 0 , η ) ) be a bounded function. Then,
w ( x , t ) i = 0 M 1 1 j = 0 M 2 1 w i j β i , j ( x , t )
with
w i j = 1 Δ x Δ t x i x i + 1 t j t j + 1 w ( x , t ) d t d x .
Proof. 
After multiplying Equation (28) by β i , j ( x , t ) and then integrating both sides, the result of the theorem will follow directly. □
Using Theorem 4, we can approximate the functions w ( x , t ) , ν 2 ( t ) , g ( x , t ) , and r ( x ) as
w ( x , t ) i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) β t j ( t ) ,
g ( x , t ) i = 0 M 1 1 j = 0 M 2 1 g i j β x i ( x ) β t j ( t ) ,
ν 2 ( t ) j = 0 M 2 1 v 2 j β t j ( t ) ,
r ( x ) i = 0 M 1 1 r i β i x ( x ) .
Hence, the sequential derivatives of the function w can be expanded as
D t 2 q w ( x , t ) = i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) D t 2 q β t j ( t ) ,
D t q w ( x , t ) = i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) D t q β t j ( t ) .
In addition, the product and orthogonality relations of BPF imply that
ν 2 ( t ) w ( x , t ) = j = 0 M 2 1 v 2 j β t j ( t ) i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) β t j ( t ) = i = 0 M 1 1 j = 0 M 2 1 v 2 j w i j β x i ( x ) β t j ( t ) .
If we substitute these approximations in Equation (21), we obtain
i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) D t 2 q β t j ( t ) + ν 1 i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) D t q β t j ( t ) + i = 0 M 1 1 j = 0 M 2 1 v 2 j w i j β x i ( x ) β t j ( t ) + ν 3 i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) I t q β t j ( t ) = i = 0 M 1 1 j = 0 M 2 1 g i , j β x i ( x ) β t j ( t ) w x x ( x , t ) w x ( x , t ) .
To be able to write Equation (37) in matrix form, let us define the following matrices:
β x = β x 0 ( x ) β x 1 ( x ) β x ( M 1 1 ) ( x ) , W = w 00 w 01 w 0 ( M 2 1 ) w 10 w 11 w 1 ( M 2 1 ) w ( M 1 1 ) 0 w ( M 1 1 ) 1 w ( M 1 1 ) ( M 2 1 ) ,
β t = β t 0 ( t ) β t 1 ( t ) β t ( M 2 1 ) ( t ) , V 2 = ν 20 0 0 0 ν 21 0 0 0 ν 2 ( M 2 1 ) , R = r 0 r 1 r M 1 1 ,
H = g 00 g 01 g 0 ( M 2 1 ) g 10 g 11 g 1 ( M 2 1 ) g ( M 1 1 ) 0 g ( M 1 1 ) 1 g ( M 1 1 ) ( M 2 1 ) .
Thus, Equation (37) becomes
β x * ( x ) W D t 2 q β t ( t ) + ν 1 β x * ( x ) W D t q β t ( t ) + β x * ( x ) W V 2 β t ( t ) + ν 3 β x * ( x ) W I t q β t ( t ) = β x * ( x ) G β t ( t ) w x x ( x , t ) w x ( x , t )
where * mean the transpose of a matrix. The initial conditions yield to
D t q w ( x , 0 ) = i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) D t q β t j ( 0 ) = β x ( x ) * W D t q β t ( 0 ) = 0 ,
w ( x , 0 ) = i = 0 M 1 1 j = 0 M 2 1 w i j β x i ( x ) β t j ( 0 ) = β x * ( x ) W β t ( 0 ) = β x * R ,
which give
W D t q β t ( 0 ) = 0
and
W β t ( 0 ) = R .
The next step is to find the operational integration matrix of the Riemann integral.
Theorem 5.
Let
I ( r ( x ) ) = 0 x r ( t ) d t .
Then, the operational matrix of I is
O I = Δ x 2 1 2 2 2 2 0 1 2 2 2 0 0 1 2 2 0 0 0 1 2 0 0 0 0 1 .
Proof. 
Let 0 k M 1 1 and k N . Then,
I ( β x k ( x ) ) = 0 x β x k ( t ) d t
= 0 , x < x k x x k , x k x < x k + 1 Δ x , x k + 1 x < 1 .
Therefore, if
I ( β x k ( x ) ) = j = 0 M 1 1 γ j , k β x j ( x ) ,
then
γ j , k = 1 Δ x 0 1 I ( β x k ( t ) β x j ( t ) d t
= 1 Δ x x j x j + 1 I ( β x k ( t ) d t
= Δ x 2 , 0 j = k M 1 1 Δ x , 0 j < k M 1 1 0 , 0 k < j M 1 1
which completes the proof. □
Theorem 5 and Equation (41) give the following equation when we take the Riemann integral
β x * ( x ) O I * W D t 2 q β t ( t ) + ν 1 β x * ( x ) O I * W D t q β t ( t ) + β x * ( x ) O I * W V 2 β t ( t ) + ν 3 β x * ( x ) O I * W I t q β t ( t ) = β x * ( x ) O I * G β t ( t ) w x ( x , t ) + w x ( 0 , t ) w ( x , t ) + w ( x , 0 ) .
Assume that w x ( 0 , t ) = λ ( t ) . Then,
λ ( t ) = k = 0 M 2 1 λ k β t k ( t ) = i = 0 M 1 1 k = 0 M 2 1 λ k β x i ( x ) β t k ( t ) = β x * ( x ) Λ β t ( t )
where
Λ = λ 0 λ 1 λ M 2 1 λ 0 λ 1 λ M 2 1 λ 0 λ 1 λ M 2 1 .
This is true since
k = 0 M 1 1 β x k ( x ) = 1 .
Hence, Equation (54) and the boundary condition w ( x , 0 ) = 0 give
β x * ( x ) O I * W D t 2 q β t ( t ) + ν 1 β x * ( x ) O I * W D t q β t ( t ) + β x * ( x ) O I * W V 2 β t ( t ) + ν 3 β x * ( x ) O I * W I t q β t ( t ) = β x * ( x ) O I * G β t ( t ) w x ( x , t ) + β x * ( x ) Λ β t ( t ) β x * ( x ) W β t ( t ) .
Taking the Riemann integral for both sides of Equation (58), we obtain
β x * ( x ) O I * O I * W D t 2 q β t ( t ) + ν 1 β x * ( x ) O I * O I * W D t q β t ( t ) + β x * ( x ) O I * O I * W V 2 β t ( t ) + ν 3 β x * ( x ) O I * O I * W I t q β t ( t ) = β x * ( x ) O I * O I * G β t ( t ) β x * ( x ) W β t ( t ) + β x * ( x ) O I * Λ β t ( t ) β x * ( x ) O I * W β t ( t ) .
Since the BPFs { β x k : k = 0 , 1 , . . . , M 1 1 } are linearly independent,
O I * O I * W D t 2 q β t ( t ) + ν 1 O I * O I * W D t q β t ( t ) + O I * O I * W V 2 β t ( t ) + ν 3 O I * O I * W I t q β t ( t ) = O I * O I * G β t ( t ) W β t ( t ) + O I * Λ β t ( t ) O I * W β t ( t ) .
In the next theorem, we find the operational matrix of I t q .
Theorem 6.
The operational matrix of I t q is
O q = Δ t q Γ ( q + 2 ) 1 s 1 s 2 s M 2 2 s M 2 1 0 1 s 1 s M 2 3 s M 2 2 0 0 1 s M 2 4 s M 2 3 0 0 0 1 s 1 0 0 0 0 1
where s μ = ( μ + 1 ) q + 1 2 μ q + 1 + ( μ 1 ) q + 1 , μ = 1 , 2 , . . . , M 2 1 .
Proof. 
Let 0 μ M 2 1 and μ N . Then,
I t q β t μ ( t ) = 1 Γ ( q ) 0 t ( t z ) q 1 β t μ ( z ) d z
= 0 , t < t μ ( t t μ ) q Γ ( q + 1 ) , t μ t < t μ + 1 ( t t μ ) q ( t t μ + 1 ) q Γ ( q + 1 ) , t μ + 1 t < η .
If
I t q β t μ ( t ) = k = 0 M 2 1 ξ k , μ β t k ( t ) ,
then
ξ k , μ = 1 Δ t 0 η I t q β t μ ( v ) β t k ( v ) d v
= 1 Δ t t k t k + 1 I t q β t μ ( v ) d v
= Δ t q Γ ( q + 2 ) , 0 k = μ M 2 1 Δ t q ( μ k + 1 ) q + 1 2 ( μ k ) q + 1 + ( μ k 1 ) q + 1 Γ ( q + 2 ) , 0 k < μ M 2 1 0 , 0 μ < k M 2 1 .
If we assume that s μ = ( μ + 1 ) q + 1 2 μ q + 1 + ( μ 1 ) q + 1 , and μ = 1 , 2 , . . . , M 2 1 , and the proof is completed. □
One can see that
R = R k = 0 M 2 β t k ( t ) = R O R β t ( t )
where O R = 1 1 1 is 1 × M 2 matrix since k = 0 M 2 1 β t k ( t ) = 1 . Theorem 6 and Equations (60) and (68) yield
O I * O I * W D t q β t ( t ) + ν 1 O I * O I * W β t ( t ) ν 1 O I * O I * R O R β t ( t ) + O I * O I * W V 2 O q β t ( t ) + ν 3 O I * O I * W O q O q β t ( t ) = O I * O I * G O q β t ( t ) W O q β t ( t ) + O I * Λ O q β t ( t ) O I * W O q β t ( t ) .
Apply Theorem 6 one more time on Equation (69) to obtain
O I * O I * W β t ( t ) O I * O I * R O R O q β t ( t ) + ν 1 O I * O I * W O q β t ( t ) ν 1 O I * O I * R O R O q β t ( t ) + O I * O I * W V 2 O q O q β t ( t ) + ν 3 O I * O I * W O q O q O q β t ( t ) = O I * O I * G O q O q β t ( t ) W O q O q β t ( t ) + O I * Λ O q O q β t ( t ) O I * W O q O q β t ( t ) .
Since { β t j } 0 M 2 1 are linearly independent, then
O I * O I * W O I * O I * R O R O q + ν 1 O I * O I * W O q ν 1 O I * O I * R O R O q + O I * O I * W V 2 O q O q ( t ) + ν 3 O I * O I * W O q O q O q = O I * O I * G O q O q W O q O q + O I * Λ O q O q O I * W O q O q .
Since
w ( 1 , t ) = i = 0 M 1 1 k = 0 M 2 1 w i , k β x i ( 1 ) β t k ( t ) = 0 ,
then
β x ( 1 ) * V = 0 .
Hence, we solve the algebraic system
O I * O I * W O I * O I * R O R O q + ν 1 O I * O I * W O q ν 1 O I * O I * R O R O q + O I * O I * W V 2 O q O q ( t ) + ν 3 O I * O I * W O q O q O q = O I * O I * G O q O q W O q O q + O I * Λ O q O q O I * W O q O q
β x ( 1 ) * V = 0
using Mathematica for the unknowns, which are the matrix W and the vector
λ 0 λ 1 λ M 2 .

4. Analytical Solution of a Class of Two-Dimensional Fractional Integro-Differential Equations

This section is devoted to studying the solution of a class of two-dimensional fractional IDEs analytically. First, The Laplace transform is used to find an analytical solution for the following class of fractional IDEs of the form
D 2 q r ( t ) + α 1 D q r ( t ) + α 2 r ( t ) + α 3 I q r ( t ) = 0
with
r ( 0 ) = β 1 , D q r ( 0 ) = β 2
where α 1 , α 2 , α 3 , β 1 , β 2 are constants, and 1 2 < q 1 .
Theorem 7.
The solution to the problem (77) and (78) is given as follows.
1. 
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has one real root λ and two complex roots μ ± i ν , then the solution is
r ( t ) = c 1 E q ( λ t q ) + ( c 2 + μ c 3 ) s i n q ( ( μ + i ν ) t q ) + c 3 c o s q ( ( μ + i ν ) t q ) .
2. 
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three distinct real roots λ 1 , λ 2 , and λ 3 , then
r ( t ) = c 1 E q ( λ 1 t q ) + c 2 E q ( λ 2 t q ) + c 3 E q ( λ 3 t q ) .
3. 
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three real roots λ 1 = λ 2 = λ , and λ 3 λ , then
r ( t ) = c 1 E q ( λ t q ) + c 2 t q q E q , q ( λ t q ) + c 3 E q ( λ 3 t q ) .
4. 
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three real roots λ 1 = λ 2 = λ 3 = λ , then
r ( t ) = c 1 E q ( λ t q ) + c 2 t q q E q , q ( λ t q ) + c 3 E q ( λ t q ) ( t 2 q 1 G E q ( λ t q ) )
where * means the convolution.
Proof. 
We take the Laplace transform for both sides of Equation (77) to obtain
s 2 q R ( s ) s 2 q 1 r ( 0 ) s q 1 D q r ( 0 ) + α 1 s q R ( s ) s q 1 u ( 0 ) + α 2 R ( s ) + α 3 s q R ( s ) = 0
which can be rewritten as
s 2 q R ( s ) β 1 s 2 q 1 β 2 s q 1 + α 1 s q R ( s ) β 1 s q 1 + α 2 R ( s ) + α 3 s q R ( s ) = 0 .
Simple calculation implies that
R ( s ) = β 1 s 3 q 1 + β 2 + α 1 β 1 s 2 q 1 s 3 q + α 1 s 2 q + α 2 s q + α 3
or
R ( s ) = s q 1 β 1 s 2 q + β 2 + α 1 β 1 s q s 3 q + α 1 s 2 q + α 2 s q + α 3 .
Let z = s q . Then,
β 1 s 2 q + β 2 + α 1 β 1 s q s 3 q + α 1 s 2 q + α 2 s q + α 3 = β 1 z 2 + β 2 + α 1 β 1 z z 3 + α 1 z 2 + α 2 z + α 3 .
Now, we have four cases to consider.
1.
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has one real root λ and two complex roots μ ± i ν , then simple calculations yield to
R ( s ) = c 1 s q 1 s q λ + s q 1 ( c 2 + c 3 s q ) s q μ i ν s q μ + i ν = c 1 s q 1 s q λ + s q 1 ( c 2 + c 3 s q ) s q μ 2 + ν 2 = c 1 s q 1 s q λ + ( c 2 + c 3 μ ) s q 1 s q μ 2 + ν 2 + c 3 s q 1 ( s q μ ) s q μ 2 + ν 2
where
c 1 = λ ( β 2 + ( α 1 + λ ) β 1 ) ( μ λ ) 2 + ν 2 , c 2 = ( μ 2 + ν 2 ) ( β 2 + ( α 1 + λ ) β 1 ) ( μ λ ) 2 + ν 2 ,
c 3 = β 1 λ ( β 2 + ( α 1 + λ ) β 1 ) ( μ λ ) 2 + ν 2 .
Using Theorem 2, we obtain
r ( t ) = c 1 E q ( λ t q ) + ( c 2 + μ c 3 ) s i n q ( ( μ + i ν ) t q ) + c 3 c o s q ( ( μ + i ν ) t q ) .
2.
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three distinct real roots λ 1 , λ 2 , and λ 3 , then
R ( s ) = c 1 s q 1 s q λ 1 + c 2 s q 1 s q λ 2 + c 3 s q 1 s q λ 3
where
c 1 = λ 1 ( α 1 β 1 + β 2 + λ 1 β 1 ) ( λ 2 λ 1 ) ( λ 1 λ 3 ) , c 2 = λ 2 ( α 1 β 1 + β 2 + λ 2 β 1 ) ( λ 2 λ 1 ) ( λ 2 λ 3 )
c 3 = α 1 β 1 λ 3 + β 2 λ 3 + β 1 λ 3 2 ) ( λ 1 λ 3 ) ( λ 2 λ 3 ) .
Using Theorem 2, we obtain
r ( t ) = c 1 E q ( λ 1 t q ) + c 2 E q ( λ 2 t q ) + c 3 E q ( λ 3 t q ) .
3.
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three real roots λ 1 = λ 2 = λ , and λ 3 λ , then
R ( s ) = c 1 s q 1 s q λ + c 2 s q 1 ( s q λ ) 2 + c 3 s q 1 s q λ 3
where
c 1 = β 1 α 1 β 1 λ 3 + β 2 λ 3 + β 1 λ 3 2 ( λ λ 3 ) 2 , c 2 = α 1 β 1 λ + β 2 λ + β 1 λ 2 λ λ 3
c 3 = α 1 β 1 λ 3 + β 2 λ 3 + β 1 λ 3 2 ( λ λ 3 ) 2 .
Using Theorem 2, we obtain
r ( t ) = c 1 E q ( λ t q ) + c 2 t q q E q , q ( λ t q ) + c 3 E q ( λ 3 t q ) .
4.
If z 3 + α 1 z 2 + α 2 z + α 3 = 0 has three real roots λ 1 = λ 2 = λ 3 = λ , then
R ( s ) = c 1 s q 1 s q λ + c 2 s q 1 ( s q λ ) 2 + c 3 s q 1 s q λ 1 ( s q λ ) 2 .
Using Theorem 2, we obtain
r ( t ) = c 1 E q ( λ t q ) + c 2 t q q E q , q ( λ t q ) + c 3 E q ( λ t q ) ( t 2 q 1 G E q ( λ t q ) )
where
c 1 = β 1 , c 2 = ( α 1 + 2 λ ) β 1 + β 2 , c 3 = ( α 1 λ + λ 2 ) β 1 + β 2 λ .
Now, we generalize the idea to study the two-dimensional fractional IDEs. The solution to this class is given in the following theorem.
Theorem 8.
Let ν 1 , ν 2 , and ν 3 be constants. Then, the solution to the following fractional IDE
w x x ( x , t ) + D t 2 q w ( x , t ) + ν 1 D t q w ( x , t ) + w x ( x , t ) + ν 2 w ( x , t ) + ν 3 I t q w ( x , t ) = 0
with
w ( 0 , t ) = w ( 1 , t ) = 0 , t 0
is given by
w ( x , t ) = n = 1 z n ( x ) r n ( t ) = n = 1 e 1 2 x r n ( t ) sin ( n π x )
where { r n ( t ) } are given as follows
1. 
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has one real root ξ n and two complex roots μ n ± i ν n , then
r n ( t ) = a n E q ( ξ n t q ) + ( b n + μ n c n ) s i n q ( ( μ n + i ν n ) t q ) + c n c o s q ( ( μ n + i ν n ) t q ) .
2. 
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three distinct real roots ξ 1 n , ξ 2 n , and ξ 3 n , then
r n ( t ) = a n E q ( ξ 1 n t q ) + b n E q ( ξ 2 n t q ) + c n E q ( ξ 3 n t q ) .
3. 
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three real roots ξ 1 n = ξ 2 n = ξ n , and ξ 3 n ξ n , then
r n ( t ) = a n E q ( ξ n t q ) + b n t q q E q , q ( ξ n t q ) + c n E q ( ξ 3 n t q ) .
4. 
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three real roots ξ 1 n = ξ 2 n = ξ 3 n = ξ n , then
r n ( t ) = a n E q ( ξ n t q ) + b n t q q E q , q ( ξ n t q ) + c n E q ( ξ n t q ) ( t 2 q 1 G E q ( ξ n t q ) ) .
Proof. 
Following the idea of the separation of variables, we assume that
w ( x , t ) = z ( x ) r ( t ) .
Then, substitute Equation (110) into Equation (103) and divide both sides by z ( x ) r ( t ) to obtain
z ( x ) + z ( x ) z ( x ) = D t 2 q r ( t ) + ν 1 D t q r ( t ) + ν 2 r ( t ) + ν 3 I t q r ( t ) r ( t ) = λ
where λ is a constant. The boundary conditions in Equation (104) yield to
z ( 0 ) = z ( 1 ) = 0 .
Hence, we obtain the following eigenvalue problem
z ( x ) + z ( x ) z ( x ) = λ , z ( 0 ) = z ( 1 ) = 0 .
Thus, the auxiliary equation is
α 2 + α λ = 0
which implies that
α = 1 ± 1 + 4 λ 2 .
Three cases should be considered when λ is either equal to, less than, or greater than 1 4 . These three cases give the following eigenvalues and eigenfunctions
z n ( x ) = e 1 2 x sin ( n π x ) , λ n = 4 n 2 π 2 + 1 4 , n N .
Hence, Equation (111) gives
D t 2 q r ( t ) + ν 1 D t q r ( t ) + ν 2 4 n 2 π 2 + 1 4 r ( t ) + ν 3 I t q r ( t ) = 0
Using Theorem 7, we obtain
1.
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has one real root ξ n and two complex roots μ n ± i ν n , then
r n ( t ) = a n E q ( ξ n t q ) + ( b n + μ n c n ) s i n q ( ( μ n + i ν n ) t q ) + c n c o s q ( ( μ n + i ν n ) t q ) .
2.
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three distinct real roots ξ 1 n , ξ 2 n , and ξ 3 n , then
r n ( t ) = a n E q ( ξ 1 n t q ) + b n E q ( ξ 2 n t q ) + c n E q ( ξ 3 n t q ) .
3.
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three real roots ξ 1 n = ξ 2 n = ξ n , and ξ 3 n ξ n , then
r n ( t ) = a n E q ( ξ n t q ) + b n t q q E q , q ( ξ n t q ) + c n E q ( ξ 3 n t q ) .
4.
If α 3 + ν 1 α 2 + ν 2 4 n 2 π 2 + 1 4 α + ν 3 = 0 has three real roots ξ 1 n = ξ 2 n = ξ 3 n = ξ n , then
r n ( t ) = a n E q ( ξ n t q ) + b n t q q E q , q ( ξ n t q ) + c n E q ( ξ n t q ) ( t 2 q 1 G E q ( ξ n t q ) ) .
Therefore, Equations (116) and (118)–(121) give
w ( x , t ) = n = 1 z n ( x ) r n ( t ) = n = 1 e 1 2 x r n ( t ) sin ( n π x )
where { r n ( t ) } n = 1 are given in Equations (118)–(121). □
If we have initial conditions, then we use the Fourier series expansions to find the { a n } , { b n } , and { c n } .

5. Illustrative Examples

In this section, we discuss several examples to explain the previous discussion. The first example will be about Theorem 7. We find the exact solution to the one-dimensional fractional IDEs with constant coefficients. Then, we will discuss the behavior of the solution when q approaches one.
Example 1.
Consider the following one-dimensional fractional IDE of the form
D 2 q r ( t ) + α 1 D q r ( t ) + α 2 r ( t ) + α 3 I q r ( t ) = 0
with
r ( 0 ) = 1 , D q r ( 0 ) = 1 .
We will consider four cases:
1. 
If α 1 = 3 , α 2 = 4 , and α 3 = 2 , then the roots of z 3 + α 1 z 2 + α 2 z + α 3 = 0 are 1 , 1 ± i . Using the first part of Theorem 7, c 1 = 1 , c 2 = 2 and c 3 = 2 . Then, the solution of Equations (123) and (124) will be
r q ( t ) = E q ( t q ) + 2 cos q ( ( 1 + i ) t q ) .
Then,
lim q 1 r q ( t ) = lim q 1 ( E q ( t q ) + 2 cos q ( ( 1 + i ) t q ) )
= e t + 2 e t c o s t ,
which is the solution of
r ( t ) 3 r ( t ) + 4 r ( t ) 2 0 t r ( s ) d s = 0 , r ( 0 ) = r ( 0 ) = 1 .
Note that
lim q 1 E q ( λ t q ) = e λ t , lim q 1 cos q ( ( μ + i ν ) t q ) ) = e μ t cos ( ν t ) ,
lim q 1 sin q ( ( μ + i ν ) t q ) ) = e μ t sin ( ν t ) .
2. 
If α 1 = 6 , α 2 = 11 , and α 3 = 6 , then the roots of z 3 + α 1 z 2 + α 2 z + α 3 = 0 are 1 , 2 , 3 . Using the second part of Theorem 7, we obtain the result that c 1 = 2 , c 2 = 6 , and c 3 = 3 . Then, the solution of Equations (123) and (124) will be
r q ( t ) = 2 E q ( t q ) + 6 E q ( 2 t q ) 3 E q ( 3 t q ) .
Then,
lim q 1 r q ( t ) = lim q 1 ( 2 E q ( t q ) + 6 E q ( 2 t q ) 3 E q ( 3 t q ) )
= 2 e t + 6 e 2 t 3 e 3 t ,
which is the solution of
r ( t ) 6 r ( t ) + 11 r ( t ) 6 0 t r ( s ) d s = 0 , r ( 0 ) = r ( 0 ) = 1 .
3. 
If α 1 = 4 , α 2 = 5 , and α 3 = 2 , then the roots of z 3 + α 1 z 2 + α 2 z + α 3 = 0 are 1 , 1 , 2 . Using the third part of Theorem 7, c 1 = 3 , c 2 = 2 and c 3 = 2 . Then, the solution of Equations (123) and (124) will be
r q ( t ) = 3 E q ( t q ) + 2 t q q E q , q ( t q ) 2 E q ( 2 t q ) .
Then,
lim q 1 r q ( t ) = lim q 1 ( 3 E q ( t q ) + 2 t q q E q , q ( t q ) 2 E q ( 2 t q ) )
= 3 e t + 2 t e t 2 e 2 t ,
which is the solution of
r ( t ) 4 r ( t ) + 5 r ( t ) 2 0 t r ( s ) d s = 0 , r ( 0 ) = r ( 0 ) = 1 .
4. 
If α 1 = 3 , α 2 = 3 , and α 3 = 1 , then the roots of z 3 + α 1 z 2 + α 2 z + α 3 = 0 are 1 , 1 , 1 . Using the fourth part of Theorem 7, we find that c 1 = 1 , c 2 = 0 , and c 3 = 1 . Then, the solution of Equations (123) and (124) will be
r q ( t ) = E q ( t q ) E q ( t q ) ( t 2 q 1 G E q ( t q ) ) .
Then,
lim q 1 r q ( t ) = lim q 1 ( E q ( t q ) E q ( t q ) ( t 2 q 1 G E q ( t q ) )
= e t e t t e t = e t 1 2 t 2 e t
which is the solution of
r ( t ) 3 r ( t ) + 3 r ( t ) 0 t r ( s ) d s = 0 , r ( 0 ) = r ( 0 ) = 1 .
Note that
lim q 1 G E q ( t q ) ) = j = 0 ( j + 1 ) t j Γ ( j + 2 ) = j = 0 ( j + 1 ) t j ( j + 1 ) ! = e t .
Now, we will show how we can implement Theorem 8.
Example 2.
Consider the following two-dimensional fractional IDE of the form
w x x ( x , t ) + D t 2 q w ( x , t ) 3 D t q w ( x , t ) + w x ( x , t ) + 3 w ( x , t ) I t q w ( x , t ) = 0
with
w ( 0 , t ) = w ( 1 , t ) = 0 , t 0
and
w ( x , 0 ) = e 1 2 x sin ( π x ) , D q w ( x , 0 ) = e 1 2 x sin ( 2 π x ) , 0 x 1 .
Theorem 8 gives the result that the solution of Equations (144)–(146) has the following form
w ( x , t ) = n = 1 e 1 2 x r n ( t ) s i n ( n π x )
where r n ( t ) satisfies the equation
D t 2 q r n ( t ) 3 D t q r n ( t ) + 3 4 n 2 π 2 + 1 4 r n ( t ) I t q r n ( t ) = 0 .
The initial conditions in Equation (146) give
r n ( 0 ) = 1 , n = 1 0 , n > 1 , D q r n ( 0 ) = 1 , n = 2 0 , n = 1 o r n > 2 .
Simple calculations imply that the equation z 3 3 z 2 + 3 4 n 2 π 2 + 1 4 z 1 = 0 has one real root
λ n = γ n 12 2 3 12 2 3 γ n + 1 ,
and two complex roots
μ n ± i ν n = 1 + 6 2 3 γ n + γ n 24 2 3 ± i 6 2 3 3 γ n a 8 2 3 3
where
γ n = 1728 π 2 n 2 + 1728 π 2 n 2 432 2 11943936 432 3 .
for n > 2 , β 1 , n = r n ( 0 ) = 0 , and β 2 , n = D q r n ( 0 ) = 0 . Using the first part of Theorem 8, we have
r n ( t ) = 0 , n > 2 .
for n = 1 , β 1 , 1 = r 1 ( 0 ) = 1 , and β 2 , 1 = D q r n ( 0 ) = 0 . Using the first part of Theorem 8, we have
c 1 = 0.127535 , c 2 = 0.0790346 , c 3 = 0.872465 .
Thus,
r 1 ( t ) = 0.127535 E q ( 3.61971 t q ) 0.191303 sin q ( ( 0.309854 + 1.46532 i ) t q )
+ 0.872465 cos q ( ( 0.309854 + 1.46532 i ) t q ) .
For n = 2 , β 1 , 1 = r 1 ( 0 ) = 0 and β 2 , 1 = D q r n ( 0 ) = 1 . Using the first part of Theorem 8, we have
c 1 = 0.123247 , c 2 = 0.210087 , c 3 = 0.123247 .
Thus,
r 2 ( t ) = 0.123247 E q ( 4.7046 t q ) + 0.31513 sin q ( ( 0.852301 2.70057 i ) t q )
0.123247 cos q ( ( 0.852301 2.70057 i ) t q ) .
Therefore, the solution of the Problem (144)–(146) is given by
w q ( x , t ) = e 1 2 x sin ( π x ) r 1 ( t ) + e 1 2 x sin ( 2 π x ) r 2 ( t ) .
In the next two examples, we will test our numerical approach described in Section 3.
Example 3.
Consider the following two-dimensional fractional IDE of the form
w x x ( x , t ) + D t 7 / 4 w ( x , t ) 2 D t 7 / 8 w ( x , t ) + w x ( x , t ) + 3 w ( x , t ) 4 I t 7 / 8 w ( x , t ) = g ( x , t )
with
w ( 0 , t ) = w ( 1 , t ) = 0 , t 0
and
w ( x , 0 ) = e 1 2 x sin ( 3 π x ) , D q w ( x , 0 ) = 0 , 0 x 1
where
g ( x , t ) = 11 4 π 2 t 7 / 4 8 t 7 / 8 Γ 11 4 Γ 15 8 16 t 21 / 8 Γ 11 4 Γ 29 8 + 4 Γ 11 4 sin ( π x )
+ 1 4 e x 2 16 t 7 / 8 Γ 15 8 36 π 2 + 11 sin ( 3 π x ) .
Then, the exact solution is
w ( x , t ) = t 7 / 4 e x 2 sin ( π x ) + e x 2 sin ( 3 π x ) .
Here, we use M 1 = M 2 = 40 . Then, the absolute error is defined by
ϵ ( x , t ) = w ( x , t ) w 40 , 40 ( x , t ) .
The absolute errors for different values of x and t are reported in Table 1.
We also compute the L 2 -error, which is given by
ϵ L 2 = 0 1 0 1 w ( x , t ) w 40 , 40 ( x , t ) 2 d x d t = 2.84 10 13 .
In addition, the graphs of the exact and approximate solutions for t = 0.2, 0.4, 0.6, 0.8, 1 are given in Figure 1. The dot points are the approximate solutions, and the solid lines are the exact solutions.
We discuss another example with ν 2 ( t ) is a non-constant function.
Example 4.
Consider the following two-dimensional fractional IDE of the form
w x x ( x , t ) + D t 2 q w ( x , t ) D t q w ( x , t ) + w x ( x , t ) + ( t q + 1 ) w ( x , t ) I t q w ( x , t ) = g ( x , t )
with
w ( 0 , t ) = w ( 1 , t ) = 0 , t 0
and
w ( x , 0 ) = e 1 2 x sin ( 5 π x ) , D q w ( x , 0 ) = 0 , 0 x 1
where
g ( x , t ) = 4 e x 2 t q ( Γ ( 2 q + 1 ) sin ( 2 π x ) + Γ ( 3 q + 1 ) sin ( 3 π x ) sin ( 5 π x ) ) 4 Γ ( q + 1 )
+ e x 2 4 Γ ( 2 q + 1 ) 2 sin ( 2 π x ) t 3 q Γ ( 3 q + 1 ) 4 t 2 q Γ ( 3 q + 1 ) 2 sin ( 3 π x ) 4 Γ ( 2 q + 1 ) Γ ( 3 q + 1 )
+ e x 2 Γ ( 4 q + 1 ) t 2 q 4 t q 16 π 2 + 3 sin ( 2 π x ) 4 t 4 q Γ ( 3 q + 1 ) sin ( 3 π x ) 4 Γ ( 4 q + 1 )
+ e x 2 t 3 q 4 t q 36 π 2 + 3 sin ( 3 π x ) + 4 t q 100 π 2 + 3 sin ( 5 π x ) 4 Γ ( 4 q + 1 ) .
Then, the exact solution is
w ( x , t ) = e x 2 t 2 q sin ( 2 π x ) + e x 2 t 3 q sin ( 3 π x ) + e x 2 sin ( 5 π x ) .
Here, we use M 1 = M 2 = 40 . Then, the absolute error for the fractional derivative of order q is defined by
ϵ q ( x , t ) = w ( x , t ) w 40 , 40 ( x , t ) .
The absolute errors for different values of x and t are reported in Table 2.
We also compute the L 2 -error, which is given by
ϵ L 2 ( q ) = 0 1 0 1 w ( x , t ) w 40 , 40 ( x , t ) 2 d x d t .
The L 2 - errors for different values of q are reported in Table 3.
Additionally, the graphs of the exact and approximate solutions for t = 0.2 , 0.4 , 0.6 , 0.8 , 1 for q = 9 10 are given in Figure 2. The dot points are the approximate solutions, and the solid lines are the exact solutions.
We also show the graph of exact and approximate solutions at t = 0.5 for different values of q in Figure 3.

6. Conclusions

This article is divided into two main parts. In the first part, we use OMM to numerically explore the 2D fractional IDE. We derive a numerical approach and present a proof of these operational matrices. We test the numerical method on some examples to show its efficiency. Results are shown in Examples 1 and 2. We use various methods to measure the errors of the approximation, such as absolute and L 2 errors. In addition, we explore the effect of derivative order. Note the following when using these two examples:
1.
The absolute error for different values of x and t is on the order of 10 14 , which is very small, as seen in Table 1 and Table 2.
2.
The L 2 error is also on the order of 10 13 and is very small, as shown in Table 3.
3.
Approximate and exact solutions for various values of t agree, as shown in Figure 1 and Figure 2.
4.
As shown in Figure 3, the approximate solution for the derivative order q converges to the solution for q = 1 .
Part 2 examines the exact solutions for 1D and 2D fractional IDEs. We use the fractional Laplace transform to generate these exact solutions. All homogeneous cases are considered. We provide two examples for two purposes. First, we show how the results obtained are implemented, especially the exact solution for some 1D and 2D IDE classes. We then show that the exact fractional solution converges to the exact solution of the ordinary derivative as the order of the fractional derivative approaches one. These results are reported in Examples 3 and 4. Note that this work uses sequential Caputo derivative. It helps us fluently implement numerical and analytical approaches. In this paper, we show that our approach works efficiently for these types of problems. This can be generalized to other types of fractional IDEs and other scientific and engineering applications.

Author Contributions

Methodology, S.M.S.; Software, S.M.S.; Formal analysis, S.M.S.; Investigation, S.M.S.; Writing—original draft, R.M.K.; Writing—review editing, S.H.A.; Supervision, Z.S. All authors reviewed the results and approved the final version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The Authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work (grant code: 23UQU4310382DSR001). We would like to express our gratitude to the editor and the reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 1993; Volume 6. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1999; Volume 198. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  4. Zhang, S.; Wei, T. Sequential fractional derivative: Definition and its properties. J. Math. Anal. Appl. 2009, 353, 441–446. [Google Scholar]
  5. Atanackovic, T.M.; Pilipovic, S. On sequential fractional derivatives in the Caputo sense. Appl. Math. Comput. 2010, 216, 3452–3460. [Google Scholar]
  6. Li, X.; Liu, F.; Yang, X.J.; Baleanu, D. Sequential fractional derivative models and applications. Chaos Solitons Fractals 2021, 143, 110739. [Google Scholar]
  7. Zhong, W.P.; Chen, Y.Q.; Sun, X.J.; Wei, X.J. On the sequential fractional derivative and its application to fractional partial differential equations. Commun. Nonlinear Sci. Numer. Simul. 2019, 74, 257–268. [Google Scholar]
  8. Baleanu, D.; Gülsu, M.; Mohammadi, H. On the non-commutativity of sequential fractional derivatives. Appl. Math. Lett. 2013, 26, 393–397. [Google Scholar]
  9. Zhang, H.; Li, Y.; Shen, Y. An efficient and accurate numerical method for solving multi-term fractional differential equations based on sequential fractional derivatives. J. Comput. Phys. 2020, 409, 109346. [Google Scholar]
  10. Columbu, A.; Frassu, S.; Viglialoro, G. Refined criteria toward boundedness in an attraction-repulsion chemotaxis system with nonlinear productions. Appl. Anal. 2023, 1–17. [Google Scholar] [CrossRef]
  11. Li, T.; Frassu, S.; Viglialoro, G. Combining effects ensuring boundedness in an attraction–repulsion chemotaxis model with production and consumption. Z. Angew. Math. Phys. 2023, 74, 109. [Google Scholar] [CrossRef]
  12. Volterra, V. Theory of Functionals and of Integral and Integro-Differential Equations; Dover Publications: Mineola, NY, USA, 1913. [Google Scholar]
  13. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. (Eds.) Fractional Calculus Models and Numerical Methods; World Scientific: Singapore, 2012; Volume 3. [Google Scholar]
  14. Wang, Y.; Sun, H.; Li, C. Stability analysis of fractional differential equations with Caputo derivative. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1318–1326. [Google Scholar] [CrossRef]
  15. Wu, G.C.; Liao, S.J. An efficient approach to obtain higher-order approximations of fractional derivatives by the generalized Taylor matrix method. J. Vib. Control 2012, 18, 1013–1023. [Google Scholar]
  16. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods: Fundamentals in Single Domains; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  17. Abuasbeh, K.; Kanwal, A.; Shafqat, R.; Taufeeq, B.; Almulla, M.A.; Awadalla, M. A Method for Solving Time-Fractional Initial Boundary Value Problems of Variable Order. Symmetry 2023, 15, 519. [Google Scholar] [CrossRef]
  18. Abuasbeh, K.; Shafqat, R.; Alsinai, A.; Awadalla, M. Analysis of the Mathematical Modelling of COVID-19 by Using Mild Solution with Delay Caputo Operator. Symmetry 2023, 15, 286. [Google Scholar] [CrossRef]
  19. Dallashi, Q.; Syam, M.I. An Efficient Method for Solving Second-Order Fuzzy Order Fuzzy Initial Value Problems. Symmetry 2022, 14, 1218. [Google Scholar] [CrossRef]
  20. Tayeb, M.; Boulares, H.; Moumen, A.; Imsatfia, M. Processing Fractional Differential Equations Using ψ-Caputo Derivative. Symmetry 2023, 15, 955. [Google Scholar] [CrossRef]
  21. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  22. Baleanu, D.; Mustafa, O.G. On the solutions of sequential fractional differential equations. J. Comput. Nonlinear Dyn. 2013, 8, 031012. [Google Scholar]
  23. Vatsala, A.S.; Pageni, G.; Vijesh, V.A. Analysis of Sequential Caputo Fractional Differential Equations versus Non-Sequential Caputo Fractional Differential Equations with Applications. Foundations 2022, 2, 1129–1142. [Google Scholar] [CrossRef]
  24. Alomari, A.K.; Baleanu, D.; Al-Smadi, M.O. An efficient operational matrix method for solving multi-term time-fractional diffusion equations. Appl. Math. Comput. 2020, 367, 124820. [Google Scholar]
  25. Syam, M.I.; Sharadga, M.; Hashum, I. A numerical method for solving fractional delay differential equations based on the operational matrix. Chaos Solitons Fractals 2021, 147, 110977. [Google Scholar] [CrossRef]
  26. Kashkari, B.; Syam, M. Fractional-order Legendre operational matrix of fractional integration for solving the Riccati equation with fractional order. Appl. Math. Comput. 2016, 290, 281–291. [Google Scholar] [CrossRef]
Figure 1. The exact and approximate solutions for t = 0.2 , 0.4 , 0.6 , 8 , 1 .
Figure 1. The exact and approximate solutions for t = 0.2 , 0.4 , 0.6 , 8 , 1 .
Symmetry 15 01263 g001
Figure 2. The exact and approximate solutions for t = 0.2 , 0.4 , 0.6 , 8 , 1 .
Figure 2. The exact and approximate solutions for t = 0.2 , 0.4 , 0.6 , 8 , 1 .
Symmetry 15 01263 g002
Figure 3. The exact and approximate solutions at t = 0.5 for different values of q.
Figure 3. The exact and approximate solutions at t = 0.5 for different values of q.
Symmetry 15 01263 g003
Table 1. The absolute errors for different values of x and t for Example 3.
Table 1. The absolute errors for different values of x and t for Example 3.
x ϵ ( x , 0.25 ) ϵ ( x , 0.5 ) ϵ ( x , 0.75 )
0000
0.2 3.12 10 14 3.00 10 14 3.74 10 14
0.4 3.94 10 14 3.28 10 14 3.81 10 14
0.6 4.11 10 14 3.50 10 14 3.99 10 14
0.8 4.38 10 14 3.85 10 14 4.17 10 14
1 4.97 10 14 3.22 10 14 4.31 10 14
Table 2. The absolute errors for different values of x and t for Example 4.
Table 2. The absolute errors for different values of x and t for Example 4.
x ϵ 3 4 ( x , 0.25 ) ϵ 3 4 ( x , 0.5 ) ϵ 3 4 ( x , 0.75 )
0000
0.2 4.23 10 14 4.51 10 14 4.11 10 14
0.4 4.79 10 14 4.72 10 14 4.38 10 14
0.6 4.91 10 14 4.89 10 14 4.83 10 14
0.8 5.34 10 14 5.21 10 14 5.32 10 14
1 5.80 10 14 5.42 10 14 5.55 10 14
Table 3. The L 2 - errors for different values of q for Example 4.
Table 3. The L 2 - errors for different values of q for Example 4.
q ϵ L 2 ( q )
3 4 4.23 10 13
7 8 4.48 10 13
8 9 4.51 10 13
10 11 4.42 10 13
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Syam, S.M.; Siri, Z.; Altoum, S.H.; Md. Kasmani, R. Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations. Symmetry 2023, 15, 1263. https://doi.org/10.3390/sym15061263

AMA Style

Syam SM, Siri Z, Altoum SH, Md. Kasmani R. Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations. Symmetry. 2023; 15(6):1263. https://doi.org/10.3390/sym15061263

Chicago/Turabian Style

Syam, Sondos M., Z. Siri, Sami H. Altoum, and R. Md. Kasmani. 2023. "Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations" Symmetry 15, no. 6: 1263. https://doi.org/10.3390/sym15061263

APA Style

Syam, S. M., Siri, Z., Altoum, S. H., & Md. Kasmani, R. (2023). Analytical and Numerical Methods for Solving Second-Order Two-Dimensional Symmetric Sequential Fractional Integro-Differential Equations. Symmetry, 15(6), 1263. https://doi.org/10.3390/sym15061263

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