On the concept of a conformable fractional differential equation
Soumen Shaw^{1} , Mohamed I. A. Othman^{2}
^{1}Department of mathematics, Indian Institute of Engineering Science and Technology, Shibpur, Howrah, 711103, India
^{2}Department of Mathematics, Faculty of Science, Zagazig University, P.O. Box 44519, Zagazig, Egypt
^{2}Corresponding author
Journal of Engineering and Thermal Sciences, Vol. 1, Issue 1, 2021, p. 1729.
https://doi.org/10.21595/jets.2021.22072
Received 21 May 2021; received in revised form 15 June 2021; accepted 29 June 2021; published 30 June 2021
JVE Conferences
A new simple wellbehaved definition of the fractional derivative termed as conformable fractional derivative and introducing a geometrical approach of fractional derivatives, nonintegral order initial value problems are an attempt to solve in this article. Based on the geometrical interpretation of the fractional derivatives, the solution curve is approximated numerically. Two special phenomena are employed for concave upward and downward curves. In order to obtain the solution of fractional order differential equation (FDE) with the integerorder initial condition, some new criteria on fractional derivatives are proposed.
Keywords: conformable differentiation, nonintegral differential equation, geometrical interpretation, fractional derivatives, numerical approximation.
1. Introduction
The impact of this fractional calculus in both pure and applied branches of science and engineering started to increase substantially during the last two decades apparently. Differential equations governing most of the physical system, evade closedform solutions and thereby necessitating the adoption of numerical methods to arrive at the desired solutions of the equations. Moreover, in many practical situations, some of the coefficients or functions in the differential equation may be nonlinear or are presented as a set of discrete data, and then, numerical methods become inevitable for obtaining solutions.
Since the last decade, fractionalorder differential equations have gained considerably more attention due to their applications in many engineering and scientific disciplines. As the mathematical models for the systems and processes, the fractional differential equation has been successfully applied in various fields of physics and engineering such as biophysics, bioengineering, quantum mechanics, finance, control theory, image and signal processing. A rather detailed account of diverse recent theoretical advances and applications of fractional calculus in the various fields can be found in the books of Sabatier et al. [14], Hilfer [6] and Atanackovic et al. [1]. Works such as [2, 3, 7, 8, 11, 15, 16, 2229] and the monographs [4, 13] analyzed qualitative and quantitative aspects of the solution of the fractionalorder differential equation. The methods employed in the aforementioned literature include the sequential technique of successive approximation as well as the classical fixedpoint approaches of Banach and Schauder. Othman et al. [20] studied the effect of hydrostatic initial stress and gravity on a fiberreinforced thermoelastic medium with a fractional derivative heat transfer. The effect of hydrostatic initial stress on the plane waves in a fiberreinforced magnetothermoelastic medium with fractional derivative heat transfer was explained by Sarkar et al. [21].
In the sequel, we would confine our attention to developing an algorithm for solving the initialvalue problem, namely:
Eq. (1) may be viewed as a curve in the $xy$plane, at each point on the curve the value of its $\alpha $order derivative is given in terms of $x$ and $y$, however, refers to a particular curve that passes through a given point $({x}_{0},{y}_{0})$.
In an endeavour to incorporate the curvature of the curve to arrive at an approximate numerical solution of the initialvalue problem (1); one tries to explore the possibility of alternative approaches. This leads to the following two types of procedures:
(i) Semianalytical approach based on the power series expansion (see Refs. [5, 11, 15, 17]).
(ii) Numerical approach applying the geometrical interpretation of fractional derivatives. A new geometrical approach to obtain an approximate solution of FDEs, subject to integerorder initial conditions which are significant to describe most of the physical phenomena, is the main course of study of this paper.
2. Fractional Taylor series
An extensive introduction to the use and meaning of fractional derivatives in physical and biological systems can be found in the article Metzler and Klafter [9]. The fractional Taylor series is a generalization of the Taylor series of fractional derivatives. The fractional Taylor series at the point $y=x+\mathrm{\Delta}x$ is defined by Odibat and Shawagfeh [10]:
where $\mathrm{\Gamma}\left(x\right)$ is the gamma Function, and ${D}_{x}^{\alpha}$ is the Caputo fractional derivative of order $0<\alpha <1$ with base point $x$, which is defined by:
where ${F}^{0}$^{}is the usual first derivative, and the notation $x+$ in Eq. (2) indicates the limit as we approach $x$ from the right. Note that there are other definitions of fractional derivatives, but the fractional Taylor series is valid for the Caputo form. The main distinguishing feature of the Caputo fractional derivative is that, analogous to integer order derivatives, the Caputo fractional derivative of a constant is zero. This property is very critical for a fractional Taylor series. It is also a notable fact that the third term of Eq. (2) involves the $\alpha $ order fractional derivative of $\alpha $ fractional derivative, which is not identical to the $2\alpha $ order fractional derivative.
To ensure this fact we consider the function:
where $y=x+h.$
Now the $\alpha $ order derivative of Eq. (4) is a constant when $q=\alpha $, and the $\alpha $ fractional derivative of that constant is zero. Then the coefficients of the fractional Taylor series can be found in the usual manner, by repeated fractional differentiations. The traditional integerorder Taylor series can be recovered from Eq. (2) when $\alpha =1$, applying the wellknown property of the Gamma function $\mathrm{\Gamma}(n+1)=n!$.
The fractional Taylor series gives an extremely good estimation of noninteger powerlaw functions. We now adopt Eq. (4) for a nonlinear powerlaw function, to illustrate this behaviour.
The traditional integerorder Taylor series approximation for Eq. (4) with $q=2$, expanded about $x$ and truncated at the secondorder term, is:
If $p=f\left(x\right)\text{,}$$f\text{'}\left(x\right)=2{\left.(yx{)}^{1}\right}_{y=x}=0\text{,}$$f\text{'}\left(x\right)=2{\left.(yx{)}^{1}\right}_{y=x}=0\text{,}$ and $f\text{'}\text{'}\left(x\right)=2.1{\left.(yx{)}^{0}\right}_{y=x}=2$, Eq. (5) becomes:
Hence, the secondorder Taylor series approximation is exact for $f\left(y\right),$ because the order of nonlinearity of the function matches the order of the Taylor series approximation. If $q=3$, a thirdorder Taylor series would provide an exact approximation.
However, if $q$ is a noninteger real number, no finite integer order Taylor series can give an exact match between the value of the function and its Taylor series approximation.
Now let us examine the fractional Taylor series approximation of Eq. (4), when $q>0$ is some real number. The Caputo fractional differentiation of Eq. (3) is:
Thus, firstorder fractional Taylor series for (4) expanded about $x$ is exact when $q=\alpha :$ the first term is $p=f\left(x\right).$ For the second term, use Eq. (7) along with the fact that the Caputo fractional derivative of the constant $p$ is zero to write:
Then for any $y>x$ we have (since $q=\alpha \text{,}$ the exponent $q\alpha =0\text{,}$ so the term $(yx{)}^{q\alpha}=1$ for any $y>x$):
Then of course the limit $y\to x+$ is:
Thus, the second term of fractional Taylor series is given by:
Here we have used $q=\alpha $ and $\mathrm{\Gamma}\left(1\right)=1.$ Since ${D}_{x}^{\alpha}f\left(y\right)$ is a constant, the remaining higher order Caputo derivatives are all zero. Therefore, the twoterm fractional Taylor series approximation is exact:
This is a very significant result. It proves that if we match the order of the fractional Taylor series approximation to the exponent in the powerlaw function, then the twoterm fractional Taylor series approximation to that function is exact.
3. Modified RiemannLiouville fractional derivative
In this section, we illustrate the fractional derivatives based on fractional differences which are slightly different from the classical RiemannLiouville framework. The fractional calculus so obtained is quite parallel to the classical calculus, and it involves noncommutative derivatives, which seems to be quite consistent with noncommutative geometry. The purpose herein is to contribute some new results to this approach.
3.1. Fractional derivatives via fractional differences
Definition 1. (Fractional right derivative): Let $f:R\to R$ be a continuous (not necessarily differentiable) function and let $h>0$ denotes the constant discretization span. The forward difference operator ${F}_{w}\left(h\right)$ is defined by:
the fractional difference of order $\alpha (0\alpha \le 1)$, from the right, of $f\left(x\right)$ is defined by the expression:
and the corresponding fractional derivative on the right is given by:
Definition 2. (Fractional left derivative): The left hand (or backward) fractional difference of the function $f\left(x\right)$ of order $\alpha (0\alpha \le 1)$, is defined as:
and its fractional derivative on the left is given by:
These are local definitions as compared with the standard integral approach. These are very close to the standard definitions of derivative and as a direct result, the $\alpha $order derivative of a constant is zero.
When ${f}_{+}^{\alpha}\left(x\right)={f}_{}^{\alpha}\left(x\right)$ to shorten in writing we shall set ${f}_{+}^{\alpha}={f}_{}^{\alpha}={f}^{\alpha}=\frac{{d}^{\alpha}f\left(x\right)}{d{x}^{\alpha}}$. Here we are fully in the Leibnitz framework, i.e., to say both ${d}^{\alpha}f$ and $d{x}^{\alpha}$^{}denote the finite increments.
3.2. Some important properties
Here we can mention few important properties of the aforementioned fractional derivative:
(1) The fractional derivative, as so defined, is not commutative. Clearly one has:
The equality holds if ${f}^{\alpha}\left(0\right)={f}^{\beta}\left(0\right).$
(2) Let $f\left(x\right):R\to R$ be a continuous function such that $f\left(x\right)$ has fractional derivative of $k\alpha $, where $k$ is any integer and $0<\alpha \le 1$, then:
(3) Let us consider the compound function $f\left(u\right(x\left)\right)$. Assume that $f\left(u\right)$ is $\alpha $differentiable with respect to $u$ and $u\left(x\right)$ is differentiable with respect to $x$ then:
(4) Assume that both $f\left(u\right)$ as well as $u\left(x\right)$ are $\alpha $differentiable with respect to $u$ and $x$ respectively, then:
(5) One can extent the aforementioned fractional derivatives for negative order also i.e. of order $\alpha (\alpha 0)$, can be expressed in the following integral form:
For $\alpha >0$ one can set:
The difference between Eqs. (22) and (23) is that the second one involves the constant $f\left(0\right)$ whereas the first equation does not. We shall refer to this fractional derivative as the modified RiemannLiouville derivative.
Proof. proofs of the properties (1)(4) are left for the readers. For (5); take Laplace to transform with respect to x in both sides of the Eq. (23) and denoting $L\left[f\right(x);x\to s]=\stackrel{}{f}\left(s\right)$ yields:
$=\frac{s}{\mathrm{\Gamma}(1\alpha )}{\int}_{\xi =0}^{\infty}\left\{{\int}_{x=0}^{\infty}(x\xi {)}^{\alpha}{e}^{sx}dx\right\}\left(f\right(\xi )f(0\left)\right)d\xi $
$=\frac{s}{\mathrm{\Gamma}(1\alpha )}{\int}_{\xi =0}^{\infty}\left\{{e}^{\xi s}{\int}_{t=0}^{\infty}{t}^{\alpha}{e}^{st}dt\right\}\left(f\right(\xi )f(0\left)\right)d\xi $
$=\frac{s}{\mathrm{\Gamma}(1\alpha )}{\int}_{\xi =0}^{\infty}{e}^{\xi s}\frac{\mathrm{\Gamma}(1\alpha )}{{s}^{1\alpha}}\left(f\right(\xi )f(0\left)\right)d\xi $
$={s}^{\alpha}{\int}_{\xi =0}^{\infty}{e}^{\xi s}\left(f\right(\xi )f(0\left)\right)d\xi ={s}^{\alpha}L\left[f\right(x\left)\right]{s}^{\alpha 1}f\left(0\right).$
Now according to the Eq. (15) (in definition 1):
$=\sum _{k=0}^{\infty}(1{)}^{k}{(}_{k}^{\alpha}){\int}_{0}^{\infty}f(x+(\alpha k\left)h\right){e}^{sx}dx$
$=\sum _{k=0}^{\infty}(1{)}^{k}{(}_{k}^{\alpha}\left){e}^{(\alpha k)hs}L\right[f\left(x\right)]=({e}^{hs}1{)}^{\alpha}\stackrel{}{f}\left(s\right).$
Therefore, from Eq. (15) we get:
Hence the proof.
3.3. Fractional Taylor series with MittagLeffler function
Let us consider the continuous function $f:R\to R$, $x\to f\left(x\right)$ has a fractional derivative of order $k\alpha $, for any positive integer $k$ and $\alpha $, $0<\alpha \le 1$, then the following relation holds:
$f(x+h)=\sum _{k=0}^{\infty}\frac{{h}^{k\alpha}}{k\alpha}({D}^{\alpha}{)}^{k}f\left(x\right),0\alpha \le 1,$
where ($\left({D}^{\alpha}{)}^{k}f\right(x)={D}^{\alpha}{D}^{\alpha}\cdots {D}^{\alpha}f(x)$ ($k$ times), is the derivative of order $\alpha +\alpha +\cdots +\alpha $ ($k$ times) of $f\left(x\right)$, and with the relation:
The definition of forwarding difference operator, given by the Eq. (13) and from the Eq. (14) one can show that the forward operator ${F}_{w}\left(h\right)$ satisfies the following fractional differential equation:
and the solution of this equation is ${F}_{w}\left(h\right)={E}_{\alpha}\left({h}^{\alpha}{D}_{x}^{\alpha}\right)=\sum _{k=0}^{\infty}\frac{{h}^{k\alpha}}{\left(k\alpha \right)!}.$
Therefore, the series can be expressed in the following manner:
In which ${D}_{x}$_{}is the differential operator with respect to $x$, and ${E}_{\alpha}\left(u\right)$ is the MittagLaffler function given by:
This fractional Taylor’s series does not hold with the standard RiemannLiouville derivatives, and it applies to nondifferentiable functions only. In addition, it is different from Osler’s series of fractional order [12].
4. The concept of approximation
In this section, we attempt to prove a special nature of the straight “line” (or tangent line) drawn at a point to the curve taking different order fractional derivatives as its slope at that point. To do that, firstly we demonstrate few important properties as follows:
Theorem 1. For any two real numbers $\alpha $ and $\beta (0\beta \alpha 1)$, ${D}^{\alpha}f>$ or ${<D}^{\beta}f$ accordingly $f\left(x\right)$ is a concave upward or concave downward function.
Proof. From the definition 1, for any real number $\alpha (0\alpha 1)$:
$=\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{\sum _{k=0}^{\infty}(1{)}^{k}\left({\mathrm{}}_{+}^{\alpha}\right)\{f(x+(\alpha k\left)h\right)f\left(\right(\alpha k\left)h\right)\}}{{h}^{\alpha}}={\sum}_{k=0}^{\infty}(1{)}^{k}u(x;k,\alpha ),$
where $u(x;k,\alpha )=\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{\left({\mathrm{}}_{+}^{\alpha}\right)\left\{f\right(x+(\alpha k)h)f((\alpha k)h\left)\right\}}{{h}^{\alpha}}.$
Now it is clear that, if $f\left(x\right)$ is a nondecreasing function then $u(x;k,\alpha )$ is also nondecreasing for all values for $x>0(0\alpha 1).$
Therefore:
$=\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{\left(\begin{array}{l}\sum _{k=0}^{\infty}(1{)}^{k}[\{{h}^{\alpha \beta}{(}_{k}^{\alpha}\left)f\right(x+(\alpha k)h){(}_{k}^{\beta})f(x+(\beta k\left)h\right)\}\\ \left\{{h}^{\alpha \beta}{(}_{k}^{\alpha}\right)f\left(\right(\alpha k\left)h\right){(}_{k}^{\beta}\left)f\right((\beta k)h\left)\right\}]\end{array}\right)}{{h}^{\alpha}}.$
Since, $x+(\alpha k)h>(\alpha k)h$ and $x+(\beta k)h>(\beta k)h$, then for a concave upward function $f\left(x\right)\text{,}$$\left\{{h}^{\alpha \beta}{(}_{k}^{\alpha}\right)f(x+(\alpha k\left)h\right){(}_{k}^{\beta}\left)f\right(x+(\beta k)h\left)\right\}>\left\{{h}^{\alpha \beta}{(}_{k}^{\alpha}\right)f\left(\right(\alpha k\left)h\right){(}_{k}^{\beta}\left)f\right((\beta k)h\left)\right\}$ for all positive real values of $x.$
Thus, if $f\left(x\right)$ is a concave upward function, for any two real numbers $\alpha $ and $\beta (0\beta \alpha 1)$, ${D}^{\alpha}f>{D}^{\beta}f$.
With the similar arguments, it can be prove that for a concave downward function $f\left(x\right)$, ${D}^{\alpha}f<{D}^{\beta}f$, where $(0<\beta <\alpha <1).$
Hence the Proof.
Therefore, in order to estimate the curve, for a concave upward function, the “lines” (or tangent lines) with the highest order fractional derivative ${D}_{x}^{\alpha}f\left(x\right)$ where $\alpha \in \left(\mathrm{0,1}\right]$ to the curve at the point will be a more appropriate one.
Whilst, for a concave downward curve, the “line” (or tangent) at a point having the slope ${D}_{x}^{\alpha}f\left(x\right)$ with the lowest order fractional derivative where $\alpha \in \left(\mathrm{0,1}\right]$ at that point to the curve would have been more suitable than the conventional tangent line.
5. Graphical illustration
Here we illustrate the nature of the tangent lines with different fractionalorder derivatives ${D}_{x}^{\alpha}f\left(x\right)$ where $\alpha \in \left(\mathrm{0,1}\right]$, which is taken as the slope of the line, for concave upward and downward curves. For this purpose, we have considered two simple polynomials (see Fig. 1 and Fig. 2).
Let $y\left(x\right)$ be an approximation in the interval $({x}_{0},{y}_{0})$. Expanding $y\left(x\right)$ about the point $x={x}_{m}$, we have:
$+\frac{{D}_{x}^{\alpha}{)}^{n+1}f\left(\xi \right)}{\left(\right(n+1)\alpha +1)!}(x{x}_{m}{)}^{(n+1)\alpha},$
where ${D}_{x}^{\alpha}y\left({x}_{m}\right)$ is the $\alpha $ order derivative of $y\left(x\right)$ at $x={x}_{m}$_{}and ${x}_{0}\le \xi \le {x}_{m}.$
Fig. 1. Various “lines” (or tangent lines) with different slopes for $f\left(x\right)={x}^{3}$
Fig. 2. Various “lines” (or tangent lines) with different slopes for $f\left(x\right)={x}^{3}$
Replace $x$ by ${x}_{m}+h$ in Eq. (35) we obtain:
Assuming the variation of ${D}_{x}^{\alpha}{D}_{x}^{\alpha}y$ is negligible in ${x}_{m}\le x\le {x}_{m+1}$ the truncation error can be put as ${e}_{T}\approx k{h}^{2\alpha}$.
Hence, to determine the function $y\left(x\right)$, which will be the solution of the initial value problem Eq. (1), may apply the following algorithm:
where $\alpha $ is a constant to be chosen in such a way that $\alpha =\mathrm{m}\mathrm{a}\mathrm{x}$ or $\mathrm{m}\mathrm{i}\mathrm{n}\{x,x\in (\mathrm{0,1}\left]\right\}$ accordingly the curve is concave upward or downward i.e., $\frac{{d}^{2}y}{d{x}^{2}}>$or $<0$ at $({x}_{m},{y}_{m})$.
6. Convergence criteria
A method is convergent if, for every ODE with a Lipschitz function $f$ and every fixed $T$, with $T=n.h$, it holds that:
Global error: ${\epsilon}_{n}={y}_{n}y\left({x}_{n}\right)$ and ${\epsilon}_{n+1}={y}_{n+1}y\left({x}_{n+1}\right).$
Local error: ${l}_{n+1}={\stackrel{~}{y}}_{n+1}y\left({x}_{n+1}\right)$ where ${\stackrel{~}{y}}_{n+1}=y\left({x}_{n}\right)+\frac{{h}^{\alpha}}{\alpha !}f({x}_{n},y({x}_{n}\left)\right).$
Therefore, local error definition implies the residual:
Now comparing with the Taylor series expansion Eq. (28), the local error for this proposed method is ${l}_{n+1}\approx \frac{{h}^{2\alpha}}{\left(2\alpha \right)!}{D}_{x}^{\alpha}{D}_{x}^{\alpha}y\left({x}_{n}\right)$.
Thus, the global error recursion is given by:
Take norms and use Lipschitz condition:
Lemma. If $\left\{{a}_{n}\right\},{a}_{0}=0,$ be a sequence of nonnegative numbers satisfying $\left\{{a}_{n}\right\}$, ${a}_{0}=0$, ${a}_{n+1}\le (1+h\mu ){a}_{n}+c{h}^{2}$, for $\mu \ge 0$, then ${a}_{n}\le \frac{c}{\mu}h\left[\right(1+h\mu {)}^{n}{a}_{n}1]$, $n=\mathrm{1,2},\dots $.
Proof. Standard result of sequence of real nos.
Theorem 2. Prove that the global error $\Vert {\epsilon}_{m,h}\Vert \to 0$ as $h$ approaches to zero.
Proof. Assuming that the function $f$ is sufficiently differentiable, given $h$ and a fixed $T=n.h$, let us consider ${\epsilon}_{m,h}={y}_{m,h}t\left({x}_{m}\right)$.
Applying the lemma to global error recursion, we get:
with:
Since:
We have, for $m\le n$:
Therefore:
which proves the convergence as:
Hence, the iteration converges invariably in keeping with expectation.
7. Applications
In order to assess the applicability of the proposed method, we have considered a number of fractional differential equations with initial conditions. The iteration procedure is carried out through computer simulation (MatLab 2015a) with a moderate value of the step length.
For the computational purpose we consider the following fractional order initial value problems.
Example 1. (Example – 4.1 [18]):
Example 2. (Example – 4.3 [18]) with $y\left(0\right)=0$:
Example 3. (Eigen value problem: example – 5.4 [19]):
Fig. 3. Comparison of solutions between the analytical result and present analysis (Ex. 1)
For numerical computation, in all of the above examples, we consider $h=$ 0.001, ${\alpha}_{\mathrm{m}\mathrm{i}\mathrm{n}}=$ 0.125 and ${\alpha}_{\mathrm{m}\mathrm{a}\mathrm{x}}=$ 1. The analytical solutions for the considered initial value problems (IVPs) are presented in Table 1 (for details see Refs. [18, 19]). Those results are compared with the approximate solutions obtained by using the proposed method. It has been observed that the present approximations (by using the iteration scheme Eq. (37)) are in high agreement with the results obtained by [18]. Figs. 35 represent the solution curves of the IVPs Eqs. (48)(50) by using the present analysis and analytical results. It has also been noticed that this newly adopted method is identical with the analytical solutions up to two decimal places.
Fig. 4. Comparison of solutions between the analytical result and present analysis (Ex. 2)
Fig. 5. Comparison of solutions between the analytical result and present analysis (Ex. 3)
Table 1. Analytical solutions of the IBPs
Example

