On the concept of a conformable fractional differential equation

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.


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][23][24][25][26][27][28][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 fiberreinforced 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: JOURNAL OF ENGINEERING AND THERMAL SCIENCES. JUNE 2021, VOLUME 1, ISSUE 1 Eq. (1) may be viewed as a curve in the -plane, at each point on the curve the value of its -order derivative is given in terms of and , however, refers to a particular curve that passes through a given point ( , ).
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.

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 = + Δ is defined by Odibat and Shawagfeh [10]: where Γ( ) is the gamma Function, and is the Caputo fractional derivative of order 0 < < 1 with base point , which is defined by: where is the usual first derivative, and the notation + in Eq. (2) indicates the limit as we approach 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 order fractional derivative of fractional derivative, which is not identical to the 2 order fractional derivative.
To ensure this fact we consider the function: where = + ℎ. Now the order derivative of Eq. (4) is a constant when = , and the 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 = 1, applying the well-known property of the Gamma function Γ( + 1) = !.
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 = 2, expanded about and truncated at the second-order term, is: Eq. (5) becomes: Hence, the second-order Taylor series approximation is exact for ( ), because the order of nonlinearity of the function matches the order of the Taylor series approximation. If = 3, a third-order Taylor series would provide an exact approximation.
However, if 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 > 0 is some real number. The Caputo fractional differentiation of Eq. (3) is: Thus, first-order fractional Taylor series for (4) expanded about is exact when = : the first term is = ( ). For the second term, use Eq. (7) along with the fact that the Caputo fractional derivative of the constant is zero to write: Then for any > we have (since = , the exponent − = 0, so the term ( − ) = 1 for any > ): Then of course the limit → + is: Thus, the second term of fractional Taylor series is given by: Here we have used = and Γ(1) = 1. Since ( ) is a constant, the remaining higher order Caputo derivatives are all zero. Therefore, the two-term 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 power-law function, then the two-term fractional Taylor series approximation to that function is exact. JOURNAL OF ENGINEERING AND THERMAL SCIENCES. JUNE 2021, VOLUME 1, ISSUE 1

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.

Fractional derivatives via fractional differences
Definition 1. (Fractional right derivative): Let : → be a continuous (not necessarily differentiable) function and let ℎ > 0 denotes the constant discretization span. The forward difference operator (ℎ) is defined by: the fractional difference of order (0 < ≤ 1) , from the right, of ( ) is defined by the expression: and the corresponding fractional derivative on the right is given by: 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 -order derivative of a constant is zero. . Here we are fully in the Leibnitz framework, i.e., to say both and denote the finite increments.

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 (0) = (0).
(2) Let ( ): → be a continuous function such that ( ) has fractional derivative of , where is any integer and 0 < ≤ 1, then: (3) Let us consider the compound function ( ( )). Assume that ( ) is -differentiable with respect to and ( ) is differentiable with respect to then: (4) Assume that both ( ) as well as ( ) are -differentiable with respect to and respectively, then: (5) One can extent the aforementioned fractional derivatives for negative order also i.e. of order ( < 0), can be expressed in the following integral form: For > 0 one can set: The difference between Eqs. (22) and (23) is that the second one involves the constant (0) whereas the first equation does not. We shall refer to this fractional derivative as the modified Riemann-Liouville derivative.

Fractional Taylor series with Mittag-Leffler function
Let us consider the continuous function : → , → ( ) has a fractional derivative of order , for any positive integer and , 0 < ≤ 1, then the following relation holds: where (( ) ( ) = ⋯ ( ) ( times), is the derivative of order + + ⋯ + ( times) of ( ), 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 (ℎ) satisfies the following fractional differential equation: and the solution of this equation Therefore, the series can be expressed in the following manner: In which is the differential operator with respect to , and ( ) is the Mittag-Laffler function given by: 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].

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 and (0 < < < 1) , > or < accordingly ( ) is a concave upward or concave downward function.

Now it is clear that, if ( ) is a non-decreasing function then ( ; , )
is also non-decreasing for all values for > 0 (0 < < 1).
Therefore: Thus, if ( ) is a concave upward function, for any two real numbers and (0 < < < 1), > . With the similar arguments, it can be prove that for a concave downward function ( ), < , where (0 < < < 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 ( ) where ∈ (0,1] 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 ( ) with the lowest order fractional derivative where ∈ (0,1] at that point to the curve would have been more suitable than the conventional tangent line.

Graphical illustration
Here we illustrate the nature of the tangent lines with different fractional-order derivatives ( ) where ∈ (0,1] , 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 ( ) be an approximation in the interval ( , ) . Expanding ( ) about the point = , we have: where ( ) is the order derivative of ( ) at = and ≤ ≤ .
Assuming the variation of is negligible in ≤ ≤ the truncation error can be put as ≈ ℎ . Hence, to determine the function ( ), which will be the solution of the initial value problem Eq. (1), may apply the following algorithm: where is a constant to be chosen in such a way that = max or min{ , ∈ (0,1]} accordingly the curve is concave upward or downward i.e., > or < 0 at ( , ).

Convergence criteria
A method is convergent if, for every ODE with a Lipschitz function and every fixed , with = . ℎ, it holds that: Global error: = − ( ) and = − ( ).
Local error: Therefore, local error definition implies the residual: Now comparing with the Taylor series expansion Eq. (28), the local error for this proposed Thus, the global error recursion is given by: Take norms and use Lipschitz condition: Proof. Standard result of sequence of real nos. Theorem 2. Prove that the global error , → 0 as ℎ approaches to zero. Proof. Assuming that the function is sufficiently differentiable, given ℎ and a fixed = . ℎ, let us consider , = , − ( ).
Applying the lemma to global error recursion, we get: with: Since: We have, for ≤ : Therefore: which proves the convergence as: Hence, the iteration converges invariably in keeping with expectation.

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

Conclusions
With the help of the geometrical interpretation of fractional order derivatives for different values of the parameter ∈ (0,1] 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.