# On the concept of a conformable fractional differential equation

## Soumen Shaw1, Mohamed I. A. Othman2

1Department of mathematics, Indian Institute of Engineering Science and Technology, Shibpur, Howrah, 711103, India

2Department of Mathematics, Faculty of Science, Zagazig University, P.O. Box 44519, Zagazig, Egypt

2Corresponding author

Journal of Engineering and Thermal Sciences, Vol. 1, Issue 1, 2021, p. 17-29. 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

Copyright © 2021 Soumen Shaw, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Editor's Pick
Views 115
CrossRef Citations 0
Abstract.

A new simple well-behaved definition of the fractional derivative termed as conformable fractional derivative and introducing a geometrical approach of fractional derivatives, non-integral 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 integer-order initial condition, some new criteria on fractional derivatives are proposed.

Keywords: conformable differentiation, non-integral 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 closed-form 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 non-linear or are presented as a set of discrete data, and then, numerical methods become inevitable for obtaining solutions.

Since the last decade, fractional-order 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, 22-29] and the monographs [4, 13] analyzed qualitative and quantitative aspects of the solution of the fractional-order differential equation. The methods employed in the aforementioned literature include the sequential technique of successive approximation as well as the classical fixed-point approaches of Banach and Schauder. Othman et al. [20] studied the effect of hydrostatic initial stress and gravity on a fiber-reinforced thermoelastic medium with a fractional derivative heat transfer. The effect of hydrostatic initial stress on the plane waves in a fiber-reinforced magneto-thermoelastic 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 initial-value problem, namely:

(1)

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 $\left({x}_{0},{y}_{0}\right)$.

In an endeavour to incorporate the curvature of the curve to arrive at an approximate numerical solution of the initial-value problem (1); one tries to explore the possibility of alternative approaches. This leads to the following two types of procedures:

(i) Semi-analytical 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 integer-order 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]:

(2)
$F\left(y\right)=F\left(x\right)+\frac{\left(y-x{\right)}^{\alpha }}{\mathrm{\Gamma }\left(\alpha +1\right)}{D}_{x}^{\alpha }F\left(x+\right)+\frac{\left(y-x{\right)}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}{D}_{x}^{\alpha }{D}_{x}^{\alpha }F\left(x+\right)+\cdots ,$

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:

(3)
${D}_{x}^{\alpha }F\left(y\right)=\frac{1}{\mathrm{\Gamma }\left(1-\alpha \right)}{\int }_{x}^{y}F\text{'}\left(y-x\right)\left(u-x{\right)}^{-\alpha }du,$

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:

(4)
$f\left(y\right)=p+\left(y-x{\right)}^{q},$

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 integer-order Taylor series can be recovered from Eq. (2) when $\alpha =1$, applying the well-known property of the Gamma function $\mathrm{\Gamma }\left(n+1\right)=n!$.

The fractional Taylor series gives an extremely good estimation of non-integer power-law functions. We now adopt Eq. (4) for a non-linear power-law function, to illustrate this behaviour.

The traditional integer-order Taylor series approximation for Eq. (4) with $q=2$, expanded about $x$ and truncated at the second-order term, is:

(5)
$f\left(y\right)=f\left(x\right)+f\text{'}\left(y-x\right)+f\text{'}\text{'}\frac{\left(y-x{\right)}^{2}}{2!}.$

If $p=f\left(x\right)\text{,}$$f\text{'}\left(x\right)=2{\left(y-x{\right)}^{1}|}_{y=x}=0\text{,}$$f\text{'}\left(x\right)=2{\left(y-x{\right)}^{1}|}_{y=x}=0\text{,}$ and $f\text{'}\text{'}\left(x\right)=2.1{\left(y-x{\right)}^{0}|}_{y=x}=2$, Eq. (5) becomes:

(6)
$f\left(y\right)=p+\left(y-x{\right)}^{2}.$

Hence, the second-order 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 third-order Taylor series would provide an exact approximation.

However, if $q$ is a non-integer 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:

(7)
${D}_{x}^{\alpha }f\left(y\right)={D}_{x}^{\alpha }\left[p+\left(y-x{\right)}^{q}\right]=\frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-1+\alpha \right)}\left(y-x{\right)}^{q-\alpha }.$

Thus, first-order 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:

(8)
${D}_{x}^{\alpha }f\left(y\right)=\frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-1+\alpha \right)}\left(y-x{\right)}^{q-\alpha }.$

Then for any $y>x$ we have (since $q=\alpha \text{,}$ the exponent $q-\alpha =0\text{,}$ so the term $\left(y-x{\right)}^{q-\alpha }=1$ for any $y>x$):

