In order to reduce the oscillations of the numerical solution of fractional exotic options pricing model, a class of numerical schemes are developed and well studied in this paper which are based on the 4th-order Padé approximation and 2nd-order weighted and shifted Grünwald difference scheme. Since the spatial discretization matrix is positive definite and has lower Hessenberg Toeplitz structure, we prove the convergence of the proposed scheme. Numerical experiments on fractional digital option and fractional barrier options show that the (0, 4)-Padé scheme is fast, and significantly reduces the oscillations of the solution and smooths the Delta value.
Citation: Ming-Kai Wang, Cheng Wang, Jun-Feng Yin. A class of fourth-order Padé schemes for fractional exotic options pricing model[J]. Electronic Research Archive, 2022, 30(3): 874-897. doi: 10.3934/era.2022046
[1] | Muhammad Imran Asjad, Waqas Ali Faridi, Mohammed M. Al-Shomrani, Abdullahi Yusuf . The generalization of Hermite-Hadamard type Inequality with exp-convexity involving non-singular fractional operator. AIMS Mathematics, 2022, 7(4): 7040-7055. doi: 10.3934/math.2022392 |
[2] | Muhammad Tariq, Asif Ali Shaikh, Sotiris K. Ntouyas, Jessada Tariboon . Some novel refinements of Hermite-Hadamard and Pachpatte type integral inequalities involving a generalized preinvex function pertaining to Caputo-Fabrizio fractional integral operator. AIMS Mathematics, 2023, 8(11): 25572-25610. doi: 10.3934/math.20231306 |
[3] | Eze R. Nwaeze, Muhammad Adil Khan, Ali Ahmadian, Mohammad Nazir Ahmad, Ahmad Kamil Mahmood . Fractional inequalities of the Hermite–Hadamard type for m-polynomial convex and harmonically convex functions. AIMS Mathematics, 2021, 6(2): 1889-1904. doi: 10.3934/math.2021115 |
[4] | Muhammad Tariq, Hijaz Ahmad, Abdul Ghafoor Shaikh, Soubhagya Kumar Sahoo, Khaled Mohamed Khedher, Tuan Nguyen Gia . New fractional integral inequalities for preinvex functions involving Caputo-Fabrizio operator. AIMS Mathematics, 2022, 7(3): 3440-3455. doi: 10.3934/math.2022191 |
[5] | Hong Yang, Shahid Qaisar, Arslan Munir, Muhammad Naeem . New inequalities via Caputo-Fabrizio integral operator with applications. AIMS Mathematics, 2023, 8(8): 19391-19412. doi: 10.3934/math.2023989 |
[6] | Hengxiao Qi, Muhammad Yussouf, Sajid Mehmood, Yu-Ming Chu, Ghulam Farid . Fractional integral versions of Hermite-Hadamard type inequality for generalized exponentially convexity. AIMS Mathematics, 2020, 5(6): 6030-6042. doi: 10.3934/math.2020386 |
[7] | Anjum Mustafa Khan Abbasi, Matloob Anwar . Hermite-Hadamard inequality involving Caputo-Fabrizio fractional integrals and related inequalities via s-convex functions in the second sense. AIMS Mathematics, 2022, 7(10): 18565-18575. doi: 10.3934/math.20221020 |
[8] | Manar A. Alqudah, Artion Kashuri, Pshtiwan Othman Mohammed, Muhammad Raees, Thabet Abdeljawad, Matloob Anwar, Y. S. Hamed . On modified convex interval valued functions and related inclusions via the interval valued generalized fractional integrals in extended interval space. AIMS Mathematics, 2021, 6(5): 4638-4663. doi: 10.3934/math.2021273 |
[9] | Sobia Rafeeq, Sabir Hussain, Jongsuk Ro . On fractional Bullen-type inequalities with applications. AIMS Mathematics, 2024, 9(9): 24590-24609. doi: 10.3934/math.20241198 |
[10] | Iqra Nayab, Shahid Mubeen, Rana Safdar Ali, Faisal Zahoor, Muath Awadalla, Abd Elmotaleb A. M. A. Elamin . Novel fractional inequalities measured by Prabhakar fuzzy fractional operators pertaining to fuzzy convexities and preinvexities. AIMS Mathematics, 2024, 9(7): 17696-17715. doi: 10.3934/math.2024860 |
In order to reduce the oscillations of the numerical solution of fractional exotic options pricing model, a class of numerical schemes are developed and well studied in this paper which are based on the 4th-order Padé approximation and 2nd-order weighted and shifted Grünwald difference scheme. Since the spatial discretization matrix is positive definite and has lower Hessenberg Toeplitz structure, we prove the convergence of the proposed scheme. Numerical experiments on fractional digital option and fractional barrier options show that the (0, 4)-Padé scheme is fast, and significantly reduces the oscillations of the solution and smooths the Delta value.
Fractional calculus provides a concise model for the description of the dynamic events that occur in biological tissues. Such a description is important for gaining an understanding of the underlying multiscale processes that occur when, for example, tissues are electrically stimulated or mechanically stressed. The mathematics of fractional calculus has been applied successfully in physics, chemistry, and materials science to describe dielectrics, electrodes and viscoelastic materials over extended ranges of time and frequency, see [1,2].
Fractional calculus is now a well-established tool in engineering science, with very promising applications in materials modelling. Indeed, several studies have shown that fractional operators can successfully describe complex long-memory and multiscale phenomena in materials, which can hardly be captured by standard mathematical approaches as, for instance, classical differential calculus. Furthermore, fractional calculus has recently proved to be an excellent framework for modelling non-conventional fractal and non-local media, opening valuable prospects on future engineered materials, see [3,4].
Fractional calculus become a key tool for many fields such as Biology [5,6], Economy [7], Demography [8], Geophysics [9], Medicine [10] and Bio-engineering [11]. There are many interesting controversies and generalizations for fractional calculus available to handle more and more real world problems [5,12]. This operator is significant because of their singular Kernal and hence, it is interesting to develop inequalities involving fractional order derivatives and fractional operators. Now a days, "Hermite-Hadamard type inequalities" are one of the top trend topic for the researchers of convex analysis and inequalities, which is define as:
Theorem 1.1. [13] Let ξ:J⊆R→R be a convex function and u,v∈J with u<v, then the following double inequality holds:
ξ(u+v2)≤1v−u∫vuξ(x)dx≤ξ(u)+ξ(v)2. | (1.1) |
Many fractional operators are used to generalized Hermite-Hadamard inequality, for example, Chen and Katugampola [14] presented Hermite-Hadamard and Hermite-Hadamard-Fejér type inequalities for generalized fractional integrals. In 2015, Set et al., [15] presented some new inequalities of Hermite-Hadamard-Fejer type for convex functions via fractional integrals. In 2016, Set, Sarikaya and Karakoc [16] presented Hermite-Hadamard type inequalities for convex functions via fractional integrals. Iscan in [17] presented Hermite-Hadamard type inequalities for harmonically convex functions. Gurbuz et al., in 2020, [18] studied Hermite-Hadamard inequality for fractional integrals of Caputo-Fabrizio type and presented some other related inequalities. The inequalities related to h-convexity had been studied in [19] by Varosanec. Also, there are many different extensions and versions appeared in number of papers and these extensions received many applications in different areas of mathematics, see for example [20,21,22,23,24,25,26]. For more detailed study about inequality theory and its applications, we refer to the readers [27,28,29,30,31] and the references therein.
Motivated by the work done in past years on generalizations of Hermite-Hadamard type inequalities for different convexities involving certain fractional integral operators, we developed Hermite-Hadamard type inequality for Caputo-Fabrizio fractional operator for the class of h convex functions in this paper. We also presented some other inequalities and gave applications of our results in special means. Our results are generalization and extension of many existing results in literature.
In this section, we present some known definitions and results that will help us in proving main results of this paper.
Definition 1. (Convex function) Consider an extended real valued function ξ:J→R, where J⊂Rn is any convex set, then the function ξ is convex on J if
ξ(θu+(1−θ)v)≤θξ(u)+(1−θ)ξ(v), | (2.1) |
holds for all u,v∈J and θ∈(0,1).
Definition 2. (h-convex function) [19] Let I,J be intervals in R, (0,1)⊂J and let h:J→R be a non-negative function, h≠0. A non-negative function ξ:I→R is h-convex if
ξ(θu+(1−θ)v)≤h(θ)ξ(u)+h(1−θ)ξ(v), | (2.2) |
holds for all u,v∈J and θ∈(0,1).
If the inequality (2.2) is reversed, then the function ξ is h-concave. The class of h convex functions is denoted by SX(h,J) and the class of h concave functions is denoted by SV(h,J).
In the following remark, we give the relationship between the definitions 1 and 2.
Remark 1. 1. If h(θ)=θ then (2.2) reduces to (2.1).
2. If h(θ)≥θ for all θ∈(0,1) and the function ξ is convex and non-negative than ξ∈SX(h,J).
3. If h(θ)≤θ for all θ∈(0,1)and the function ξ is convex and non-negative than ξ∈SV(h,J).
Definition 3. (Caputo-Fabrizio fractional time derivative) For any function ξ, the Caputo fractional derivative of order σ is denoted by (UFDt) and is defined as
Dσtξ(t)=1γ(1−σ)∫tuξ′(x)(t−x)σdx, | (2.3) |
with σ∈(0,1) and u∈[−∞,t),ξ∈H1(u,v),u<v, (H1(u,v) is class of first order differentiable function). By changing the kernal (t−x)−σ with the function exp(−σ(t−x)σ1−σ) and 1γ(1−σ) with B(σ)1−σ, where B(σ)>0 is a normalization function satisfying B(0)=B(1)=1, we obtained the new definition of fractional time derivative
(Dσtξ)(t)=B(σ)1−σ∫tuξ′(x)e−σ(t−x)σ1−σdx, | (2.4) |
Definition 4. [1] Let ξ∈H1(u,v),u<v,σ∈(0,1), then the left Caputo-Fabrizio fractional derivative is defined as
(CFCuDσξ)(t)=B(σ)1−σ∫tuξ′(x)e−σ(t−x)σ1−σdx, | (2.5) |
and the integral associated with this fractional derivative is
(CFuIσξ)(t)=1−σB(σ)ξ(t)+σB(σ)∫tuξ(x)dx. | (2.6) |
Here, the function B(σ)>0 is normalization that satisfy the condition B(0)=B(1)=1.
Now, the right Caputo-Fabrizio fractional derivative is defined as
(CFCDσvξ)(t)=−B(σ)1−σ∫vtξ′(x)e−σ(x−t)σ1−σdx, | (2.7) |
and the integral associated with this fractional derivative is
(CFIσvξ)(t)=1−σB(σ)ξ(t)+σB(σ)∫vtξ(x)dx. | (2.8) |
Lemma 1. Consider a differential mapping ξ:J⊆R→R defined on ˚J and u,v∈J with u<v. If ξ′∈L[u,v] then we have
ξ(u)+ξ(v)2−1v−u∫vuξ(x)dx=v−u2∫10(1−2α)ξ′(θu−(1−θ)v)dθ. |
Lemma 2. [18,Lemma 2] Consider a differential mapping ξ:J=[u,v]→R defined on ˚J, and u,v∈J with u<v. If ξ∈L1[u,v], and σ∈(0,1) then we have
v−u2∫10(1−2θ)ξ′(θu−(1−θ)v)dθ−2(1−σ)σ(v−u)ξ(k)=ξ(u)+ξ(v)2−B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)], | (2.9) |
where k∈[u,v] and B(σ)>0 is a normalization function.
In the following theorem, we present a variant of Hermite-Hadamard inequality in the setting of h-convex functions.
Theorem 3.1. Let ξ:J=[u,v]→R be an h-convex function defined on [u,v] and ξ∈L1[u,v]. If σ∈(0,1), then we have
12h(12)ξ(u+v2)≤B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)−2(1−σ)B(σ)ξ(k)]≤(ξ(u)+ξ(v))∫10h(θ)dθ, | (3.1) |
where k∈[u,v] and B(σ)>0 is as defined above.
Proof. The Hermite-Hadamard inequality for h-convex is as follows;
12h(12)ξ(u+v2)≤1v−u∫vuξ(x)dx≤(ξ(u)+ξ(v))∫10h(θ)dθ. | (3.2) |
Since ξ is a h-convex function on [u,v], we can write
22h(12)ξ(u+v2)≤2v−u∫vuξ(x)dx=2v−u(∫kuξ(x)dx+∫vkξ(x)dx). | (3.3) |
Multiplying both sides of (3.3) by σ(v−u)2B(σ) and adding 2(1−σ)B(σ)ξ(k), we get
2(1−σ)B(σ)ξ(k)+σ(v−u)B(σ)(12h(12)ξ(u+v2))≤2(1−σ)B(σ)ξ(k)+σB(σ)[∫kuξ(x)dx+∫vkξ(x)dx]=((1−σ)B(σ)ξ(k)+σB(σ)∫kuξ(x)dx)+((1−σ)B(σ)ξ(k)+σB(σ)∫vkξ(x)dx)=(CFuIσξ)(k)+(CFIσvξ)(k). | (3.4) |
After suitable rearrangement of (3.4), we arrive at the left inequality of (3.1).
Now we will prove the right side of (3.1). The Hermite-Hadamard inequality for h-convex functions is
2v−u∫vuξ(x)dx≤2[(ξ(u)+ξ(v))∫10h(θ)dθ],2v−u[∫kuξ(x)dx+∫vkξ(x)dx]≤2[(ξ(u)+ξ(v))∫10h(θ)dθ]. | (3.5) |
By using the same operator with (3.3) in (3.5), we have
(CFuIσξ)(k)+(CFIσvξ)(k)≤2(1−σ)B(σ)ξ(k)+σ(v−u)B(σ)((ξ(u)+ξ(v))∫10h(θ)dθ). | (3.6) |
After suitable rearrangement of (3.6), we get the required right side of (3.1), which complete the proof.
Following remark proofs that our result is generalization of existing result.
Remark 2. Taking h(θ)=θ in Theorem 3.1, we get [18,Theorem 2].
In the following theorem, we present another variant of Hermite-Hadamard inequality.
Theorem 3.2. Let ξ1,ξ2:J⊆R→R be a h-convex functions on J. If ξ1ξ2∈L1([u,v]), then we have the following inequality:
2B(σ)σ(v−u)[(CFuIσξ1ξ2)(k)+(CFIσvξ1ξ2)(k)−2(1−σ)B(σ)ξ1(k)ξ2(k)]≤(2∫10(h(θ))2dθ)M(u,v)+(2∫10h(θ)h(1−θ)dθ)N(u,v), | (3.7) |
where M(u,v)=ξ1(u)ξ2(u)+ξ1(v)ξ2(v), N(u,v)=ξ1(u)ξ2(v)+ξ1(v)ξ2(u) and k∈[u,v] and B(σ)>0.
Proof. By definition of h-convexity of ξ1 and ξ2, we have
ξ1(θu+(1−θ)v)≤h(θ)ξ1(u)+h(1−θ)ξ1(v),∀θ∈[0,1], | (3.8) |
and
ξ2(θu+(1−θ)v)≤h(θ)ξ2(u)+h(1−θ)ξ2(v),∀θ∈[0,1]. | (3.9) |
Multiplying both sides of (3.8) and (3.9), we have
ξ1(θu+(1−θ)v)ξ2(θu+(1−θ)v)≤(h(θ))2ξ1(u)ξ2(u)+(h(1−θ))2ξ1(v)ξ2(v)+h(θ)h(1−θ)[ξ1(u)ξ2(v)+ξ1(v)ξ2(u)]. | (3.10) |
Integrating (3.10) with respect to θ over [0,1] and changing variables, we obtain
1v−u∫vuξ1(x)ξ2(x)dx≤ξ1(u)ξ2(u)∫10(h(θ))2dθ+ξ1(v)ξ2(v)∫10(h(1−θ))2dθ+[ξ1(u)ξ2(v)+ξ1(v)ξ2(u)]∫10h(θ)h(1−θ)dθ, |
which implies
2v−u[∫tuξ1(x)ξ2(x)dx+∫vtξ1(x)ξ2(x)dx]≤2[∫10(h(θ))2dθ[ξ1(u)ξ2(u)+ξ1(v)ξ2(v)]+∫10h(θ)h(1−θ)dθ+[ξ1(u)ξ2(v)+ξ1(v)ξ2(u)]].≤2[(∫10(h(θ))2dθ)M(u,v)+(∫10h(θ)h(1−θ)dθ)N(u,v)]. | (3.11) |
Multiplying both sides of (3.11) by σ(v−u)2B(σ) and adding 2(1−σ)B(σ)ξ1(k)ξ2(k), we get
σB(σ)[∫kuξ1(x)ξ2(x)dx+∫vkξ1(x)ξ2(x)dx]+2(1−σ)B(σ)ξ1(k)ξ2(k)≤σ(v−u)B(σ)[2(∫10(h(θ))2dθ)M(u,v)+2(∫10h(θ)h(1−θ)dθ)N(u,v)]+2(1−σ)B(σ)ξ1(k)ξ2(k). |
Thus
(CFuIσξ1ξ2)(k)+(CFIσvξ1ξ2)(k)≤σ(v−u)B(σ)[2(∫10(h(θ))2dθ)M(u,v)+2(∫10h(θ)h(1−θ)dθ)N(u,v)]+2(1−σ)B(σ)ξ1(k)ξ2(k). | (3.12) |
The suitable arrangement of (3.12) completes the proof.
The following remark shows that our result is generalization of existing result.
Remark 3. Taking h(θ)=θ in Theorem 3.2 we obtain [18,Theorem 3].
Following theorem contains another variant of Hermite-hadamard inequality.
Theorem 3.3. Let ξ1 and ξ2 are h-convex functions on J. If ξ1ξ2∈L([u,v]), then we have
12[h(12)]2ξ1(u+v2)ξ2(u+v2)−B(σ)σ(v−u)[(CFuIσξ1ξ2)(k)+(CFIσvξ1ξ2)(k)]+2(1−σ)σ(v−u)ξ1(k)ξ2(k)≤M(u,v)∫10h(θ)h(1−θ)dθ+12N(u,v)∫10[(h(θ))2+(h(1−θ))2]dθ, | (3.13) |
where M(u,v)=ξ1(u)ξ1(u)+ξ1(v)ξ1(v), N(u,v)=ξ2(u)ξ2(v)+ξ2(v)ξ2(u) and k∈[u,v] and B(σ)>0.
Proof. By using h-convexity of ξ1 and ξ2 and taking θ=12, we have
ξ1(u+v2)≤[h(12)]2ξ1((1−θ)u+θv)+[h(12)]2ξ1(θu+(1−θ)v), |
and
ξ2(u+v2)≤[h(12)]2ξ2((1−θ)u+θv)+[h(12)]2ξ2(θu+(1−θ)v). |
Multiplying the above inequalities at both sides, we get
ξ1(u+v2)ξ2(u+v2)≤[h(12)]2[ξ1((1−θ)u+θv)ξ2((1−θ)u+θv)+ξ1(θu+(1−θ)v)ξ2(θu+(1−θ)v)+ξ1((1−θ)u+θv)ξ2(θu+(1−θ)v)+ξ1(θu+(1−θ)v)ξ2((1−θ)u+θv)]≤[h(12)]2[ξ1((1−θ)u+θv)ξ2((1−θ)u+θv)+ξ1(θu+(1−θ)v)ξ2(θu+(1−θ)v)+2h(θ)h(1−θ){ξ1(u)ξ2(u)+ξ1(v)ξ2(v)}{(h(θ))2+(h(1−θ))2}{ξ1(u)ξ2(v)+ξ1(v)ξ2(u)}]. | (3.14) |
Integrating (3.14) with respect to θ over [0,1] and changing variables, we obtain
ξ1(u+v2)ξ2(u+v2)≤[h(12)]2[2v−u∫vuξ1(x)ξ2(x)dx+2M(u,v)∫10h(θ)h(1−θ)dθ+N(u,v)∫10[(h(θ))2+(h(1−θ))2]dθ]. | (3.15) |
Multiplying both sides of (3.15) by σ(v−u)2B(σ) and subtracting 2(1−σ)B(σ)ξ1(k)ξ2(k), we obtain
σ(v−u)2B(σ)[h(12)]2ξ1(u+v2)ξ2(u+v2)−2(1−σ)B(σ)ξ1(k)ξ2(k)≤σB(σ)∫vuξ1(x)ξ2(x)dx+σ(v−u)2B(σ)[2M(u,v)∫10h(θ)h(1−θ)dθ+N(u,v)∫10[(h(θ))2+(h(1−θ))2]dθ]−2(1−σ)B(σ)ξ1(k)ξ2(k). |
Thus
σ(v−u)2B(σ)[h(12)]2ξ1(u+v2)ξ2(u+v2)−2(1−σ)B(σ)ξ1(k)ξ2(k)−σB(σ)[∫kuξ1(x)ξ2(x)dx+∫ukξ1(x)ξ2(x)dx+]≤σ(v−u)2B(σ)[2M(u,v)∫10h(θ)h(1−θ)dθ+N(u,v)∫10[(h(θ))2+(h(1−θ))2]dθ]−2(1−σ)B(σ)ξ1(k)ξ2(k). |
This implies that
σ(v−u)2B(σ)[h(12)]2ξ1(u+v2)ξ2(u+v2)−(CFuIσξ1ξ2)(k)+(CFIσvξ1ξ2)(k)≤σ(v−u)2B(σ)[2M(u,v)∫10h(θ)h(1−θ)dθ+N(u,v)∫10[(h(θ))2+(h(1−θ))2]dθ]−2(1−σ)B(σ)ξ1(k)ξ2(k). | (3.16) |
Multiplying (3.16) 2B(σ)σ(v−u), we obtained the required inequality (3.13).
Remark 4. If we take h(θ)=θ in Theorem 3.3 we obtain [18,Theorem 4].
In the following theorem, we present an inequality concerning Caputo-Fabrizio fractional operator in the setting of h-convexity.
Theorem 4.1. Consider a differentiable function ξ:J→R defined on ˚J such that the function |ξ′| is h-convex on [u,v], where u,v∈J with u<v. If ξ′∈L1[u,v] and σ∈(0,1), then we have
|ξ(u)+ξ(v)2+2(1−σ)σ(v−u)ξ(k)−B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)]|≤v−u2[E1|ξ′(u)|+E2|ξ′(v)|], | (4.1) |
where
E1=(∫120|1−2θ|h(θ)dθ+∫112|2θ−1|h(θ)dθ), | (4.2) |
E2=(∫120|1−2θ|h(1−θ)dθ+∫112|2θ−1|h(1−θ)dθ), | (4.3) |
where k∈[u,v] and B(σ)>0 is a normalization function.
Proof. In the light of the Lemma 2 and the fact that |ψ′| is h-convex, we get
|ξ(u)+ξ(v)2+2(1−σ)σ(v−u)ξ(k)−B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)]|≤v−u2∫10|1−2θ||ξ′(θu+(1−θ)v)|dθ≤v−u2∫10|1−2θ|[h(θ)|ξ′(u)|+h(1−θ)|ξ′(v)|]dθ=v−u2(∫120|1−2θ|[h(θ)|ξ′(u)|+h(1−θ)|ξ′(v)|]dθ+∫112|2θ−1|[h(θ)|ξ′(u)|+h(1−θ)|ξ′(v)|]dθ)=v−u2[|ξ′(u)|(∫120|1−2θ|h(θ)dθ+∫112|2θ−1|h(θ)dθ)+|ξ′(v)|(∫120|1−2θ|h(1−θ)dθ+∫112|2θ−1|h(1−θ)dθ)]=v−u2[E1|ξ′(u)|+E2|ξ′(v)|]. | (4.4) |
Which completes the proof.
Remark 5. If we take h(θ)=θ in Theorem 4.1, we obtain [18,Theorem 5].
Theorem 4.2. Consider a differentiable funciton ξ:J→R defined on ˚J such that the function |ξ′|q is h-convex on [u,v], u,v∈˚J with u<v,q>1,1p+1q=1 where u,v∈J with u<v. If ξ′∈L1[u,v], and σ∈(0,1), then we have
|ξ(u)+ξ(v)2+2(1−σ)σ(v−u)ξ(k)−B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)]|≤v−u2(1p+1)1p[|ξ′(u)|q∫10h(θ)dθ+|ξ′(v)|q∫10h(1−θ)dθ]1q, | (4.5) |
where k∈[u,v] and B(σ)>0 is a normalization function.
Proof. In the light of Lemma 2, Hölder's inequality and the fact that |ξ′|q is h-convex, we get
|ξ(u)+ξ(v)2+2(1−σ)σ(v−u)ξ(k)−B(σ)σ(v−u)[(CFuIσξ)(k)+(CFIσvξ)(k)]|≤v−u2∫10|1−2θ||ξ′(θu+(1−θ)v)|dθ |
≤v−u2[(∫10|1−2θ|pdθ)1p(∫10|ξ′(θu+(1−θ)v)|qdθ)1q]=v−u2(1p+1)1p(|ξ′(u)|q∫10h(θ)dθ+|ξ′(v)|q∫10h(1−θ)dθ)1q. |
Which completes the proof.
Remark 6. If we take h(θ)=θ in Theorem 4.2 we obtain [18,Theorem 6].
Means has significant important in applied and pure mathematics, especially the accuracy of a results can be confirmed by the application to special means for real numbers u,v such that u≠v. They are in the following order;
H≤G≤L≤I≤A | (5.1) |
For two positive numbers u>0 and v>0, the arithmetic mean is define as
A(u,v)=u+v2,u,v∈R |
The generalized logarithmic mean is defined as
L=Lrr(u,v)=vr+1−ur+1(r+1)(v−u),r∈R−[−1,0],u,v∈R,u≠v. | (5.2) |
Proposition 5.1. Let u,v∈[0,∞) with u<v, we have
|A(u2,v2)−L22(u,v)|≤(v−u)[E1|u|+E2|v|]. | (5.3) |
Proof. In Theorem 4.1, if we set ξ(x)=x2, with σ=1 and B(σ)=B(1)=1, we obtained the required result.
Remark 7. If we put h(θ)=θ, in preposition (5.1) then we will obtain [18,preposition 1].
Proposition 5.2. Let u,v∈[0,∞) with u<v, then we have
|A(eu,ev)−L(eu,ev)|≤(v−u)2[E1eu+E2ev]. | (5.4) |
Proof. In Theorem 4.1, if we set ξ(x)=ex, with σ=1 and B(σ)=B(1)=1, we obtained the required result.
Remark 8. If we put h(θ)=θ, in preposition (5.1) then we will obtain [18,preposition 2].
Proposition 5.3. Let u,v∈R+, u<v, then
|A(un,vn)−Lnn(u,v)|≤n(v−u)2[E1|un−1|+E2|vn−1|]. | (5.5) |
Proof. In Theorem 4.1, if we set ξ(x)=xn, where n is an even number with σ=1 and B(σ)=B(1)=1, we obtained the required result.
Remark 9. If we put h(θ)=θ, in preposition (5.1) then we will obtain [18,preposition 3].
The class of convex functions has special place in the theory of optimization problems due to their special properties, for example, any convex function defined on an open domain has exactly one minimum. An another important tool of mathematics is Caputo-Fabrizio integral operator that attract attentions of many mathematician and researcher working in other fields. This operator is very helpful for sake of modeling of problems of applies sciences and engineering. The fractional derivatives contribute significantly in modeling almost every phenomena in the field of Engineering and applied physics. In this paper, Hermite-Hadamard type inequalities for h-convex functions via Caputo-Fabrizio integral operator are derived. Some new and important integral inequalities involving Caputo-Fabrizio fractional integral operator are also obtained for h-convex functions. Many existing results in literature become the particular cases for these results as mentioned in remarks.
The work is partially supported by the Natural Science Foundation of Hebei Province of China (No. A2019106037), the Doctoral Scientific Research Foundation of Shijiazhuang University (No. 18BS013), and the Science and Technology R & D Program of Handan (No. 1721211052-5).
The authors are thankful to all the reviewers for their nice comments that help us to improve the quality of the paper. The present form of the paper is due to constructive comments of all five reviewers.
The authors declare that they do not have any conflict of interests.
[1] |
F. Black, M. Scholes, The pricing of options and corporate liabilities, J. Political Economy, 81 (1973), 637–654. https://doi.org/10.1086/260062 doi: 10.1086/260062
![]() |
[2] |
R. C. Merton, Theory of rational option pricing, Bell J. Econ. Manage. Sci., 4 (1973), 141–183. https://doi.org/10.2307/3003143 doi: 10.2307/3003143
![]() |
[3] |
R. C. Merton, Option pricing when underlying stock returns are discontinuous, J. Financ. Econ., 3 (1976), 125–144. https://doi.org/10.1016/0304-405X(76)90022-2 doi: 10.1016/0304-405X(76)90022-2
![]() |
[4] | S. Heston, A closed-form solution for options with stochastic volatility with applications to bond and currency options, Review of Financ. Stud., 6 (1993), 327–343. Available from: https://www.jstor.org/stable/2962057. |
[5] |
A. White, J. Hull, The pricing of options on assets with stochastic volatilities, J. Finance, 42 (1987), 281–300. https://doi.org/10.1111/j.1540-6261.1987.tb02568.x doi: 10.1111/j.1540-6261.1987.tb02568.x
![]() |
[6] |
P. Carr, L. Wu, The finite moment log stable process and option pricing, J. Finance, 58 (2003), 753–777. https://doi.org/10.1111/1540-6261.00544 doi: 10.1111/1540-6261.00544
![]() |
[7] |
I. Koponen, Analytic approach to the problem of convergence of truncated Lévy flights towards the Gaussian stochastic process, Phys. Rev. E, 52 (1995), 1197. https://doi.org/10.1103/PhysRevE.52.1197 doi: 10.1103/PhysRevE.52.1197
![]() |
[8] |
L. X. Zhang, R. F. Peng, J. F. Yin, A second order numerical scheme for fractional option pricing models, East Asian J. Appl. Math., 11 (2021), 326–348. https://doi.org/10.4208/eajam.020820.121120 doi: 10.4208/eajam.020820.121120
![]() |
[9] |
P. Carr, H. Geman, D. B. Madan, M. Yor, Stochastic volatility for Lévy processes, Math. Finance, 13 (2003), 345–382. https://doi.org/10.1111/1467-9965.00020 doi: 10.1111/1467-9965.00020
![]() |
[10] | R. Cont, E. Voltchkova, A finite difference scheme for option pricing in jump diffusion and exponential lévy models, SIAM J. Numer. Anal., 43 (2005), 1596–1626. https://doi.org/10.1137/S0036142903436186 |
[11] |
Á. Cartea, D. del Castillo-Negrete, Fractional diffusion models of option prices in markets with jumps, Physica A, 374 (2007), 749–763. https://doi.org/10.1016/j.physa.2006.08.071 doi: 10.1016/j.physa.2006.08.071
![]() |
[12] |
O. Marom, E. Momoniat, A comparison of numerical solutions of fractional diffusion models in finance, Nonlinear Anal. Real World Appl., 10 (2009), 3435–3442. https://doi.org/10.1016/j.nonrwa.2008.10.066 doi: 10.1016/j.nonrwa.2008.10.066
![]() |
[13] |
W. Chen, S. Wang, A finite difference method for pricing european and american options under a geometric Lévy process, J. Ind. Manage. Optim., 11 (2014), 241–264. https://doi.org/10.3934/jimo.2015.11.241 doi: 10.3934/jimo.2015.11.241
![]() |
[14] |
H. Zhang, F. Liu, I. Turner, S. Chen, Q. Yang, Numerical simulation of a finite moment log stable model for a European call option, Numerical Algorithms, 75 (2017), 569–585. https://doi.org/10.1007/s11075-016-0212-x doi: 10.1007/s11075-016-0212-x
![]() |
[15] |
A. Q. M. Khaliq, D. A. Voss, M. Yousuf, Pricing exotic options with L-stable Padé schemes, J. Banking Finance, 31 (2007), 3438–3461. https://doi.org/10.1016/j.jbankfin.2007.04.012 doi: 10.1016/j.jbankfin.2007.04.012
![]() |
[16] |
B. A. Wade, A. Q. M. Khaliq, M. Yousuf, J. Vigo-Aguiar, R. Deininger, On smoothing of the Crank-Nicolson scheme and higher order schemes for pricing barrier options, J. Comput. Appl. Math., 204 (2007), 144–158. https://doi.org/10.1016/j.cam.2006.04.034 doi: 10.1016/j.cam.2006.04.034
![]() |
[17] |
M. K. Kadalbajoo, A. Kumar, L. P. Tripathi, Application of radial basis function with L-stable Padé time marching scheme for pricing exotic option, Comput. Math. Appl., 66 (2013), 500–511. https://doi.org/10.1016/j.camwa.2013.06.002 doi: 10.1016/j.camwa.2013.06.002
![]() |
[18] |
A. Q. M. Khaliq, B. A. Wade, M. Yousuf, J. Vigo-Aguiar, High order smoothing schemes for inhomogeneous parabolic problems with applications in option pricing, Numer. Methods Partial Differ. Equations, 23 (2010), 1249–1276. https://doi.org/10.1002/num.20228 doi: 10.1002/num.20228
![]() |
[19] |
M. Yousuf, A fourth-order smoothing scheme for pricing barrier options under stochastic volatility, Int. J. Comput. Math., 86 (2009), 1054–1067. https://doi.org/10.1080/00207160802681653 doi: 10.1080/00207160802681653
![]() |
[20] | I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Academic Press, London, 1999. Available from: https://www.sciencedirect.com/bookseries/mathematics-in-science-and-engineering/vol/198/supp l/C. |
[21] | W. Y. Tian, H. Zhou, W. H. Deng, A class of second order difference approximation for solving space fractional diffusion equations, Math. Comput., 84 (2012), 1703–1727. Available from: https://www.jstor.org/stable/24489172. |
[22] |
Y. Liu, M. Zhang, H. Li, J. Li, High-order local discontinuous galerkin method combined with wsgd-approximation for a fractional subdiffusion equation, Comput. Math. Appl., 73 (2017), 1298–1314. https://doi.org/10.1016/j.camwa.2016.08.015 doi: 10.1016/j.camwa.2016.08.015
![]() |
[23] |
Z. B. Wang, S. W. Vong, Compact difference schemes for the modified anomalous fractional sub-diffusion equation and the fractional diffusion-wave equation, J. Comput. Phys., 277 (2014), 1–15. https://doi.org/10.1016/j.jcp.2014.08.012 doi: 10.1016/j.jcp.2014.08.012
![]() |
[24] |
Y. Liu, Y. W. Du, H. Li, F. W. Liu, Y. J. Wang, Some second-order θ schemes combined with finite element method for nonlinear fractional cable equation, Numerical Algorithms, 80 (2019), 533–555. https://doi.org/10.1007/s11075-018-0496-0 doi: 10.1007/s11075-018-0496-0
![]() |
[25] |
W. F. Wang, X. Chen, D. Ding, S. L. Lei, Circulant preconditioning technique for barrier options pricing under fractional diffusion models, Int. J. Comput. Math., 92 (2015), 2596–2614. https://doi.org/10.1080/00207160.2015.1077948 doi: 10.1080/00207160.2015.1077948
![]() |
[26] | V. Thomée, Galerkin Finite Element Methods for Parabolic Problems, Springer, Berlin, Heidelberg, 1986. https://doi.org/10.1007/3-540-33122-0 |
[27] |
A. H. Armstrong, G. D. Smith, Numerical solution of partial differential equations, Am. Math. Mon., 70 (1980), 330–332. https://doi.org/10.2307/3616228 doi: 10.2307/3616228
![]() |
[28] | J. D. Lambert, Numerical Methods for Ordinary Differential Systems: The Initial Value Problem, Mathematics of Computation, 61 (1991). Available from: https://dl.acm.org/doi/book/10.5555/129839. |
[29] |
A. Q. M. Khaliq, E. H. Twizell, D. A. Voss, On parallel algorithms for semidiscretized parabolic partial differential equations based on subdiagonal Padé approximations, Numer. Methods Partial Differ. Equations, 9 (2010), 107–116. https://doi.org/10.1002/num.1690090202 doi: 10.1002/num.1690090202
![]() |
[30] | E. Gallopoulos, Y. Saad, On the parallel solution of parabolic equations, Proc. ACM SIGARCH-89, ACM press, (1989), 17–28. https://doi.org/10.1145/318789.318793 |
[31] | A. Quarteroni, R. Sacco, F. Saleri, Numerical Mathematics, Springer, Berlin, Heidelberg, 2007. https://doi.org/10.1007/b98885 |
[32] |
B. A. Wade, A. Q. M. Khaliq, M. Siddique, M. Yousuf, Smoothing with positivity‐preserving Padé schemes for parabolic problems with nonsmooth data, Numer. Methods Partial Differ. Equations, 21 (2005), 553–573. https://doi.org/10.1002/num.20039 doi: 10.1002/num.20039
![]() |
[33] |
R. Zvan, K. R. Vetzal, P. F. Forsyth, PDE methods for pricing barrier options, J. Econ. Dyn. Control, 24 (2000), 1563–1590. https://doi.org/10.1016/S0165-1889(00)00002-6 doi: 10.1016/S0165-1889(00)00002-6
![]() |
[34] |
P. A. Forsyth, R. Zvan, K. R. Vetzal, Robust numerical methods for pde models of asian options, J. Comput. Finance, 1 (1997), 39–78. https://doi.org/10.21314/JCF.1997.006 doi: 10.21314/JCF.1997.006
![]() |
[35] | N. Zheng, J. F. Yin, High order compact schemes for variable coefficient parabolic partial differential equations with non-smooth boundary conditions, Math. Numer. Sin., 35 (2013). |
[36] | N. Zheng, J. F. Yin, C. L. Xu, Projected triangular decomposition method for pricing american option under stochastic volatility model, Commun. Appl. Math. Comput., 27 (2013). https://doi.org/10.3969/j.issn.1006-6330.2013.01.011 |
1. | Muhammad Samraiz, Maria Malik, Kanwal Saeed, Saima Naheed, Sina Etemad, Manuel De la Sen, Shahram Rezapour, A Study on the Modified Form of Riemann-Type Fractional Inequalities via Convex Functions and Related Applications, 2022, 14, 2073-8994, 2682, 10.3390/sym14122682 | |
2. | Vaijanath L. Chinchane, Asha B. Nale, Satish K. Panchal, Christophe Chesneau, Certain Weighted Fractional Inequalities via the Caputo–Fabrizio Approach, 2022, 6, 2504-3110, 495, 10.3390/fractalfract6090495 | |
3. | Mohamed Jleli, Bessem Samet, Calogero Vetro, Nonexistence Results for Higher Order Fractional Differential Inequalities with Nonlinearities Involving Caputo Fractional Derivative, 2021, 9, 2227-7390, 1866, 10.3390/math9161866 | |
4. | Jamshed Nasir, Shahid Qaisar, Ather Qayyum, Hüseyin Budak, New results on Hermite-Hadamard type inequalities via Caputo-Fabrizio fractional integral for s-convex function, 2023, 37, 0354-5180, 4943, 10.2298/FIL2315943N |