IVP

Solution

Ex. 1

${D}^{1/2}y+y={x}^{2}+2{x}^{3/2},y\left(0\right)=0$

$y={x}^{2}$

Ex. 2

${D}^{1/2}y=x{e}^{x}\sqrt{xy},y\left(0\right)=0$

$y=\frac{2}{3}{x}^{2/3}\mathrm{e}\mathrm{x}\mathrm{p}(x)$

Ex. 3

$\left({T}_{\alpha}^{a}y\right)\left(t\right)=\lambda y\left(t\right),y\left(a=0\right)=1,t0,\lambda =0.5,\alpha =0.5$

$y={y}_{0}{{e}^{\lambda (ta)}}^{{}^{\alpha}/\alpha}$

8. Conclusions
With the help of the geometrical interpretation of fractional order derivatives for different values of the parameter $\alpha \in \left(\mathrm{0,1}\right]$ at a fixed point on a curve, in this paper, we developed a new numerical iterative method that can give an approximate solution curve of the fractional initial value problems. Here we use two different values of fractional derivatives to approximate the solution curve for the concave and convex parts. Purely depending upon geometrical interpretation, this type of works have not yet been reported on a fractional differential equation.
During the decade, it has been established that several biological and physical aspects could be encompassed in the framework of a suitable fractional differential equation. We assert that the approach to the fractional calculus via fractional difference in combination with fractional Taylor’s series provides a method of approximation of the solution curve for IVPs in Leibniz’s sense. The proposed method could be fairly competent in this kind of problems.
References
 T. M. Atanackovic, S. Pilipovic, B. Stankovic, and D. Zorica, Fractional Calculus with Application in Mechanics. London: Wiley, 2014. [Publisher]
 K. Diethelm and N. J. Ford, “Analysis of fractional differential equations,” Journal of Mathematical Analysis and Applications, Vol. 265, No. 2, pp. 229–248, Jan. 2002, https://doi.org/10.1006/jmaa.2000.7194 [Publisher]
 K. Diethelm and N. J. Ford, “Multiorder fractional differential equations and their numerical solution,” Applied Mathematics and Computation, Vol. 154, No. 3, pp. 621–640, Jul. 2004, https://doi.org/10.1016/s00963003(03)007392 [Publisher]
 K. Diethelm, The Analysis of Fractional Differential Equations. Heidelberg: Springer, 2010. [Publisher]
 A. ElAjou, O. Arqub, Z. Zhour, and S. Momani, “New Results on Fractional Power Series: Theories and Applications,” Entropy, Vol. 15, No. 12, pp. 5305–5323, Dec. 2013, https://doi.org/10.3390/e15125305 [Publisher]
 R. Hilfer, Application of Fractional Calculus to Physics. Singapore: World Scientific, 2000. [Publisher]
 V. Lakshmikantham and A. S. Vatsala, “Basic theory of fractional differential equations,” Nonlinear Analysis: Theory, Methods and Applications, Vol. 69, No. 8, pp. 2677–2682, Oct. 2008, https://doi.org/10.1016/j.na.2007.08.042 [Publisher]
 V. Lakshmikantham and A. S. Vatsala, “General uniqueness and monotone iterative technique for fractional differential equations,” Applied Mathematics Letters, Vol. 21, No. 8, pp. 828–834, Aug. 2008, https://doi.org/10.1016/j.aml.2007.09.006 [Publisher]
 R. Metzler and J. Klafter, “The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics,” Journal of Physics A: Mathematical and General, Vol. 37, No. 31, pp. R161–R208, Aug. 2004, https://doi.org/10.1088/03054470/37/31/r01 [Publisher]
 Z. M. Odibat and N. T. Shawagfeh, “Generalized Taylor’s formula,” Applied Mathematics and Computation, Vol. 186, No. 1, pp. 286–293, Mar. 2007, https://doi.org/10.1016/j.amc.2006.07.102 [Publisher]
 Z. Odibat, S. Momani, and V. S. Erturk, “Generalized differential transform method: Application to differential equations of fractional order,” Applied Mathematics and Computation, Vol. 197, No. 2, pp. 467–477, Apr. 2008, https://doi.org/10.1016/j.amc.2007.07.068 [Publisher]
 T. J. Osler, “Taylor’s series generalized for fractional derivatives and applications,” SIAM Journal on Mathematical Analysis, Vol. 2, No. 1, pp. 37–48, Feb. 1971, https://doi.org/10.1137/0502004 [Publisher]
 Igor Podlubny, Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. San Diego: Academic Press, 1999. [Search CrossRef]
 J. Sabatier, O. P. Agrawal, and J. A. Tenreiro Machado, Advances in Fractional Calculus.Theoretical Developments and Applications in Physics and Engineering. Dordrecht: Springer, 2007. [Publisher]
 C. C. Tisdell, “On the application of sequential and fixedpoint methods to fractional differential equations of arbitrary order,” Journal of Integral Equations and Applications, Vol. 24, No. 2, pp. 283–319, Jun. 2012, https://doi.org/10.1216/jie2012242283 [Publisher]
 C. C. Tisdell, “Solutions to fractional differential equations that extend,” Journal of Classical Analysis, Vol. 5, No. 2, pp. 129–136, 2014, https://doi.org/10.7153/jca0511 [Publisher]
 X.F. Zhou, F. Yang, and W. Jiang, “Analytic study on linear neutral fractional differential equations,” Applied Mathematics and Computation, Vol. 257, pp. 295–307, Apr. 2015, https://doi.org/10.1016/j.amc.2014.12.056 [Publisher]
 R. Khalil, M. Al Horani, A. Yousef, and M. Sababheh, “A new definition of fractional derivative,” Journal of Computational and Applied Mathematics, Vol. 264, pp. 65–70, Jul. 2014, https://doi.org/10.1016/j.cam.2014.01.002 [Publisher]
 T. Abdeljawad, “On conformable fractional calculus,” Journal of Computational and Applied Mathematics, Vol. 279, pp. 57–66, May 2015, https://doi.org/10.1016/j.cam.2014.10.016 [Publisher]
 M. Ibrahim Othman, S. M. Said, and N. Sarker, “Effect of hydrostatic initial stress on a fiberreinforced thermoelastic medium with fractional derivative heat transfer,” Multidiscipline Modeling in Materials and Structures, Vol. 9, No. 3, pp. 410–426, Sep. 2013, https://doi.org/10.1108/mmms1120120026 [Publisher]
 N. Sarkar, S. Y. Atwa, and M. I. A. Othman, “The Effect of Hydrostatic Initial Stress on the Plane Waves in a FiberReinforced MagnetoThermoelastic Medium with Fractional Derivative Heat Transfer,” International Applied Mechanics, Vol. 52, No. 2, pp. 203–216, Mar. 2016, https://doi.org/10.1007/s1077801607484 [Publisher]
 O. Brandibur and E. Kaslik, “Stability analysis of multiterm fractionaldifferential equations with three fractional derivatives,” Journal of Mathematical Analysis and Applications, Vol. 495, No. 2, p. 124751, Mar. 2021, https://doi.org/10.1016/j.jmaa.2020.124751 [Publisher]
 N. D. Cong, H. T. Tuan, and H. Trinh, “On asymptotic properties of solutions to fractional differential equations,” Journal of Mathematical Analysis and Applications, Vol. 484, No. 2, p. 123759, Apr. 2020, https://doi.org/10.1016/j.jmaa.2019.123759 [Publisher]
 F. Alzahrani and I. A. Abbas, “Generalized Thermoelastic Interactions in a Poroelastic Material Without Energy Dissipations,” International Journal of Thermophysics, Vol. 41, No. 7, Jul. 2020, https://doi.org/10.1007/s10765020026730 [Publisher]
 T. Saeed, I. Abbas, and M. Marin, “A GL Model on ThermoElastic Interaction in a Poroelastic Material Using Finite Element Method,” Symmetry, Vol. 12, No. 3, p. 488, Mar. 2020, https://doi.org/10.3390/sym12030488 [Publisher]
 I. Abbas, “Natural frequencies of a poroelastic hollow cylinder,” Acta Mechanica, Vol. 186, No. 14, pp. 229–237, Oct. 2006, https://doi.org/10.1007/s007070060314y [Publisher]
 A. M. ElNaggar, Z. Kishka, A. M. AbdAlla, I. A. Abbas, S. M. AboDahab, and M. Elsagheer, “On the Initial Stress, Magnetic Field, Voids and Rotation Effects on Plane Waves in Generalized Thermoelasticity,” Journal of Computational and Theoretical Nanoscience, Vol. 10, No. 6, pp. 1408–1417, Jun. 2013, https://doi.org/10.1166/jctn.2013.2862 [Publisher]
 G. Palani and I. A. Abbas, “Free Convection MHD Flow with Thermal Radiation from an ImpulsivelyStarted Vertical Plate,” Nonlinear Analysis: Modelling and Control, Vol. 14, No. 1, pp. 73–84, Jan. 2009, https://doi.org/10.15388/na.2009.14.1.14531 [Publisher]
 R. Kumar and I. Abbas, “Deformation Due to Thermal Source in Micropolar Thermoelastic Media with Thermal and Conductive Temperatures,” Journal of Computational and Theoretical Nanoscience, Vol. 10, No. 9, pp. 2241–2247, Sep. 2013, https://doi.org/10.1166/jctn.2013.3193 [Publisher]