(9)
${D}_{x}^{\alpha }f\left(y\right)=\frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-1+\alpha \right)}.1.$

Then of course the limit $y\to x+$ is:

(10)
${D}_{x}^{\alpha }f\left(x+\right)=\frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-1+\alpha \right)}.1.$

Thus, the second term of fractional Taylor series is given by:

(11)
${D}_{x}^{\alpha }f\left(x+\right).\frac{\left(y-x{\right)}^{\alpha }}{\mathrm{\Gamma }\left(\alpha +1\right)}=\left[\frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-1+\alpha \right)}.1\right]\frac{\left(y-x{\right)}^{\alpha }}{\mathrm{\Gamma }\left(\alpha +1\right)}=\left(y-x{\right)}^{q}.$

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 two-term fractional Taylor series approximation is exact:

(12)
$f\left(y\right)=f\left(x\right)+\frac{\left(y-x{\right)}^{\alpha }}{\mathrm{\Gamma }\left(\alpha +1\right)}{D}_{x}^{\alpha }f\left(x+\right)+\frac{\left(y-x{\right)}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}{D}_{x}^{\alpha }{D}_{x}^{\alpha }f\left(x+\right)+\cdots =p+\left(y-x{\right)}^{q}+0+\cdots .$

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 power-law function, then the two-term fractional Taylor series approximation to that function is exact.

#### 3. Modified Riemann-Liouville fractional derivative

In this section, we illustrate the fractional derivatives based on fractional differences which are slightly different from the classical Riemann-Liouville framework. The fractional calculus so obtained is quite parallel to the classical calculus, and it involves non-commutative derivatives, which seems to be quite consistent with non-commutative 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:

(13)
${F}_{w}\left(h\right)f\left(x\right)=f\left(x+h\right),$

the fractional difference of order , from the right, of $f\left(x\right)$ is defined by the expression:

(14)
${\mathrm{\Delta }}_{+}^{\alpha }f\left(x\right)=\left({F}_{w}-1{\right)}^{\alpha }f\left(x\right)=\sum _{k=0}^{\infty }\left(-1{\right)}^{k}{\left(}_{k}^{\alpha }\right)f\left(x+\left(\alpha -k\right)h\right),$

and the corresponding fractional derivative on the right is given by:

(15)
${f}_{+}^{\alpha }\left(x\right)=\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{{\mathrm{\Delta }}_{+}^{\alpha }\left[f\left(x\right)-f\left(0\right)\right]}{{h}^{\alpha }}.$

Definition 2. (Fractional left derivative): The left hand (or backward) fractional difference of the function $f\left(x\right)$ of order , is defined as:

(16)
${\mathrm{\Delta }}_{-}^{\alpha }f\left(x\right)=\left(1-{F}_{w}{\right)}^{\alpha }f\left(x\right)=\sum _{k=0}^{\infty }\left(-1{\right)}^{k}{\left(}_{k}^{\alpha }\right)f\left(x-kh\right),$

and its fractional derivative on the left is given by:

(17)
${f}_{-}^{\alpha }\left(x\right)=\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}\frac{{\mathrm{\Delta }}_{-}^{\alpha }\left[f\left(x\right)-f\left(0\right)\right]}{{h}^{\alpha }}.$

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 . 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:

(18)
${D}^{\alpha }{D}^{\beta }f\left(x\right)\ne {D}^{\beta }{D}^{\alpha }f\left(x\right).$

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:

(19)

