Next Article in Journal
Flavor Symmetry of Hydrogen Atoms Potentially Affecting the Proton Radius Deduced from the Electron-Hydrogen Scattering
Next Article in Special Issue
Editorial for the Special Issue of “Fractional Differential and Fractional Integro-Differential Equations: Qualitative Theory, Numerical Simulations, and Symmetry Analysis”
Previous Article in Journal
Analysis of Implicit Solutions for a Coupled System of Hybrid Fractional Order Differential Equations with Hybrid Integral Boundary Conditions in Banach Algebras
Previous Article in Special Issue
Systems of Sequential ψ1-Hilfer and ψ2-Caputo Fractional Differential Equations with Fractional Integro-Differential Nonlocal Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Semi-Discretization Method Based on Finite Difference and Differential Transform Methods to Solve the Time-Fractional Telegraph Equation

Department of Mathematics, University of Sistan and Baluchestan, Zahedan 98155-987, Iran
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(9), 1759; https://doi.org/10.3390/sym15091759
Submission received: 12 February 2023 / Revised: 15 April 2023 / Accepted: 12 September 2023 / Published: 13 September 2023

Abstract

:
The telegraph equation is a hyperbolic partial differential equation that has many applications in symmetric and asymmetric problems. In this paper, the solution of the time-fractional telegraph equation is obtained using a hybrid method. The numerical simulation is performed based on a combination of the finite difference and differential transform methods, such that at first, the equation is semi-discretized along the spatial ordinate, and then the resulting system of ordinary differential equations is solved using the fractional differential transform method. This hybrid technique is tested for some prominent linear and nonlinear examples. It is very simple and has a very small computation time; also, the obtained results demonstrate that the exact solutions are exactly symmetric with approximate solutions. The results of our scheme are compared with the two-dimensional differential transform method. The numerical results show that the proposed method is more accurate and effective than the two-dimensional fractional differential transform technique. Also, the implementation process of this method is very simple, so its computer programming is very fast.

1. Introduction

The differential transform method (DTM) is an iterative method based on Taylor’s series. DTM has been used to solve various differential equations. It was first applied to solve electrical circuit problems. After that, it was used to solve ordinary differential equations (ODEs), partial differential equations (PDEs), fuzzy PDEs, fractional-order ODEs and PDEs, systems of ODEs, systems of PDEs, differential-algebraic equations, and eigenvalue problems [1,2,3,4,5,6,7]. Additionally, fractional DTM, which is based on a generalized Taylor’s series, has been applied to solve various differential, differential-algebraic, and integral equations of fractional order [8,9,10,11,12,13,14]. In this paper, we intend to apply a combination of the finite difference (FD) and fractional differential transform (FDT) methods (FD-FDTM) to solve the one-dimensional time-fractional telegraph equation (FTE).
We consider the FTE in the following form [15,16,17]:
2 γ v ( x , t ) t 2 γ + 2 λ γ v ( x , t ) t γ + μ v ( x , t ) = ν 2 v ( x , t ) x 2 + q ( x , t ) , a < x < b , t 0 ,
with initial conditions
v ( x , 0 ) = f 1 ( x ) , v t ( x , 0 ) = f 2 ( x ) , a < x < b ,
and boundary conditions
v ( a , t ) = g 1 ( t ) , v ( b , t ) = g 2 ( t ) , t 0 .
where 0 < γ < 1 , λ , μ , and ν are arbitrary positive constants and f 1 , f 2 , g 1 , and g 2 are known functions.
Also, 2 γ v ( x , t ) t 2 γ and γ v ( x , t ) t γ denote the Caputo fractional derivative of order 2 γ and γ , respectively. The γ -order Caputo fractional derivative of the function f for γ > 0 , n 1 < γ < n , is defined as follows:
D 0 γ f ( x ) = 1 Γ ( n γ ) 0 x ( x t ) n γ 1 f ( n ) ( t ) d t .
The telegraph equation is a hyperbolic PDE that has many applications in physics and engineering, for example, in signal analysis, random walk theory, anomalous diffusion processes, wave phenomena, and wave propagation of the electrical signal in the cable of a transmission line. Different numerical and analytical techniques have been used to solve fractional-order telegraph equations [15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33].
This work aims to obtain an approximate solution to the FTE (1) using a hybrid method. In 2008, a hybrid method based on the combination of DTM and FDM was presented to solve a nonlinear heat conduction differential equation [34]. Also, in 2012, some nonlinear PDEs were solved with the hybrid method [35]. Arsalan (2020) applied a hybrid scheme to solve the one-dimensional integer-order telegraph equations [36,37].
We organize the rest of the paper as follows. In Section 2, we present the FDTM and related theorems. We propose our hybrid method to solve FTE in Section 3 and prove its convergence in Section 4. Also, in Section 5, we give some examples and solve them with the proposed method, we draw a conclusion in Section 6.

2. Fractional Differential Transform Method

