On Caputo modification of Hadamard-type fractional derivative and fractional Taylor series

*Correspondence: monibashams@gmail.com 1Department of Mathematics, School of Natural Sciences, National University of Sciences and Technology, Islamabad, Pakistan Abstract In this paper a general framework is presented on some operational properties of Caputo modification of Hadamard-type fractional differential operator along with a new algorithm proposed for approximation of Hadamard-type fractional integral using Haar wavelet method. Moreover, a generalized Taylor expansion based on Caputo–Hadamard-type fractional differential operator is also established, and an example is presented for illustration.


Introduction
The history of fractional operators dates back to the seventeenth century [1]. However, interest in this theory gained momentum after the publication of [2]. Later, Riemann constructed an integral-based fractional operator [3] which is the well-known Riemann-Liouville (R-L) integral. It proved to be a valuable research in fractional calculus and is being used extensively in this field ever since. The readers are referred to a comprehensive background history of fractional calculus and details on the R-L operators given in [4].
In the last two decades scientists and researchers presented findings dealing with the fractional calculus theory and found that integer order derivatives do not work adequately in many cases. Fractional calculus has played a very important role in fields such as turbulent fluid flows, blood flow through biological tissues, control theory, and signal and image processing, fractional Brownian motion, power law models, history-dependent process, porous media, fractional filters, etc. The detailed literature including theory and applications of fractional calculus can be found in [5][6][7][8][9][10][11].
There are different definitions of fractional derivatives available in the literature. However, the most commonly used is the Hadamard fractional derivative given by Hadamard [12]. Butzer et al. [13][14][15][16] studied various properties of Hadamard-type derivatives which are more generalized than the Hadamard fractional derivatives. In this context, the readers are also referred to [1] for a detailed study on generalized fractional derivatives and references therein. Caputo introduced another type of fractional derivative [17] which has an advantage over R-L derivative in a differential equation since it does not require to define the fractional order initial conditions (see, for example, [18][19][20][21]). Recently the authors in [22,23] utilized Caputo-type modification on Hadamard factional derivatives. Moreover Trujillo et al. in [24] derived Taylor's formula with Riemann-Liouville derivatives, and Odibat and Shawagfeh [25] derived the same based on Caputo fractional derivative. Gulsu et al. [26] extended the work of previous authors and proposed a numerical scheme to approximate solutions of relaxation-oscillation equation by using the fractional Taylor series. Finally, Fernandez and Baleanu [27] developed mean value theorem and Taylor's theorem for certain fractional differential operators.
The present research presents an extension of Haar wavelet method to approximate Hadamard-type fractional integrals for a wide range of functions, and a new operational matrix of fractional integration of Haar wavelets is derived for Hadamard-type fractional integrals. A Caputo modification of Hadamard-type fractional derivative, motivated through [22,23], is also carried out in this work, and composition properties of Caputo modification of Hadamard-type fractional derivative are discussed in detail. Lastly, motivated by [28] and the aforementioned studies, we also establish a generalized Taylor expansion based on Caputo-Hadamard-type fractional differential operators and present an example for illustration.

Preliminaries
In this section some definitions and basic concepts needed in the sequel are presented precisely. We refer the interested reader to Samko et al. [3] and Kilbas et al. [7] for more details.
The left-sided Hadamard fractional integral is defined as follows.

Definition 2.1
For f ∈ L p [a, b], 1 ≤ p < ∞, the Hadamard fractional integral of order α > 0 is defined as Hadamard fractional derivative of order α > 0 is given by such that m = α and δ = x d dx .
We introduce the notation δ = x d dx and D c = x -c δ, where c ∈ R. Define a space of functions The Hadamard-type fractional integral was introduced in [15]. This operator is formally defined as follows.

Definition 2.4
For α > 0, c ∈ R, the Hadamard-type fractional integral is defined as The Hadamard-type fractional differential operator is defined as follows.
Then the Caputo modification of Hadamard-type derivative is given by

Lemma 2.9
For α ≥ 0, m = α , and μ > -1, consider the function the Caputo-Hadamard-type derivative of f is given by

Approximation of Hadamard-type fractional integral
where c ∈ R, It is to be noted that J is the maximal resolution level. It is observed that the Haar wavelets and modified Haar wavelets are related as is the ith Haar wavelet [32].
To approximate Hadamard-type integral by Haar wavelets, the Hadamard-type integrals of Haar wavelets is required. Using the definition of H i (x), these integrals can be analytically evaluated using Lemma 2.6.
For x ∈ [a, b], For x ∈ [η 1 (i), η 2 (i)), Finally, for x ≥ n 3 (i), In summary, we have A function f ∈ L 1 [a, b] is represented by Haar-type wavelet as In applications, a function f is approximated by the kth partial sum of the Haar series as We now introduce an operational matrix of Hadamard-type integration using the modified Haar wavelets. Let M(x) = [H 0 (x), H 1 (x), H 2 (x), . . . , H k-1 (x)], then at collocation points are

The Hadamard-type fractional integration of Haar wavelets yields
where P k×k is the operational matrix of Hadamard-type integration for Haar type wavelets. The approximation of the fractional Hadamard-type integral is given as follows: The numerical and exact values of the integral I α a,c f (x) are shown in Table 1. Plots for different values of c along with graphs of absolute error are shown in Figs. 1 and 2.
Example 3.2 Let g(x) = x α-c sin(3αx). The Hadamard-type fractional integral of g is computed for 1 ≤ α ≤ 2, 1 ≤ x ≤ 8. The graph for numerical values is presented in Fig. 3. The resolution level is taken as J = 6.  Proof The proof follows from Eq. (4.5) and Theorem 2.1 in [22].
For composition properties of Hadamard-type integral and Riemann-Liouville Hadamard-type derivative, the reader is referred to [13][14][15][16]. In the following section, we present an important property for Caputo-Hadamard-type operators, which plays an important role in the qualitative analysis of fractional differential equations.

(4.2)
Proof For m = 1, by definitions of I 1 a,c and D 1 c , we have Assume that identity (4. Finally, using (4.2), we get It is a common practice to transform an initial value problem or boundary value problem to an equivalent Volterra integral equation to perform qualitative analysis. In this context, the following property for Caputo-Hadamard-type operator, which is very useful for the qualitative analysis of Caputo-Hadamard differential equations, is established.
Proof The proof follows from the definition of C D α c , the semigroup property of I α a,c , and Lemma 4.
For convenience, the Leibniz rule for Hadamard-type calculus is stated. The proof of the following lemma follows from the standard Leibniz rule, where the notation is used to represent the partial derivative.

Lemma 4.4 Assume that h(x, t) and
Proof By Definition 2.8, Lemma 4.4, and the semigroup property of I α c , we have Therefore which establishes the required identity.
Remark 4.7 From Definition 2.4, we notice an obvious relation between Hadamard-type integral I α a,c and Hadamard integral I α a . Also, from Definition 2.7 and Eq. A similar relation holds between a Caputo modification of Hadamard-type derivative C D α a,c and a Caputo modification of Hadamard derivative C D α a , that is,

Taylor's formula
In this section, Taylor's formula involving Caputo-Hadamard-type derivative is discussed.