(3) Let us consider the compound function $f\left(u\left(x\right)\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:

(20)
$\left(f\left(u\left(x\right)\right){\right)}^{\alpha }={f}_{u}^{\left(\alpha \right)}\left(u\right)\left({u\text{'}}_{x}{\right)}^{\alpha }.$

(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:

(21)

(5) One can extent the aforementioned fractional derivatives for negative order also i.e. of order , can be expressed in the following integral form:

(22)

For $\alpha >0$ one can set:

(23)
(24)

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 Riemann-Liouville 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\left(x\right);x\to s\right]=\stackrel{-}{f}\left(s\right)$ yields:

(25)

Now according to the Eq. (15) (in definition 1):

(26)

Therefore, from Eq. (15) we get:

(27)

Hence the proof.

#### 3.3. Fractional Taylor series with Mittag-Leffler 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:

(28)

where ($\left({D}^{\alpha }{\right)}^{k}f\left(x\right)={D}^{\alpha }{D}^{\alpha }\cdots {D}^{\alpha }f\left(x\right)$ ($k$ times), is the derivative of order $\alpha +\alpha +\cdots +\alpha$ ($k$ times) of $f\left(x\right)$, and with the relation:

(29)
$\mathrm{\Gamma }\left(1+k\alpha \right)=\left(k\alpha \right)!.$

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:

(30)
$\left({D}_{h}{\right)}^{\alpha }{F}_{w}\left(h\right)=\left({D}_{x}{\right)}^{\alpha }{F}_{w}\left(h\right),$

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:

(31)
$f\left(x+h\right)={E}_{\alpha }\left({h}^{\alpha }{D}_{x}^{\alpha }\right)f\left(x\right).$

In which ${D}_{x}$is the differential operator with respect to $x$, and ${E}_{\alpha }\left(u\right)$ is the Mittag-Laffler function given by:

(32)

This fractional Taylor’s series does not hold with the standard Riemann-Liouville derivatives, and it applies to non-differentiable 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 , ${D}^{\alpha }f>$ or ${ accordingly $f\left(x\right)$ is a concave upward or concave downward function.

Proof. From the definition 1, for any real number :

(33)

where

Now it is clear that, if $f\left(x\right)$ is a non-decreasing function then $u\left(x;k,\alpha \right)$ is also non-decreasing for all values for

Therefore:

(34)

Since, $x+\left(\alpha -k\right)h>\left(\alpha -k\right)h$ and $x+\left(\beta -k\right)h>\left(\beta -k\right)h$, then for a concave upward function $f\left(x\right)\text{,}$$\left\{{h}^{\alpha -\beta }{\left(}_{k}^{\alpha }\right)f\left(x+\left(\alpha -k\right)h\right)-{\left(}_{k}^{\beta }\right)f\left(x+\left(\beta -k\right)h\right)\right\}>\left\{{h}^{\alpha -\beta }{\left(}_{k}^{\alpha }\right)f\left(\left(\alpha -k\right)h\right)-{\left(}_{k}^{\beta }\right)f\left(\left(\beta -k\right)h\right)\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 , ${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 $\left(0<\beta <\alpha <1\right).$

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(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(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 fractional-order derivatives ${D}_{x}^{\alpha }f\left(x\right)$ where $\alpha \in \left(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 $\left({x}_{0},{y}_{0}\right)$. Expanding $y\left(x\right)$ about the point $x={x}_{m}$, we have:

(35)
$y\left(x\right)=y\left({x}_{m}\right)+\frac{\left(x-{x}_{m}{\right)}^{\alpha }}{\alpha !}{D}_{x}^{\alpha }y\left({x}_{m}\right)+\frac{\left(x-{x}_{m}{\right)}^{2\alpha }}{\left(2\alpha \right)!}{D}_{x}^{\alpha }{D}_{x}^{\alpha }y\left({x}_{m}\right)+\cdots$

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:

(36)

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:

(37)

where $\alpha$ is a constant to be chosen in such a way that $\alpha =\mathrm{m}\mathrm{a}\mathrm{x}$ or accordingly the curve is concave upward or downward i.e., or $<0$ at $\left({x}_{m},{y}_{m}\right)$.

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

(38)
$\underset{n\to \infty }{\mathrm{l}\mathrm{i}\mathrm{m}}‖{y}_{n,h}-y\left(T\right)‖=0.$

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\left({x}_{n},y\left({x}_{n}\right)\right).$

Therefore, local error definition implies the residual:

(39)
$y\left({x}_{n+1}\right)=y\left({x}_{n}\right)+\frac{{h}^{\alpha }}{\alpha !}f\left({x}_{n},y\left({x}_{n}\right)\right)-{l}_{n+1}.$

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:

(40)
${\epsilon }_{n+1}={\epsilon }_{n}+\frac{{h}^{\alpha }}{\alpha !}\left[f\left({x}_{n},y\left({x}_{n}\right)\right)+{\epsilon }_{n}\right]-\frac{{h}^{\alpha }}{\alpha !}f\left({x}_{n},y\left({x}_{n}\right)\right)+{l}_{n+1}.$

Take norms and use Lipschitz condition:

(41)
$‖{\epsilon }_{n+1}‖\le ‖{\epsilon }_{n}‖+\frac{{h}^{\alpha }}{\alpha !}L\left[f\right].‖{\epsilon }_{n}‖+‖{l}_{n+1}‖.$

Lemma. If be a sequence of non-negative numbers satisfying $\left\{{a}_{n}\right\}$, ${a}_{0}=0$, ${a}_{n+1}\le \left(1+h\mu \right){a}_{n}+c{h}^{2}$, for $\mu \ge 0$, then ${a}_{n}\le \frac{c}{\mu }h\left[\left(1+h\mu {\right)}^{n}{a}_{n}-1\right]$, $n=1,2,\dots$.

Proof. Standard result of sequence of real nos.

Theorem 2. Prove that the global error $‖{\epsilon }_{m,h}‖\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:

(42)

with:

(43)
$c=\underset{m}{\mathrm{m}\mathrm{a}\mathrm{x}}‖{l}_{m}‖{h}^{2\alpha }\approx \underset{m}{\mathrm{m}\mathrm{a}\mathrm{x}}‖{D}_{x}^{\alpha }{D}_{x}^{\alpha }y‖/\left(2\alpha \right)!.$

Since:

(44)
${\left(1+\frac{{h}^{\alpha }}{\alpha !}L\left[f\right]\right)}^{m}\le {e}^{m\frac{{h}^{\alpha }}{\alpha !}L\left[f\right]}\le {e}^{TL\left[f\right]}.$

We have, for $m\le n$:

(45)
$‖{\epsilon }_{m,h}‖\le \frac{c}{L\left[f\right]}\frac{{h}^{\alpha }}{\alpha !}\left({e}^{TL\left[f\right]}-1\right).$

Therefore:

(46)
$‖{\epsilon }_{m,h}‖\le C\left(T\right).\frac{{h}^{\alpha }}{\alpha !},$

which proves the convergence as:

(47)
$\underset{h\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}‖{\epsilon }_{m,h}‖=0.$

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]):

(48)

Example 2. (Example – 4.3 [18]) with $y\left(0\right)=0$:

(49)

Example 3. (Eigen value problem: example – 5.4 [19]):

(50)

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. 3-5 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 $y={x}^{2}$ Ex. 2 $y=\frac{2}{3}{x}^{2/3}\mathrm{e}\mathrm{x}\mathrm{p}\left(-x\right)$ Ex. 3

#### 8. Conclusions

With the help of the geometrical interpretation of fractional order derivatives for different values of the parameter $\alpha \in \left(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

1. T. M. Atanackovic, S. Pilipovic, B. Stankovic, and D. Zorica, Fractional Calculus with Application in Mechanics. London: Wiley, 2014. [Publisher]
2. 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]
3. K. Diethelm and N. J. Ford, “Multi-order 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/s0096-3003(03)00739-2 [Publisher]
4. K. Diethelm, The Analysis of Fractional Differential Equations. Heidelberg: Springer, 2010. [Publisher]
5. A. El-Ajou, 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]
6. R. Hilfer, Application of Fractional Calculus to Physics. Singapore: World Scientific, 2000. [Publisher]
7. 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]
8. 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]
9. 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/0305-4470/37/31/r01 [Publisher]
10. 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]
11. 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]
12. 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]
13. 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]
14. 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]
15. C. C. Tisdell, “On the application of sequential and fixed-point 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/jie-2012-24-2-283 [Publisher]
16. 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/jca-05-11 [Publisher]
17. 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]
18. 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]
19. 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]
20. M. Ibrahim Othman, S. M. Said, and N. Sarker, “Effect of hydrostatic initial stress on a fiber-reinforced 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/mmms-11-2012-0026 [Publisher]
21. N. Sarkar, S. Y. Atwa, and M. I. A. Othman, “The Effect of Hydrostatic Initial Stress on the Plane Waves in a Fiber-Reinforced Magneto-Thermoelastic Medium with Fractional Derivative Heat Transfer,” International Applied Mechanics, Vol. 52, No. 2, pp. 203–216, Mar. 2016, https://doi.org/10.1007/s10778-016-0748-4 [Publisher]
22. O. Brandibur and E. Kaslik, “Stability analysis of multi-term fractional-differential 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]
23. 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]
24. 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/s10765-020-02673-0 [Publisher]
25. T. Saeed, I. Abbas, and M. Marin, “A GL Model on Thermo-Elastic 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]
26. I. Abbas, “Natural frequencies of a poroelastic hollow cylinder,” Acta Mechanica, Vol. 186, No. 1-4, pp. 229–237, Oct. 2006, https://doi.org/10.1007/s00707-006-0314-y [Publisher]
27. A. M. El-Naggar, Z. Kishka, A. M. Abd-Alla, I. A. Abbas, S. M. Abo-Dahab, 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]
28. G. Palani and I. A. Abbas, “Free Convection MHD Flow with Thermal Radiation from an Impulsively-Started 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]
29. 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]