The fractional differential transform method (or generalized differential transform method) is based on the fractional Taylor’s formula. The α -order fractional Taylor expansion of function u ( t ) about point t = t 0 is defined as [38]
u ( t ) = k = 0 ( t t 0 ) k α Γ ( k α + 1 ) ( d α d t α ) k u ( t ) t = t 0 ,
where d α d t α is the α -order Caputo fractional derivative and ( d α d t α ) k = d α d t α d α d t α k t i m e s .
The α -order FDT of the function u ( t ) about t = t 0 is denoted by U α ( k ) and defined as U α ( k ) = 1 Γ ( k α + 1 ) ( d α d t α ) k u ( t ) t = t 0 , and the inverse transform is u ( t ) = k = 0 U α ( k ) ( t t 0 ) k α [9]. Therefore, at t = 0 , we have
u ( t ) = k = 0 U α ( k ) t k α = k = 0 φ k ( t ) ,
where φ k ( t ) = U α ( k ) t k α .
Also, the m-approximation fractional differential transform of u ( t ) is defined as
U m ( t ) = k = 0 m U α ( k ) t k α = k = 0 m φ k ( t ) .
Theorem 1 
([39]). Suppose that F α ( k ) , G α ( k ) , and H α ( k ) are the differential transformations of the functions f ( t ) , g ( t ) , and h ( t ) , respectively. Then we have
(a) 
if f ( t ) = g ( t ) ± h ( t ) , then F α ( k ) = G α ( k ) ± H α ( k ) ,
(b) 
if f ( t ) = ( t t 0 ) q , then F α ( k ) = δ ( k q α ) , where δ ( k ) = { 1 , i f k = 0 0 . i f k 0
(c) 
if f ( t ) = g ( t ) h ( t ) , then F α ( k ) = l = 0 k G α ( l ) H α ( k l ) ,
Theorem 2 
([39]). Suppose that f ( t ) = t λ g ( t ) , where λ > 1 and g ( t ) has the generalized power series expansion g ( t ) = n = 0 a n ( t t 0 ) n α with radius of convergence R > 0 , 0 < α 1 . Then
D a γ D a β f ( t ) = D a γ + β f ( t ) ,
for all t ( 0 , R ) if:
(a) 
β < λ + 1 and α arbitrary or
(b) 
β λ + 1 , γ arbitrary, and a k = 0 for k = 0 , 1 , , m 1 , where m 1 < β m .
Theorem 3 
([39]). If f ( t ) = D t 0 γ g ( t ) , m 1 < γ m , and the function g ( t ) satisfies the conditions in theorem (2), then
F α ( k ) = Γ ( k α + γ + 1 ) Γ ( k α + 1 ) G α ( k + γ α ) .

3. FD-FDTM for Solving the FTE

Consider the FTE (1). If the x-derivative at ( x , t ) is replaced by 1 h 2 v ( x h , t ) 2 v ( x , t ) + v ( x + h , t ) + O ( h 2 ) and x is considered as a constant, Equation (1) can be written as the following ordinary differential equation
d 2 γ v ( t ) d t 2 γ + 2 λ d γ v ( t ) d t γ + μ v ( t ) = ν 1 h 2 v ( x h , t ) 2 v ( x , t ) + v ( x + h , t ) + O ( h 2 ) + q ( x , t ) .
We subdivide the interval [ a , b ] into N equal subintervals of step-length h = b a N . Thus, the mesh points x i = a + i h , i = 0 , 1 , , N are obtained. Now, we write Equation (6) at the mesh point x i , i = 1 , , N 1 , along with time level t. If we discard the local truncation error O ( h 2 ) and denote u i ( t ) as the approximate solution of v i ( t ) = v ( x i , t ) , we have the following system of ODEs:
d 2 γ u i ( t ) d t 2 γ + 2 λ d γ u i ( t ) d t γ + μ u i ( t ) = ν h 2 u i 1 ( t ) 2 u i ( t ) + u i + 1 ( t ) + q i ( t ) , i = 1 , , N 1 .
We solve the system (7) using FDTM. For this purpose, we consider the solution of equation, u i ( t ) , as follows:
u i ( t ) = k = 0 U i α ( k ) t k α ,
where the unknown coefficients U i α ( k ) are the FDT of u i ( t ) and should be obtained.
By choosing a suitable value for α , assuming q i α ( k ) as the α -fractional differential transform of q i ( t ) , and by using Theorems 2 and 3, the fractional differential transform of Equation (7) leads to the following relation:
Γ ( k α + 2 γ + 1 ) Γ ( k α + 1 ) U i α ( k + 2 γ α ) + 2 λ Γ ( k α + γ + 1 ) Γ ( k α + 1 ) U i α ( k + γ α ) + μ U i α ( k ) = ν h 2 { U i 1 α ( k ) 2 U i α ( k ) + U i + 1 α ( k ) } + q i α ( k ) .
Now, suppose that G 1 α ( k ) and G 2 α ( k ) are the fractional differential transforms of the functions g 1 ( t ) and g 2 ( t ) , respectively. Therefore, by applying the FDTM to conditions (2) and (3), we have the initial conditions
U i α ( 0 ) = f 1 ( x i ) , U i α ( 1 α ) = f 2 ( x i ) ,
and the boundary conditions
U 0 α ( k ) = G 1 α ( k ) , U N α ( k ) = G 2 α ( k ) ,
We rewrite relations (9)–(11) as follows:
U i α ( k + 2 γ α ) = Γ ( k α + 1 ) Γ ( k α + 2 γ + 1 ) ν U i 1 α ( k ) 2 U i α ( k ) + U i + 1 α ( k ) h 2 μ U i α ( k ) + q i α ( k )
2 λ Γ ( k α + γ + 1 ) Γ ( k α + 2 γ + 1 ) U i α ( k + γ α ) , k 0
U i α ( 0 ) = f 1 ( i h ) , U i α ( 1 α ) = f 2 ( i h ) ,
U 0 α ( k ) = G 1 α ( k ) , U N α ( k ) = G 2 α ( k ) .
Also, according to [9], the unknown coefficients U i α ( 1 ) , U i α ( 2 ) , , U i α ( 2 γ α 1 ) , will be available as follows:
U i α ( k ) = 1 Γ ( k α + 1 ) [ d k α d t k α u i ( t ) ] t = 0 , if k α Z + k = 0 , 1 , ( 2 γ α 1 ) . 0 , if k α Z +
Therefore, all the unknown coefficients U i α ( k ) , i = 0 , 1 , 2 , , N , k 0 are calculated according to the recursive formula (12) and relations (13)–(15).

4. Convergence of FD-FDTM for FTE

Here, we discuss the convergence of FD-FDTM for solving the FTE (1). First, we present the following lemma [40].
Lemma 1. 
Suppose that for some k 0 N 0 and for every j k 0 , there exist 0 < δ j < 1 such that φ j + 1   δ j + 1 φ j , ( φ j   = m a x t | φ j ( t ) | ) . Then the series k = 0 φ k ( t ) converges to u ( t ) .
Proof. 
Consider the sequence s 0 , s 1 , s 2 , , where s n = k = 0 n φ k ( t ) . To prove the lemma, we show the sequence s n is a Cauchy sequence in ( C [ 0 , 1 ] , . ) .
For 0 < δ j < 1 , we can write
s j s j 1   =   φ j   δ j φ j 1   δ j δ j 1 φ j 2     δ j δ j 1 δ k 0 φ k 0 ,
Thus, for n m k 0 , we have
s n s m   =   j = m + 1 n ( s j s j 1 )   j = m + 1 n s j s j 1   j = m + 1 n δ j δ j 1 δ k 0 φ k 0 ,
If we let δ = m a x { δ k 0 , δ k 0 + 1 , , δ m , δ m + 1 , , δ n } , the following relation is obtained:
s n s m   j   =   m + 1 n δ j k 0 + 1 φ k 0     φ k 0 δ m k 0 + 2 1 + δ + δ 2 + + δ n m 1 = 1 δ n m 1 δ δ m k 0 + 2 φ k 0 .
Since 0 δ < 1 , we can derive lim n , m s n s m = 0 , which means { s n } n = 0 is a Cauchy sequence in ( C [ I ] , . ) . Since the space C [ I ] with . is a Banach space, we can derive that the series k = 0 φ k ( t ) is convergent to u ( t ) . □
Theorem 4. 
Suppose that v ( x i , t ) is the exact solution of the FTE at point ( x i , t ) , u i ( t ) is the exact solution of Equation (7), and U i m ( t ) = k = 0 m φ k ( t ) , is the m-approximation of u i ( t ) as the approximate solution of FTE at point ( x i , t ) . Also, suppose for some k 0 N 0 and for every n m k 0 , 0 < δ i < 1 , such that φ i + 1   δ i + 1 φ i , where φ i   = m a x t | φ i ( t ) | . Then the solution U i m ( t ) converges to the exact solution, v i ( t ) , as m . Furthermore, for some a < ξ < b the maximum absolute error of the m-series, U i m ( t ) , as an approximation of the FTE’s exact solution satisfies the following relation:
v ( x i , t ) U i m ( t )   h 2 12 4 u ( ξ , t ) x 4 + 1 1 δ δ m k 0 + 2 φ k 0 ,
where δ = m a x { δ k 0 , δ k 0 + 1 , , δ n } .
Proof. 
We can write
v ( x i , t ) U i m ( t )   =   v ( x i , t ) u i ( t ) + u i ( t ) U i m ( t )     v ( x i , t ) u i ( t )   +   u i ( t ) U i m ( t ) ,
where v ( x i , t ) and u i ( t ) are the solutions of Equations (6) and (7), respectively.
Also, Equation (6) has been obtained by replacing the second x-derivative of v ( x , t ) with the central FD formula in Equation (1). Therefore, for some a < ξ < b we can write
v ( x i , t ) u i ( t )   h 2 12 4 u ( ξ , t ) x 4
From relation (16), for n m k 0 , we have
s n s m   1 δ n m 1 δ δ m k 0 + 2 φ k 0 ,
and since 0 δ < 1 , then 1 δ n m < 1 , so we have
s n s m   1 1 δ δ m k 0 + 2 φ k 0 .
If n approaches , then s n u i ( t ) and we have
u i ( t ) s m   1 1 δ δ m k 0 + 2 φ k 0 ,
in the other words,
u i ( t ) U i m ( t )   1 1 δ δ m k 0 + 2 φ k 0 .
By replacing relations (18) and (19) in relation (17), the theorem is proved. □

5. Numerical Examples

In this section, we give some examples to show the efficiency and convenience of the mentioned method. The examples include linear and non-linear FTEs. We present the results of FD-FDTM for solving the examples and calculate the maximum absolute error (MAE) for different values of N using the following formula:
E ( N ) = max 1 i N | v ( x i , t ) U i m ( t ) | .
Also, we compare the results of FD-FDTM with two-dimensional FDTM (2D-FDTM). Moreover, we obtain the rate of convergence (ROC) of FD-FDTM with the following formula:
R O C = log 2 ( E ( N ) E ( 2 N ) ) .
Example 1. 
Consider the following linear FTE [30]:
2 γ v ( x , t ) t 2 γ + 2 γ v ( x , t ) t γ + v ( x , t ) = 2 v ( x , t ) x 2 , 0 < γ < 1 , 0 < x < 1 , t > 0 ,
with the following conditions
v ( x , 0 ) = e x , v t ( x , 0 ) = 2 e x ,
v ( 0 , t ) = e 2 t , v ( 1 , t ) = e 1 2 t ,
which has the exact solution v ( x , t ) = e x 2 t for γ = 1 .
We describe the mentioned method to solve Equations (20)–(22) for γ = 0.75 , α = 0.25 , and h = 0.1 ( N = 10 ). Using Theorem 3, we can obtain the differential transform of the derivatives in Equation (20) as follows:
1.5 v ( x , t ) t 1.5 Γ ( 0.25 k + 2.5 ) Γ ( 0.25 k + 1 ) U i 0.25 ( k + 6 ) , 0.75 v t 0.75 Γ ( 0.25 k + 1.75 ) Γ ( 0.25 k + 1 ) U i 0.25 ( k + 3 ) , v ( x , t ) U i 0.25 ( k ) , 2 v x 2 U i 1 0.25 ( k ) 2 U i 0.25 ( k ) + U i + 1 0.25 ( k ) h 2 .
According to relation (13), for the initial conditions we have
v ( x , 0 ) = e x U i 0.25 ( 0 ) = e x i , i = 0 , 1 , , 10 , v t ( x , 0 ) = 2 e x U i 0.25 ( 4 ) = 2 e x i , i = 0 , 1 , , 10 ,
and according to relation (15), we have
U i 0.25 ( 1 ) = U i 0.25 ( 2 ) = U i 0.25 ( 3 ) = U i 0.25 ( 5 ) = 0 , i = 0 , 1 , 2 , , 10 .
Also, we use relation (14) for the boundary conditions, and for the right side of these relations, we use the FDT of the e 2 t , e 1 2 t functions, so we obtain
v ( 0 , t ) = e 2 t U 0 0.25 ( k ) = ( 2 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 . i f k 4 Z +
v ( 1 , t ) = e 1 2 t U 10 0.25 ( k ) = e × ( 2 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 . i f k 4 Z +
By replacing the above relations in Equation (20), we obtain the following recursive relationship:
U i 0.25 ( k + 6 ) = Γ ( 0.25 k + 1 ) Γ ( 0.25 k + 2.5 ) U i 1 0.25 ( k ) 2 U i 0.25 ( k ) + U i + 1 0.25 ( k ) h 2 U i 0.25 ( k ) 2 Γ ( 0.25 k + 1.75 ) Γ ( 0.25 k + 2.5 ) U i 0.25 ( k + 3 ) . i 1
Thus, the solution in ( x i , t ) , i = 0 , 1 , , 10 is obtained as follows:
x 0 = 0 , u 0 ( t ) = k = 0 U 0 0.25 ( k ) t 0.25 k = U 0 0.25 ( 0 ) + U 0 0.25 ( 1 ) t 0.25 + U 0 0.25 ( 2 ) t 0.5 + U 0 0.25 ( 3 ) t 0.75 + U 0 0.25 ( 4 ) t + , = 1 + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 t + , x 1 = 0.1 , u 1 ( t ) = k = 0 U 1 0.25 ( k ) t 0.25 k = U 1 0.25 ( 0 ) + U 1 0.25 ( 1 ) t 0.25 + U 1 0.25 ( 2 ) t 0.5 + U 1 0.25 ( 3 ) t 0.75 + U 1 0.25 ( 4 ) t + = e 0.1 + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 e 0.1 t + , x 10 = 1 , u 10 ( t ) = k = 0 U 10 0.25 ( k ) t 0.25 k = U 10 0.25 ( 0 ) + U 10 0.25 ( 1 ) t 0.25 + U 10 0.25 ( 2 ) t 0.5 + U 10 0.25 ( 3 ) t 0.75 + U 10 0.25 ( 4 ) t + , = e + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 e t + .
We show the results of our method for solving Example 1 in Table 1 and Table 2. Table 1 contains the maximum absolute error of the obtained solution using the FD-FDTM for γ = 1 , m = 3 , and different values of N at t = 0.001 . Also, we compared the results of FD-FDTM with two-dimensional FDTM. Table 1 shows the our method is more accurate than the two-dimensional DTM. Also, we can see that as N increases, the error decreases and the numerical ROC confirms the theoretical ROC. Table 2 compares the approximate solution of FD-FDTM and 2D-FDTM at t = 0.01 , for γ = 0.75 , m = 10 , and N = 10 . Figure 1 shows the comparison between the exact solution for γ = 1 and the results of FD-FDTM for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Example 2. 
In this example, we consider a non-homogeneous FTE [17]
2 γ v ( x , t ) t 2 γ + 40 γ v ( x , t ) t γ + 100 v ( x , t ) = 2 v ( x , t ) x 2 + 23 e 2 t s i n h ( x ) , 0 < γ < 1 , 0 < x < 1 , t > 0 ,
with the following conditions:
v ( x , 0 ) = s i n h ( x ) , v t ( x , 0 ) = 2 s i n h ( x ) ,
v ( 0 , t ) = 0 , v ( 1 , t ) = e 2 t s i n h ( 1 ) .
The exact solution of (23) for γ = 1 with conditions (24) and (25) is v ( x , t ) = e 2 t s i n h ( x ) .
For example, we put γ = 0.6 in (23) and consider h = 0.1 and α = 0.1 . For these values of γ , h , and α, we can write
1.2 v ( x , t ) t 1.2 Γ ( 0.1 k + 2.2 ) Γ ( 0.1 k + 1 ) U i 0.1 ( k + 12 ) , 0.6 v t 0.6 Γ ( 0.1 k + 1.6 ) Γ ( 0.1 k + 1 ) U i 0.1 ( k + 6 ) , v ( x , t ) U i 0.1 ( k ) , 2 v x 2 U i 1 0.1 ( k ) 2 U i 0.1 ( k ) + U i + 1 0.1 ( k ) h 2 .
According to relation (13), for the initial conditions we have:
v ( x , 0 ) = s i n h ( x ) U i 0.1 ( 0 ) = s i n h ( x i ) , i = 0 , 1 , 2 , , 10 , v t ( x , 0 ) = 2 s i n h ( x ) U i 0.1 ( 10 ) = 2 s i n h ( x i ) , i = 0 , 1 , 2 , , 10 ,
and according to (15), we have
U i 0.1 ( 1 ) = U i 0.1 ( 2 ) = = U i 0.1 ( 9 ) = U i 0.1 ( 11 ) = 0 , i = 0 , 1 , , 10 .
We use relation (14) for the boundary conditions, and for the right side of the boundary conditions, we use the FDT of the e 2 t s i n h ( 1 ) function. Therefore, we have
v ( 0 , t ) = 0 U 0 0.1 ( k ) = 0 , k 0
U 10 0.1 ( k ) = s i n h ( 1 ) × ( 2 ) k 10 Γ ( k 10 + 1 ) , i f k 10 Z + k 0 0 . i f k 10 Z +
By putting the above relations in Equation (23), we conclude the following recursive relationship:
U i 0.1 ( k + 12 ) = Γ ( 0.1 k + 1 ) Γ ( 0.1 k + 2.2 ) 23 E 0.1 ( k ) s i n h ( x i ) 100 U i 0.1 ( k ) + U i 1 0.1 ( k ) 2 U i 0.1 ( k ) + U i + 1 0.1 ( k ) h 2 40 Γ ( 0.1 k + 1.6 ) Γ ( 0.1 k + 2.2 ) U i 0.1 ( k + 6 ) ,
where E 0.1 ( k ) is the fractional differential transform of e 2 t and can be obtained as follows:
E 0.1 ( k ) = ( 2 ) k 10 Γ ( k 10 + 1 ) , i f k 10 Z + k 0 0 . i f k 10 Z +
Thus, in the points ( x i , t ) , i = 0 , 1 , , N , we can write:
x 0 = 0 , u 0 ( t ) = k = 0 U 0 0.1 ( k ) t 0.1 k = U 0 0.1 ( 0 ) + U 0 0.1 ( 1 ) t 0.1 + U 0 0.1 ( 2 ) t 0.2 + U 0 0.1 ( 3 ) t 0.3 + U 0 0.1 ( 4 ) t 0.4 + , = 0 + 0 t 0.1 + 0 t 0.2 + 0 t 0.3 + 0 t 0.4 + , x 1 = 0.1 , u 1 ( t ) = k = 0 U 1 0.1 ( k ) t 0.1 k = U 1 0.1 ( 0 ) + U 1 0.1 ( 1 ) t 0.1 + U 1 0.1 ( 2 ) t 0.2 + + U 1 0.1 ( 9 ) t 0.9 + U 1 0.1 ( 10 ) t + , = s i n h ( 0.1 ) + 0 t 0.1 + 0 t 0.2 + + 0 t 0.9 2 s i n h ( 0.1 ) t + , x 10 = 1 , u 10 ( t ) = k = 0 U 10 0.1 ( k ) t 0.1 k = U 10 0.1 ( 0 ) + U 10 0.1 ( 1 ) t 0.1 + U 10 0.1 ( 2 ) t 0.2 + + U 10 0.1 ( 9 ) t 0.9 + U 10 0.1 ( 10 ) t + , = s i n h ( 1 ) + 0 t 0.1 + 0 t 0.2 + + 0 t 0.9 2 s i n h ( 1 ) t + .
To compare the exact and numerical solution of Example 2 obtained using FD-FDTM and 2D-FDTM, we compute the MAE of the obtained solution at t = 0.001 , and present them in Table 3. Also, we show the results of the FD-FDTM and 2D-FDTM at t = 0.001 , for γ = 0.6 , m = 15 , and N = 10 in Table 4. Figure 2 shows the comparison between the exact solution for γ = 1 and the numerical solution for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Example 3. 
In this example, we solve the following nonlinear FTE [23]:
2 γ v ( x , t ) t 2 γ + 2 γ v ( x , t ) t γ + u 2 ( x , t ) = 2 v ( x , t ) x 2 + e 2 x 4 t e x 2 t , 0 < γ < 1 , 0 < x < 1 , t > 0 ,
with the following conditions
v ( x , 0 ) = e x , v t ( x , 0 ) = 2 e x ,
v ( 0 , t ) = e 2 t , v ( 1 , t ) = e 1 2 t .
The exact solution of Equation (26) for γ = 1 with conditions (27) and (28) is v ( x , t ) = e x 2 t .
To explain the FD-FDTM, we put γ = 0.75 in (26) and consider h = 0.1 and α = 0.25 . For these values of γ , h , and α, we can write
1.5 v ( x , t ) t 1.5 Γ ( 0.25 k + 2.5 ) Γ ( 0.25 k + 1 ) U i 0.25 ( k + 6 ) , 0.75 v t 0.75 Γ ( 0.25 k + 1.75 ) Γ ( 0.25 k + 1 ) U i 0.25 ( k + 3 ) , v ( x , t ) U i 0.25 ( i , k ) , 2 v x 2 U i 1 0.25 ( k ) 2 U i 0.25 ( k ) + U i + 1 0.25 ( k ) h 2 .
According to relation (13), for the initial conditions we have
v ( x , 0 ) = e x U i 0.25 ( 0 ) = e x i , i = 0 , 1 , , 10 , v t ( x , 0 ) = 2 e x U i 0.25 ( 4 ) = 2 e x i , i = 0 , 1 , , 10 .
From relationship (15), we can write
U i 0.25 ( 1 ) = U i 0.25 ( 2 ) = U i 0.25 ( 3 ) = U i 0.25 ( 5 ) = 0 , i = 0 , 1 , 2 , , 10 .
We use relation (14) for the boundary conditions and an FDT of e 2 t , e 1 2 t for the right side of the boundary conditions. Therefore, we have
v ( 0 , t ) = e 2 t U 0 0.25 ( k ) = ( 2 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 . i f k 4 Z +
v ( 1 , t ) = e 1 2 t U 10 0.25 ( k ) = e × ( 2 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 . i f k 4 Z +
In Example 3, we have q ( x , t ) = e 2 x 4 t e x 2 t , so by using FD-FDTM, we have:
q i 0.25 ( k ) = e 2 x i × E 1 0.25 ( k ) e x i × E 2 0.25 ( k ) , i = 0 , 1 , , 10 ,
where E 1 0.25 ( k ) and E 2 0.25 ( k ) are the FDT of the e 4 t and e 2 t functions, respectively, and are obtained as follows:
E 1 0.25 ( k ) = ( 4 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 , i f k 4 Z +
and
E 2 0.25 ( k ) = ( 2 ) k 4 Γ ( k 4 + 1 ) , i f k 4 Z + k 0 0 . i f k 4 Z +
By substituting the above relations in Equation (26), we obtain the following recursive formula
U i 0.25 ( k + 6 ) = Γ ( 0.25 k + 1 ) Γ ( 0.25 k + 2.5 ) q i 0.25 ( k ) + U i 1 0.25 ( k ) 2 U i 0.25 ( k ) + U i + 1 0.25 ( k ) h 2 j = 0 k U i 0.25 ( j ) U i 0.25 ( k j ) 2 Γ ( 0.25 k + 1.75 ) Γ ( 0.25 k + 2.5 ) U i 0.25 ( k + 3 ) . i 1
In the points ( x i , t ) , i = 0 , 1 , , N , we have:
x 0 = 0 , u 0 ( t ) = k = 0 U 0 0.25 ( k ) t 0.25 k = U 0 0.25 ( 0 ) + U 0 0.25 ( 1 ) t 0.25 + U 0 0.25 ( 2 ) t 0.5 + U 0 0.25 ( 3 ) t 0.75 + U 0 0.25 ( 4 ) t + = 1 + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 t + , x 1 = 0.1 , u 1 ( t ) = k = 0 U 1 0.25 ( k ) t 0.25 k = U 1 0.25 ( 0 ) + U 1 0.25 ( 1 ) t 0.25 + U 1 0.25 ( 2 ) t 0.5 + U 1 0.25 ( 3 ) t 0.75 + U 1 0.25 ( 4 ) t + = e 0.1 + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 e 0.1 t + , x 10 = 1 , u 10 ( t ) = k = 0 U 10 0.25 ( k ) t 0.25 k = U 10 0.25 ( 0 ) + U 10 0.25 ( 1 ) t 0.25 + U 10 0.25 ( 2 ) t 0.5 + U 10 0.25 ( 3 ) t 0.75 + U 10 0.25 ( 4 ) t + = e + 0 t 0.25 + 0 t 0.5 + 0 t 0.75 2 e t + .
We show the results of our method for solving Example 3 in Table 5. Table 5 contains the MAE of the solution obtained using FD-FDTM, for γ = 1 , m = 3 , and different values of N at t = 0.01 . Also, we compared the results of FD-FDTM with two-dimensional FDTM. Table 5 shows that the FD-FDTM is more accurate than the two-dimensional FDTM. Figure 3 compares the exact solution for γ = 1 with the numerical solution for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Example 4. 
Consider the following time-fractional multi-term wave equation [32]:
1.7 v ( x , t ) t 1.7 1 2 v ( x , t ) t = 2 v ( x , t ) x 2 , 0 < x < 4 , t > 0 ,
with initial conditions
v ( x , 0 ) = s i n ( π x ) , v t ( x , 0 ) = 0 ,
and boundary conditions
v ( 0 , t ) = 0 , v ( 4 , t ) = 0 ,
We describe the mentioned method to solve Equations (29)–(31) for α = 0.1 By using theorem (3), we can obtain the differential transform of the derivatives in Equation (29) as follows:
1.7 v ( x , t ) t 1.7 Γ ( 0.1 k + 2.7 ) Γ ( 0.1 k + 1 ) U i 0.1 ( k + 17 ) , v t Γ ( 0.1 k + 2 ) Γ ( 0.1 k + 1 ) U i 0.1 ( k + 10 ) , 2 v x 2 U i 1 0.1 ( k ) 2 U i 0.1 ( k ) + U i + 1 0.1 ( k ) h 2 .
According to relation (13), for the initial conditions we have
v ( x , 0 ) = s i n ( π x ) U i 0.1 ( 0 ) = s i n ( π x i ) , i = 0 , 1 , , N , v t ( x , 0 ) = 0 U i 0.1 ( 10 ) = 0 , i = 0 , 1 , , N ,
and according to relation (15) we have
U i 0.1 ( 1 ) = U i 0.1 ( 2 ) = = U i 0.1 ( 15 ) = U i 0.1 ( 16 ) = 0 , i = 0 , 1 , , N ,
and
v ( 0 , t ) = 0 U 0 0.1 ( k ) = 0 , v ( 4 , t ) = 0 U N 0.1 ( k ) = 0 .
By replacing the above relations in Equation (29), we obtain the following recursive relationship:
U i 0.1 ( k + 17 ) = Γ ( 0.1 k + 1 ) Γ ( 0.1 k + 2.7 ) U i 1 0.1 ( k ) 2 U i 0.1 ( k ) + U i + 1 0.1 ( k ) h 2 + 1 2 Γ ( 0.1 k + 2 ) Γ ( 0.1 k + 2.5 ) U i 0.1 ( k + 10 ) . i 1
We show the results of our method for solving Example 4 in Figure 4 and Figure 5 for t = 1 .

6. Conclusions

In this work, a hybrid method has been used to solve the linear and non-linear time-fractional telegraph equation approximately. The central difference method has been applied to discretize the spatial derivative and the fractional differential transform method has been used to solve the obtained system of fractional ordinary equations. A convergence analysis of the mentioned method has been conducted. The numerical results show that the proposed hybrid method is more accurate and effective than two-dimensional FDTM. Also, the implementation process of this method is very simple, so its computer programming is very fast.

Author Contributions

Conceptualization, M.A. and Z.S.; methodology, M.A. and Z.S.; software, Z.S.; validation, M.A. and Z.S.; formal analysis, M.A. and Z.S.; investigation, M.A. and Z.S.; resources, Z.S.; writing—original draft preparation, M.A. and Z.S.; writing—review and editing, M.A.; supervision, M.A.; project administration, M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, C.K.; Ho, S.H. Application of differential transformation to eigenvalue problems. Appl. Math. Comput. 1996, 79, 173–188. [Google Scholar] [CrossRef]
  2. Zhou, J.K. Differential Transformation and Its Applications for Electrical Circuits; Huazhong University Press: Wuhan, China, 1986. [Google Scholar]
  3. Jang, M.J.; Chen, C.L.; Liy, Y.C. On solving the initial-value problems using the differential transformation method. Appl. Math. Comput. 2000, 115, 145–160. [Google Scholar] [CrossRef]
  4. Ayaz, F. Applications of differential transform method to differential-algebraic equations. Appl. Math. Comput. 2004, 152, 649–657. [Google Scholar] [CrossRef]
  5. Hassan, I.A.H. On solving some eigenvalue problems by using a differential transformation. Appl. Math. Comput. 2002, 127, 1–22. [Google Scholar]
  6. Alquran, M.T. Applying differential transform method to nonlinear partial differential equations: A modified approach. Appl. Appl. Math. Int. J. (AAM) 2012, 7, 10. [Google Scholar]
  7. Mirzaee, F.; Yari, M.K. A novel computing three-dimensional differential transform method for solving fuzzy partial differential equations. Ain Shams Eng. J. 2016, 7, 695–708. [Google Scholar] [CrossRef]
  8. 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]
  9. Odibat, Z.; Momani, S. A generalized differential transform method for linear partial differential equations of fractional order. Appl. Math. Lett. 2008, 21, 194–199. [Google Scholar] [CrossRef]
  10. Ibis, B.; Bayram, M.; Agargun, A.G. Applications of fractional differential transform method to fractional differential-algebraic equations. Eur. J. Pure Appl. Math. 2011, 4, 129–141. [Google Scholar]
  11. Ertürk, V.S. Computing eigenelements of Sturm–Liouville problems of fractional order via fractional differential transform method. Math. Comput. Appl. 2011, 16, 712–720. [Google Scholar] [CrossRef]
  12. Secer, A.; Akinlar, M.A.; Cevikel, A. Efficient solutions of systems of fractional PDEs by the differential transform method. Adv. Differ. Equ. 2012, 2012, 1–7. [Google Scholar] [CrossRef]
  13. Zain, F.A.S. Comparison study between differential transform method and Adomian decomposition method for some delay differential equations. Int. J. Phys. Sci. 2013, 8, 744–749. [Google Scholar]
  14. Rahimi, E.; Taghvafard, H.; Erjaee, G.H. Fractional differential transform method for solving a class of weakly singular Volterra integral equations. Iran. J. Sci. Technol. 2014, 38, 69. [Google Scholar]
  15. Hosseini, V.R.; Chen, W.; Avazzadeh, Z. Numerical solution of fractional telegraph equation by using radial basis functions. Eng. Anal. Bound. Elem. 2014, 38, 31–39. [Google Scholar] [CrossRef]
  16. Asghari, M.; Ezzati, R.; Allahviranloo, T. Numerical solutions of time-fractional order telegraph equation by Bernstein polynomials operational matrices. Math. Probl. Eng. 2016, 2016, 1683849. [Google Scholar] [CrossRef]
  17. Akram, T.; Abbas, M.; Ismail, A.I.; Ali, N.H.M.; Baleanu, D. Extended cubic B-splines in the numerical solution of time fractional telegraph equation. Adv. Differ. Equ. 2019, 2019, 365. [Google Scholar] [CrossRef]
  18. Orsingher, E.; Beghin, L. Time-fractional telegraph equations and telegraph processes with Brownian time. Probab. Theory Relat. Fields 2004, 128, 141–160. [Google Scholar] [CrossRef]
  19. Chen, J.; Liu, F.; Anh, V. Analytical solution for the time-fractional telegraph equation by the method of separating variables. J. Math. Anal. Appl. 2008, 338, 1364–1377. [Google Scholar] [CrossRef]
  20. Biazar, J.; Eslami, M. Analytic solution for telegraph equation by differential transform method. Phys. Lett. A 2010, 374, 2904–2906. [Google Scholar] [CrossRef]
  21. Garg, M.; Manohar, P.; Kalla, S.L. Generalized differential transform method to space-time fractional telegraph equation. Int. J. Differ. Equ. 2011, 2011, 548982. [Google Scholar] [CrossRef]
  22. Zhao, Z.; Li, C. Fractional difference/finite element approximations for the time-space fractional telegraph equation. Appl. Math. Comput. 2012, 219, 2975–2988. [Google Scholar] [CrossRef]
  23. Srivastava, V.K.; Awasthi, M.K.; Tamsir, M. RDTM solution of Caputo time fractional-order hyperbolic telegraph equation. AIP Adv. 2013, 3, 032142. [Google Scholar] [CrossRef]
  24. Kumar, S. A new analytical modelling for fractional telegraph equation via Laplace transform. Appl. Math. Model. 2014, 38, 3154–3163. [Google Scholar] [CrossRef]
  25. Joice Nirmala, R.; Balachandran, K. Analysis of solutions of time fractional telegraph equation. J. Korean Soc. Ind. Appl. Math. 2014, 18, 209–224. [Google Scholar] [CrossRef]
  26. Saadatmandi, A.; Mohabbati, M. Numerical solution of fractional telegraph equation via the tau method. Math. Rep. 2015, 17, 155–166. [Google Scholar]
  27. Li, H. A new analytical modelling for fractional telegraph equation via Elzaki transform. J. Adv. Math. 2015, 11, 5617–5625. [Google Scholar] [CrossRef]
  28. Dhunde, R.R.; Waghmare, G.L. Double Laplace transform method for solving space and time fractional telegraph equations. Int. J. Math. Math. Sci. 2016, 2016, 1414595. [Google Scholar] [CrossRef]
  29. Eltayeb, H.; Abdalla, Y.T.; Bachar, I.; Khabir, M.H. Fractional telegraph equation and its solution by natural transform decomposition method. Symmetry 2019, 11, 334. [Google Scholar] [CrossRef]
  30. Khan, H.; Shah, R.; Kumam, P.; Baleanu, D.; Arif, M. An efficient analytical technique, for the solution of fractional-order telegraph equations. Mathematics 2019, 7, 426. [Google Scholar] [CrossRef]
  31. Abdel-Rehim, E.A.; El-Sayed, A.M.A.; Hashem, A.S. Simulation of the approximate solutions of the time-fractional multi-term wave equations. Comput. Math. Appl. 2017, 73, 1134–1154. [Google Scholar] [CrossRef]
  32. Abdel-Rehim, E.A.; Hashem, A.S. Simulation of the Space-Time-Fractional Ultrasound Waves with Attenuation in Fractal Media. In Fractional Calculus, Proceedings of the ICFDA 2018, Amman, Jordan, 16–18 July 2019; Springer: Singapore, 2019; pp. 173–197. [Google Scholar]
  33. Abdel-Rehim, E.A. The approximate and analytic solutions of the time-fractional intermediate diffusion wave equation associated with the fokker–planck operator and applications. Axioms 2021, 10, 230. [Google Scholar] [CrossRef]
  34. Chu, H.-P.; Chen, C.-L. Hybrid differential transform and finite difference method to solve the nonlinear heat conduction problem. Commun. Nonlinear Sci. Numer. Simul. 2008, 13, 1605–1614. [Google Scholar] [CrossRef]
  35. Süngü, I.Ç.; Demir, H. Application of the hybrid differential transform method to the nonlinear equations. Sci. Res. 2012, 3, 246–250. [Google Scholar] [CrossRef]
  36. Arsalan, D. The numerical study of a hybrid method for solving telegraph equation. Appl. Math. Nonlinear Sci. 2020, 5, 293–302. [Google Scholar] [CrossRef]
  37. Du, Y.; Qin, B.; Zhao, C.; Zhu, Y.; Cao, J.; Ji, Y. A novel spatio-temporal synchronization method of roadside asynchronous MMW radar-camera for sensor fusion. IEEE Trans. Intell. Transp. Syst. 2021, 23, 22278–22289. [Google Scholar] [CrossRef]
  38. Odibat, Z.; Shawagfeh, N. Generalized Taylor’s formula. Appl. Math. Comput. 2007, 186, 286–293. [Google Scholar] [CrossRef]
  39. Erturk, V.S.; Momani, S.; Odibat, Z. Application of generalized differential transform method to multi-order fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2008, 13, 1642–1654. [Google Scholar] [CrossRef]
  40. Odibat, Z.; Kumar, S.; Shawagfeh, N.; Alsaedi, A.; Hayat, T. A study on the convergence conditions of generalized differential transform method. Math. Methods Appl. Sci. 2016, 40, 40–48. [Google Scholar] [CrossRef]
Figure 1. Comparison between the exact solution of Example 1 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Figure 1. Comparison between the exact solution of Example 1 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Symmetry 15 01759 g001
Figure 2. Comparison between the exact solution of Example 2 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Figure 2. Comparison between the exact solution of Example 2 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Symmetry 15 01759 g002
Figure 3. Comparison between the exact solution of Example 3 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Figure 3. Comparison between the exact solution of Example 3 for γ = 1 and numerical solutions for γ = 0.5 , 0.7 , 0.8 , 0.9 .
Symmetry 15 01759 g003
Figure 4. Solution of Example 4 for N = 40 , m = 20 .
Figure 4. Solution of Example 4 for N = 40 , m = 20 .
Symmetry 15 01759 g004
Figure 5. Solution of Example 4 for N = 80 , m = 20 .
Figure 5. Solution of Example 4 for N = 80 , m = 20 .
Symmetry 15 01759 g005
Table 1. The MAE for Example 1 at t = 0.001 , for γ = 1 , m = 3 , and different values of N.
Table 1. The MAE for Example 1 at t = 0.001 , for γ = 1 , m = 3 , and different values of N.
N MAE of FD - FDTM ( γ = 1 ) ROC CPU Time MAE of 2 D - DTM ( γ = 1 ) CPU Time
for α = 1 for α = 1
10 1 × 10 9 0.11 3.3 × 10 2 0.09
20 2.7 × 10 10 1.89 0.11 4.1 × 10 2 0.09
40 6.7 × 10 11 2.01 0.11 4.6 × 10 2 0.09
80 1.5 × 10 11 2.16 0.11 4.8 × 10 2 0.11
160 2.6 ×   10 12 2.51 0.16 5.01 × 10 2 0.14
Table 2. The approximate solution of Example 1 obtained with FD-FDTM and 2D-FDTM at t = 0.01 , for γ = 0.75 , α = 0.25 , m = 10 , and N = 10 .
Table 2. The approximate solution of Example 1 obtained with FD-FDTM and 2D-FDTM at t = 0.01 , for γ = 0.75 , α = 0.25 , m = 10 , and N = 10 .
x i FD-FDTM2D-FDTM Exact   Solution
γ = 0.75 γ = 0.75 γ = 1
0.1 1.08391073378 1.08391006915 1.08328706767
0.2 1.19790662077 1.19790588623 1.19721736312
0.3 1.32389155984 1.32389074805 1.32312981233
0.4 1.46312645063 1.46312555346 1.46228458943
0.5 1.617004802 1.61700381117 1.61607440219
0.6 1.7870666823 1.7870655864 1.78603843075
0.7 1.9750141259 1.9750129144 1.9738777322
0.8 2.1827281748 2.1827268341 2.1814722654
0.9 2.41228770087 2.412286213 2.4108997064
Table 3. The MAE for Example 2 at t = 0.001 , for γ = 1 , m = 3 , and different values of N.
Table 3. The MAE for Example 2 at t = 0.001 , for γ = 1 , m = 3 , and different values of N.
N MAE of FD - FDTM ROC CPU Time MAE of 2 D - FDTM CPU Time
for α = 1 for α = 1
10 4.2 × 10 10 0.11 9.1 × 10 1 0.14
20 1.1 × 10 10 1.8 0.12 9.5 × 10 1 0.12
40 2.8 × 10 11 1.97 0.11 9.7 × 10 1 0.11
80 6.7 × 10 12 2.06 0.14 9.8 × 10 1 0.11
160 1.09 × 10 12 2.61 0.17 9.9 × 10 1 0.17
Table 4. The approximate solution for Example 2 at t = 0.001 , for γ = 0.6 , α = 0.1 , m = 15 , and values of N = 10 .
Table 4. The approximate solution for Example 2 at t = 0.001 , for γ = 0.6 , α = 0.1 , m = 15 , and values of N = 10 .
x i FD-FDTM2D-FDTM Exact   Solution
γ = 0.6 γ = 0.6 γ = 1
0.1 0.097705 0.0164227 0.099966
0.2 0.19633389 0.0330099 0.2009337
0.3 0.297038 0.049763 0.30391186
0.4 0.40066 0.0666829 0.409931642
0.5 0.508292 0.083771 0.52005415
0.6 0.621011 0.101028 0.63538154
0.7 0.739945 0.118456 0.757068
0.8 0.866285 0.136056 0.8863315
0.9 1.0013 0.153829 1.024465743
Table 5. The MAE for Example 3 at t = 0.01 , for γ = 1 , m = 3 , and different values of N.
Table 5. The MAE for Example 3 at t = 0.01 , for γ = 1 , m = 3 , and different values of N.
N MAE of FD - FDTM CPU Time MAE of 2 D - DTM CPU Time
for α = 1 for α = 1
10 7.04 × 10 9 0.11 3.1 × 10 1 0.11
20 8.64 × 10 9 0.12 4.1 × 10 2 0.11
40 9.3 × 10 9 0.12 4.62 × 10 2 0.11
80 9.59 × 10 9 0.16 4.88 × 10 2 0.11
160 9.72 × 10 9 0.17 5.01 × 10 2 0.12
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

Sahraee, Z.; Arabameri, M. A Semi-Discretization Method Based on Finite Difference and Differential Transform Methods to Solve the Time-Fractional Telegraph Equation. Symmetry 2023, 15, 1759. https://doi.org/10.3390/sym15091759

AMA Style

Sahraee Z, Arabameri M. A Semi-Discretization Method Based on Finite Difference and Differential Transform Methods to Solve the Time-Fractional Telegraph Equation. Symmetry. 2023; 15(9):1759. https://doi.org/10.3390/sym15091759

Chicago/Turabian Style

Sahraee, Zahra, and Maryam Arabameri. 2023. "A Semi-Discretization Method Based on Finite Difference and Differential Transform Methods to Solve the Time-Fractional Telegraph Equation" Symmetry 15, no. 9: 1759. https://doi.org/10.3390/sym15091759

APA Style

Sahraee, Z., & Arabameri, M. (2023). A Semi-Discretization Method Based on Finite Difference and Differential Transform Methods to Solve the Time-Fractional Telegraph Equation. Symmetry, 15(9), 1759. https://doi.org/10.3390/sym15091759

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