
In the present paper, we investigated some conditions to be in the class of Carathéodory functions by using the concept of the first-order differential subordinations. Moreover, various interesting special cases were considered in the geometric function theory as applications of main results presented here.
Citation: Inhwa Kim, Young Jae Sim, Nak Eun Cho. First-order differential subordinations associated with Carathéodory functions[J]. AIMS Mathematics, 2024, 9(3): 5466-5479. doi: 10.3934/math.2024264
[1] | Humaira Kalsoom, Muhammad Amer Latif, Muhammad Idrees, Muhammad Arif, Zabidin Salleh . Quantum Hermite-Hadamard type inequalities for generalized strongly preinvex functions. AIMS Mathematics, 2021, 6(12): 13291-13310. doi: 10.3934/math.2021769 |
[2] | Asghar Ahmadkhanlu, Hojjat Afshari, Jehad Alzabut . A new fixed point approach for solutions of a p-Laplacian fractional q-difference boundary value problem with an integral boundary condition. AIMS Mathematics, 2024, 9(9): 23770-23785. doi: 10.3934/math.20241155 |
[3] | Saad Ihsan Butt, Muhammad Nasim Aftab, Hossam A. Nabwey, Sina Etemad . Some Hermite-Hadamard and midpoint type inequalities in symmetric quantum calculus. AIMS Mathematics, 2024, 9(3): 5523-5549. doi: 10.3934/math.2024268 |
[4] | Changlong Yu, Jing Li, Jufang Wang . Existence and uniqueness criteria for nonlinear quantum difference equations with p-Laplacian. AIMS Mathematics, 2022, 7(6): 10439-10453. doi: 10.3934/math.2022582 |
[5] | Muhammad Uzair Awan, Muhammad Aslam Noor, Tingsong Du, Khalida Inayat Noor . On M-convex functions. AIMS Mathematics, 2020, 5(3): 2376-2387. doi: 10.3934/math.2020157 |
[6] | Reny George, Sina Etemad, Fahad Sameer Alshammari . Stability analysis on the post-quantum structure of a boundary value problem: application on the new fractional (p,q)-thermostat system. AIMS Mathematics, 2024, 9(1): 818-846. doi: 10.3934/math.2024042 |
[7] | Suphawat Asawasamrit, Muhammad Aamir Ali, Hüseyin Budak, Sotiris K. Ntouyas, Jessada Tariboon . Quantum Hermite-Hadamard and quantum Ostrowski type inequalities for s-convex functions in the second sense with applications. AIMS Mathematics, 2021, 6(12): 13327-13346. doi: 10.3934/math.2021771 |
[8] | Donny Passary, Sotiris K. Ntouyas, Jessada Tariboon . Hilfer fractional quantum system with Riemann-Liouville fractional derivatives and integrals in boundary conditions. AIMS Mathematics, 2024, 9(1): 218-239. doi: 10.3934/math.2024013 |
[9] | Nusrat Raza, Mohammed Fadel, Kottakkaran Sooppy Nisar, M. Zakarya . On 2-variable q-Hermite polynomials. AIMS Mathematics, 2021, 6(8): 8705-8727. doi: 10.3934/math.2021506 |
[10] | Chanon Promsakon, Muhammad Aamir Ali, Hüseyin Budak, Mujahid Abbas, Faheem Muhammad, Thanin Sitthiwirattham . On generalizations of quantum Simpson's and quantum Newton's inequalities with some parameters. AIMS Mathematics, 2021, 6(12): 13954-13975. doi: 10.3934/math.2021807 |
In the present paper, we investigated some conditions to be in the class of Carathéodory functions by using the concept of the first-order differential subordinations. Moreover, various interesting special cases were considered in the geometric function theory as applications of main results presented here.
Fractional calculus has been a popular topic among scholars for almost three centuries. It is a branch of mathematical analysis that deviates from classical calculus and extends the concepts of derivatives and integrals to fractional-order. Differential equations under fractional-order are known as fractional differential equations (FDEs), and they are used in various domains of real-world phenomena problems. It can be distinguished based on various criteria, such as the linearity of the equation (linear or nonlinear). Its distinguishing traits, including flexibility, memory, and hereditary properties, make it a powerful tool for modeling, evaluating, and regulating complex systems in many disciplines, including applied science and engineering (see the books of Podlubny [1], Hilfer [2], Kilbas, et al. [3], and Diethelm [4]), and a comprehensive work on boundary value problems with fractional-order [5].
It is well known that fractional derivative operators (FDOs) are often defined in the sense of fractional integral operators (FIOs) under the gamma function. In the field of fractional calculus, various definitions of fractional derivatives have emerged, such as Riemann-Liouville (RL), Caputo, Hadamard, Erdélyi-Kober, Katugampola, Hilfer, proportional, and so on, and each has different uses. The RL-FIOs and RL-FDOs concerning a function ψ are studied as in [3]. These were used to develop the ψ-Hilfer-FDO [6]. One FDO gaining attention is the proportional FDO [7]. This operator has exponential functions in the kernels, which are more advantageous than other fractional operators. For work regarding qualitative theory on proportional FDEs, see [8,9,10,11,12,13,14] and references cited therein. In parallel, the k-gamma function defines and develops fractional calculus concepts, especially the k-RL-FIOs and the k-RL-FDOs [15,16]. After that, in 2018, the (k,ψ)-RL fractional operators were proposed by [17]. Later, in 2021, Kucche and Mali proposed FDO under the (k,ψ)-Hilfer type [18], which attracted many scholars' curiosity. There are some other intriguing pieces, including those of Aljaaidi et al. [19], which proved some properties of the (k,ψ)-proportional fractional operator ((k,ψ)-PFO) and presented a new technique to investigate the k-Pólya-Szegö integral inequalities in 2021. In the same year, Boucenna et al. [20] studied the existence and uniqueness of results and a numerical technique for the Caputo proportional fractional Cauchy-type problem. While a numerical technique based on the decomposition formula was used to solve the fractional Caputo-Katugampola derivative operator [21], the analytical solutions of fractional differential equations are difficult and complex, so numerical techniques are popular for solving these equations [22,23,24]. Sudsutad et al. [25] developed the (k,ψ)-Hilfer-PFO, which cooperates the (k,ψ)-proportional fractional derivative of RL's and Caputo's senses ((k,ψ)-RL-PFDO / (k,ψ)-Caputo-PFDO). They constructed the Laplace transform concerning a function ψ of the proposed operator and applied it to solving the initial value Cauchy-type problems. In the extension, a Cauchy-type problem for FDEs is commonly defined as an initial value problem in which the goal is to find a function corresponding to the equation and the provided beginning conditions. Furthermore, numerous researchers have popularly expressed their robust findings by using a realistic application of fractional calculus in the context of diverse operators for various real-world issues such as medicine, engineering, electrical, science, and finance. For example, in 2019, Qureshi et al. [26] proposed the fractional modeling of the blood ethanal concentration using three fractional operators such as Caputo, Atangana-Baleanu, and Caputo-Fabrizio with the real data. In 2021, Norouzi and N'Guérékata [27] studied FDEs in the sense of ψ-Hilfer-FDO and used the financial crisis as an application. In 2022, Awadalla et al. [28] used a FDE in the sense of ψ-Caputo-FDO to study the model of drug concentration. In 2024, Wanassi and Torres [29] utilized the blood alcohol model as an application of the fractional analysis based on the ψ-Caputo-FDO. We refer the reader for more works in [30,31,32,33,34,35,36]. On the other hand, Ulam stability is a popularly efficient tool for ensuring that the approximate solutions generated by numerical methods remain close to the exact solutions. Numerous researchers have provided this to analyze the mathematical stability of solutions in various fields, including fractional calculus and control theory. It was initially created by Ulam in 1940 [37]. The following year, Hyer designed the Ulam-Hyers stability [38]. Presently, stability in the context of the Ulam's type has developed into various stabilities, such as Ulam-Hyers-Rassias stability [39] and Ulam-Hyer-Mittag-Leffler (UH-ML) stability, which are extensions of Ulam stability that incorporates the Mittag-Leffler function [40], and so on.
Motivated by the works [9,20,25], our major goals are to investigate the qualitative results and numerical approximations for the following Cauchy-type problem under (k,ψ)-Caputo-PFDO
{Ca,kDα,ρ;ψu(τ)=f(τ,u(τ)),0<α≤1,0<ρ≤1,u(a)=ua,ua∈R,τ∈[a,b],0≤a<b<+∞, | (1.1) |
where Ca,kDα,ρ;ψ is the (k,ψ)-Caputo-PFDO of α and f∈C([a,b]×R,R). To the best of our knowledge, this problem has yet to be considered. Specifically, we provide an extended Gronwall inequality under the (k,ψ)-PFOs to establish bounds on solutions, which is crucial in stability analysis. We investigate the existence and uniqueness of results utilizing the standard fixed point theory of Banach's and Leray-Schauder's types. Moreover, various UH–MT stability results are studied utilizing nonlinear functional analysis techniques. In addition, we develop a novel numerical approach based on a decomposition formula for solving the Cauchy-type problems. In the end, a blood alcohol level problem is presented as an application to the proposed system. The following structure governs the remaining sections of this work: In Sect. 2, we introduce the fundamental principles and some properties of the (k,ψ)-PFDO and the (k,ψ)-PFIO. The required lemmas utilized throughout this paper are also presented. In addition, we look at the extended Gronwall inequality under the aforementioned operators. The qualitative results for a Cauchy-type problem are accomplished in the remaining sections. In Section 3, we investigate the existence of the solution using a fixed point theory of Leray-Schauder's type, while the uniqueness of the solution is proved using Banach's contraction mapping principle. Then, a variety of UH-ML stability results are established to ensure the results. In Section 4, we demonstrate the numerical approach based on a decomposition formula for solving the Cauchy-type problems under the (k,ψ)-Caputo-PFDO. Some illustrative examples, particularly the blood alcohol level problem, are provided to help the accuracy of the theoretical results found in Section 5. In the final section, we outline the course of our work.
Suppose that C([a,b],R) is the Banach space of the continuous function u on [a,b] supplemented with the supremum norm ‖u‖=supτ∈[a,b]{|u(τ)|}. The space of the n-times absolutely continuous function u on [a,b] is defined by ACn([a,b],R)={u:[a,b]→R;u(n−1)∈AC([a,b],R)} and Lp([a,b],R) is the Banach space of all Lebesgue measurable g:[a,b]→R supplemented with ‖g‖Lp<+∞. Assume that ψ:[a,b]→R is a strictly increasing continuous function under ψ′(τ) not equal zero. For ease of calculation through this work, we provide a notation as follows:
ρkΨαk−1ψ(τ,s)=eρ−1kρ(ψ(τ)−ψ(s))(ψ(τ)−ψ(s))αk−1. | (2.1) |
Now, we recall some definitions and lemmas of the (k,ψ)-PFDO and the (k,ψ)-PFIO that will be applied in this work.
Definition 2.1. ([25]). Let α, k∈R+, 0<ρ≤1, and f∈L1([a,b],R). Then, the (k,ψ)-RL-PFIO of α of f is defined by
a,kIα,ρ;ψf(τ)=1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s)ds, |
where Γk(z)=∫∞0sz−1e−skkds, z∈C, Re(z)>0 and
Γk(z+k)=zΓk(z),Γk(k)=1,Γk(z)=kzk−1Γ(zk),Γ(z)=limk→1Γk(z). | (2.2) |
Definition 2.2. ([25]). Let α, k∈R+, 0<ρ≤1, f∈C([a,b],R), ψ(τ)∈Cn([a,b],R) with ψ′(τ)≠0, and n=1,2,…, such that n=⌊α/k⌋+1. Then, the (k,ψ)-RL-PFDO for α of f is defined by
RLa,kDα,ρ;ψf(τ)=kDn,ρ;ψ(a,kInk−α,ρ;ψf(τ))=kDn,ρ;ψρnk−αkkΓk(nk−α)∫τaρkΨnk−αk−1ψ(τ,s)ψ′(τ)f(s)ds, |
where kDn,ρ;ψ=kDρ;ψkDρ;ψ⋯kDρ;ψ⏟ntimes and kD1,ρ;ψf(τ)=kDρ;ψf(τ)=(1−ρ)f(τ)+kρf′(τ)ψ′(τ).
Definition 2.3. ([25]). Let α, k∈R+, 0<ρ≤1, f∈Cn(J,R), ψ(τ)∈Cn([a,b],R) with ψ′(τ)≠0, and n∈1,2,…, such that n=⌊α/k⌋+1. Then, the (k,ψ)-Caputo-PFDO for α of f is defined by
Ca,kDα,ρ;ψf(τ)=a,kInk−α,ρ;ψ(kDn,ρ;ψf(τ))=1ρnk−αkkΓk(nk−α)∫τaρkΨnk−αk−1ψ(τ,s)ψ′(s)(kDn,ρ;ψf(s))ds. |
Next, we give some important properties that are applied in this work.
Lemma 2.4. ([25]). Let α, δ∈R+∪{0}, k, η∈R+, 0<ρ≤1, ω∈R, ω/k>−1, and n=⌊ω/k⌋+1. Then,
(i) a,kIα,ρ;ψ[ρkΨωk−1ψ(τ,a)]=Γk(ω)ραkΓk(ω+α)ρkΨω+αk−1ψ(τ,a).
(ii) Ca,kDα,ρ;ψ[ρkΨωk−1ψ(τ,a)]=ραkΓk(ω)Γk(ω−α)ρkΨω−αk−1ψ(τ,a). Particularly, for m=0,1,…,n−1, we obtain Ca,kDα,ρ;ψ[ρkΨmψ(τ,a)]=0.
(iii) a,kIα,ρ;ψ[a,kIδ,ρ;ψf(τ)]=a,kIδ+α,ρ;ψf(τ)=a,kIδ,ρ;ψ[a,kIα,ρ;ψf(τ)].
(iv) Ca,kDω,ρ;ψ(a,kIη,ρ;ψf(τ))=a,kIη−ω,ρ;ψf(τ), where ω∈R+, ω<k, and η>nk.
(v) a,kIα,ρ;ψ(Ca,kDα,ρ;ψf(τ))=f(τ)−n∑i=1ρkΨn−ik−1ψ(τ,a)(ρk)n−i(n−i)!kDn−i,ρ;ψf(a), α∈(n−1,n].
Now, we prove an extended Gronwall inequality under the (k,ψ)-PFOs, and its properties are analyzed in the below Theorem.
Theorem 2.5. (An Extended (k,ψ)-Proportional Fractional Gronwall Inequality). Suppose that α>0, k>0, ρ∈(0,1], and ψ∈C1([a,b],R) is an increasing function so that ψ′(τ) not equal zero for every τ∈[a,b]. Assume that the following hypotheses hold:
(H1) The two non-negative functions u(τ) and v(τ) are locally integrable on [a,b];
(H2) The function w(τ) is a non-decreasing, non-negative, and continuous function defined on [a,b] such that w(τ)≤w∗∈R.
If
u(τ)≤v(τ)+Γk(α)kw(τ)a,kIα,ρ;ψu(τ), | (2.3) |
then, for any τ∈[a,b], we obtain
u(τ)≤v(τ)+∫τa[∞∑n=1[Γk(α)w(τ)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)v(s)]ds. | (2.4) |
Proof. First, we define an operator
Au(τ)=w(τ)ραkk2∫τaρkΨαk−1ψ(τ,s)ψ′(s)u(s)ds. | (2.5) |
From (2.5), we have u(τ)≤v(τ)+Au(τ). By applying the monotonicity of the operator A, we obtain inequalities as follows:
u(τ)≤v(τ)+Au(τ)≤v(τ)+Av(τ)+A2u(τ)≤2∑i=0Aiv(τ)+A3u(τ). | (2.6) |
By applying iterative technique, n=1,2,…, which yields that u(τ)≤∑n−1i=0Aiv(τ)+Anu(τ), τ∈[a,b], where A0v(τ)=v(τ). Next, we claim that
Anu(τ)≤∫τa[Γk(α)w(τ)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)u(s)ds, | (2.7) |
and Anu(τ)→0 as n→∞ for any τ∈[a,b]. If n=1, we get that the inequality (2.7) holds. Next, assume that the inequality (2.7) is held under n=m, that is
Amu(τ)≤∫τa[Γk(α)w(τ)]mρmαkkm+1Γk(mα)ρkΨmαk−1ψ(τ,s)ψ′(s)u(s)ds. |
If n=m+1, using the induction procedure, we have
Am+1u(τ)=A(Amu(τ))≤A(∫τa[Γk(α)w(τ)]mρmαkkm+1Γk(mα)ρkΨmαk−1ψ(τ,s)ψ′(s)u(s)ds)≤w(τ)ραkk2∫τaρkΨαk−1ψ(τ,s)ψ′(s)(∫sa[Γk(α)w(s)]mρmαkkm+1Γk(mα)ρkΨmαk−1ψ(s,r)ψ′(r)u(r)dr)ds. |
Since the function w(τ) is a non-decreasing, w(s)≤w(τ), for every s≤τ and taking z=(ψ(s)−ψ(r))/(ψ(τ)−ψ(r)), then
Am+1u(τ)≤(wm+1(τ)(Γk(α))mρ(m+1)αkkm+3Γk(mα))∫τa∫saeρ−1kρ(ψ(τ)−ψ(s))eρ−1kρ(ψ(s)−ψ(r))×(ψ(τ)−ψ(s))αk−1(ψ(s)−ψ(r))mαk−1ψ′(s)ψ′(r)u(r)drds=(wm+1(τ)(Γk(α))mρ(m+1)αkkm+3Γk(mα))∫τaeρ−1kρ(ψ(τ)−ψ(r))ψ′(r)u(r)×(∫τr(ψ(τ)−ψ(s))αk−1(ψ(s)−ψ(r))mαk−1ψ′(s)ds)dr=(wm+1(τ)(Γk(α))mρ(m+1)αkkm+2Γk(mα))∫τaρkΨ(m+1)αk−1ψ(τ,s)ψ′(r)u(r)(1k∫10(1−z)αk−1zmαk−1dz)dr=(wm+1(τ)(Γk(α))m+1ρ(m+1)αkkm+2Γk((m+1)α))∫τaρkΨ(m+1)αk−1ψ(τ,s)ψ′(r)u(r)dr. |
Since w∈C([a,b],R), then there is w∗∈R such that w(τ)≤w∗, for all τ∈[a,b], one has
Anu(τ)≤(w∗Γk(α)ραkk)n1kΓk(nα)∫τaρkΨnαk−1ψ(τ,s)ψ′(r)u(r)dr. | (2.8) |
Since 0≤eρ−1kρ(ψ(τ)−ψ(r))≤1, and u is non-negative and locally integrable on [a,b]. Hence, u is bounded on [a,b], and there is Mu∈R such that |u(τ)|≤Mu. The inequality (2.8) can be obtained
Anu(τ)≤(w∗Γk(α)(ψ(τ)−ψ(a))αkραkk)nMunαΓk(nα). | (2.9) |
Applying the Stirling's formula, that is n!∼√2πn(n/e)n, into (2.9) with (2.2), we obtain that
Anu(τ)≤Mu√k2πα(θnnnαk+12),θ:=w∗Γ(αk)(ψ(τ)−ψ(a))αkραkk2(keα)αk. | (2.10) |
where θ, ρ, k, α∈R+. This yields that, if n→∞, we obtain that Anu(τ)→0. Then,
u(τ)≤v(τ)+∫τa[∞∑n=1[Γk(α)w(τ)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)v(s)]ds. |
The inequality (2.4) is achieved.
Corollary 2.6. Let α, k∈R+, 0<ρ≤1, and ψ∈C1([a,b],R) be an increasing function such that ψ′(τ)≠0, τ∈[a,b]. Let u(τ) and v(τ) be two non-negative locally integrable functions on [a,b], and w(τ)≡M≥0. If
u(τ)≤v(τ)+MΓk(α)ka,kIα,ρ;ψu(τ), | (2.11) |
then,
u(τ)≤v(τ)+∫τa[∞∑n=1[MΓk(α)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)v(s)]ds. | (2.12) |
Corollary 2.7. Assume all conditions in Theorem 2.5 are held, and the function v(τ) is non-decreasing on τ∈[a,b]. Hence, we obtain the following inequality:
u(τ)≤v(τ)Ek,α,k((ραkk)−1Γk(α)w(τ)(ψ(τ)−ψ(s))αk), | (2.13) |
where
Ek,α,β(z)=∞∑n=0znΓk(nα+β),z∈R,α,β∈C,Re(α)>0,k>0. |
Proof. From the inequality (2.12) with the help of the non-decreasing property of v(τ), that is v(s)≤v(τ) for any τ∈[a,b], we obtain that
u(τ)≤v(τ)+∫τa[∞∑n=1[Γk(α)w(τ)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)v(s)]ds≤v(τ)(1+∫τa[∞∑n=1[Γk(α)z(τ)]nρnαkkn+1Γk(nα)ρkΨnαk−1ψ(τ,s)ψ′(s)]ds). |
Since 0≤eρ−1kρ(ψ(τ)−ψ(s))≤1, which yields that
u(τ)≤v(τ)(1+∞∑n=1[Γk(α)w(τ)]nρnαkkn+1Γk(nα)∫τa(ψ(τ)−ψ(s))nαk−1ψ′(s)ds)=v(τ)(1+∞∑n=1[Γk(α)z(τ)]nρnαkknΓk(nα+k)(ψ(τ)−ψ(a))nαk)=v(τ)∞∑n=0[Γk(α)w(τ)]nρnαkknΓk(nα+k)(ψ(τ)−ψ(a))nαk. |
Applying Definition 2.7, the inequality (2.13) is obtained.
Here, we analyze the qualitative results for the proposed problem, including existence and uniqueness results and various results for Ulam's stability. First, we prove an integral equation is equivalent to the proposed problem (1.1).
Lemma 3.1. Suppose that f∈C([a,b]×R,R), α, ρ∈(0,1] and k∈R+. Then, the proposed problem (1.1) can be stated equivalently as
u(τ)=uaeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,u(s))ds. | (3.1) |
Proof. Let u be a solution to the problem (1.1). Taking the operator a,kIα,ρ;ψ into the proposed problem (1.1) and applying (v) in Lemma 2.4, which implies that
u(τ)=u(a)eρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,u(s))ds. | (3.2) |
By using u(a)=ua into (3.2) and inserting the obtained value into itself, we get the achieved (3.1).
On the other hand, by applying (ii) and (iv) in Lemma 2.4 into (3.1), it follows that
Ca,kDα,ρ;ψu(τ)=Ca,kDα,ρ;ψ[uaeρ−1kρ(ψ(τ)−ψ(a))]+Ca,kDα,ρ;ψ(a,kIα,ρ;ψf(τ,u(τ)))=uaCa,kDα,ρ;ψ[ρkΨ0ψ(τ,s)]+a,kIα−α,ρ;ψf(τ,u(τ))=f(τ,u(τ)). |
The proof is done.
From Lemma 3.1, we provide the operator Q:C([a,b],R)→C([a,b],R) by
(Qu)(τ)=uaeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,u(s))ds. | (3.3) |
The first result is based on the Banach's contraction mapping [41].
Theorem 3.2. Suppose that α, ρ∈(0,1], k∈R+, and f∈C([a,b]×R,R). Suppose that
(A1) There is a positive real constant L such that
|f(τ,u)−f(τ,v)|≤L|u(τ)−v(τ)|,τ∈[a,b],u,v∈R. |
Then, the proposed problem (1.1) has a unique solution on [a,b], provided that
L(ψ(b)−ψ(a))αkραkΓk(α+k)<1. | (3.4) |
Proof. Assume that Br1:={u∈C([a,b],R):‖u‖≤r1} is a closed bounded and convex subset of C([a,b],R), where the radius r1 corresponds to the following condition
r1≥(ua+F(ψ(b)−ψ(a))αkραkΓk(α+k))(1−L(ψ(b)−ψ(a))αkραkΓk(α+k))−1,F=supτ∈[a,b]|f(τ,0)|<+∞. |
We prove that (ⅰ) QBr1 is bounded, i.e. QBr1⊆Br1 and (ⅱ) Q is a contraction.
Step (i). We show that QBr1⊆Br1.
For any u∈Br1 with the property of 0≤eρ−1kρ(ψ(τ)−ψ(a))≤1 and the assumption (A1), we have
|(Qu)(τ)|=supτ∈J{uaeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,u(s))ds}≤ua+1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)|f(s,u(s))|ds≤ua+1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)[|f(s,u(s))−f(s,0)|+|f(s,0)|]ds≤ua+1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)[L‖u‖+F]ds≤ua+(ψ(b)−ψ(a))αkραkΓk(α+k)[Lr1+F]≤r1. |
This implies that QBr1⊆Br1.
Step (ii). We show that Q is a contraction.
Let u, v∈C([a,b],R). Then, for every τ∈[a,b], we get
|(Qu)(τ)−(Qv)(τ)|≤supτ∈J{1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)[f(s,u(s))−f(s,v(s))]ds}≤1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)|f(s,u(s))−f(s,v(s))|ds≤1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)L|u(s)−v(s)|ds≤L‖u−v‖ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)ds=(ψ(b)−ψ(a))αkραkΓk(α+k)L‖u−v‖. |
Then, ‖Qu−Qv‖≤(ψ(b)−ψ(a))αkραkΓk(α+k)L‖u−v‖. Since the condition (3.4) holds, which implies that Q has a contraction property. Hence, by [41], Q has a unique fixed point. Then, there is a unique solution of the proposed problem (1.1) on [a,b].
Next, the existence result is achieved by utilizing the fixed point theory of Leray-Schauder's [42].
Theorem 3.3. Suppose that
(A2) There are a continuous non-decreasing function h:[a,∞)→(0,∞) and a function ϕ∈C([a,b],R+) such that the following condition:
|f(τ,u(τ))|≤ϕ(τ)h(‖u‖),∀(τ,u)∈[a,b]×R. |
(A3) There is a number N∈R+ such that the following inequality:
N(‖ϕ‖h(N)(ψ(b)−ψ(a))αkραkΓk(α+k)+ua)−1>1. |
Then the proposed problem (1.1) has at least one solution on J.
Proof. Assume that Q is given by (3.3). The technique is done in three phases.
Step (i). We show that Q maps bounded sets (balls) into bounded sets in C([a,b],R).
Assume that Br2={u∈C([a,b],R):‖u‖≤r2} is a bounded ball in C([a,b],R) and ‖ϕ‖=supτ∈[a,b]|ϕ(τ)|. By applying the fact of 0≤eρ−1kρ(ψ(τ)−ψ(a))≤1 and (A2), for every τ∈[a,b], which yields that
|(Qu)(τ)‖≤uaeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)|f(s,u(s))|ds≤ua+1ραkkΓk(α)∫τa(ψ(τ)−ψ(s))αk−1ψ′(s)ϕ(s)h(‖u‖)ds≤ua+(ψ(b)−ψ(a))αkραkΓk(α+k)‖ϕ‖h(r2):=K. |
Step (ii). We show that Q maps bounded sets into equicontinuous sets of C([a,b],R).
Assume that τ1, τ2∈[a,b] with τ1<τ2, and u∈Br2. Hence
|(Qu)(τ2)−(Qu)(τ1)|≤ua|eρ−1kρ(ψ(τ2)−ψ(a))−eρ−1kρ(ψ(τ1)−ψ(a))|+|1ραkkΓk(α)∫τ2aρkΨαk−1ψ(τ2,s)ψ′(s)f(s,u(s))ds−1ραkkΓk(α)∫τ1aρkΨαk−1ψ(τ1,s)ψ′(s)f(s,u(s))ds|≤ua|eρ−1kρ(ψ(τ2)−ψ(a))−eρ−1kρ(ψ(τ1)−ψ(a))|+1ραkkΓk(α)∫τ2τ1ρkΨαk−1ψ(τ2,s)ψ′(s)|f(s,u(s))|ds+1ραkkΓk(α)∫τ1a|ρkΨαk−1ψ(τ2,s)−ρkΨαk−1ψ(τ1,s)|ψ′(s)|f(s,u(s))|ds≤ua|eρ−1kρ(ψ(τ2)−ψ(a))−eρ−1kρ(ψ(τ1)−ψ(a))|+‖ϕ‖h(r2)ραkΓk(α+k)(ψ(τ2)−ψ(τ1))αk+‖ϕ‖h(r2)ραkΓk(α+k)|(ψ(τ2)−ψ(a))αk−(ψ(τ2)−ψ(τ1))αk−(ψ(τ1)−ψ(a))αk|. | (3.5) |
Clearly, the right-hand side of (3.5) tends to zero independently of u∈Br2 as τ2→τ1. Then, by utilizing the Arzelá-Ascoli theorem, we get that Q:C([a,b],R)→C([a,b],R) is completely continuous.
Step (iii). We show that there exists an open set D⊆C([a,b],R) with u≠λQ(u), 0<λ<1 and u∈∂D.
Assume that u∈C([a,b],R) is the solution of u=λQu, 0<λ<1. For any τ∈[a,b], we obtain that
|u(τ)|=|λ(Qu)(τ)|≤ua+(ψ(b)−ψ(a))αkραkΓk(α+k)‖ϕ‖h(‖u‖). |
Then
‖u‖≤ua+(ψ(b)−ψ(a))αkραkΓk(α+k)‖ϕ‖h(‖u‖). |
Consequently,
‖u‖(‖ϕ‖h(‖u‖)(ψ(b)−ψ(a))αkραkΓk(α+k)+ua)−1≤1. |
In view of (A3), there exists N such that ‖u‖≠N. Let us take D={u∈C([a,b],R):‖u‖<N} and U=D∩Br2.
Finally, Q:¯U→C([a,b],R) is continuous and completely continuous. Using the choice of U, there is no 0∈∂U so that u=λ(Qu) for some 0<λ<1. Hence, by [42], we obtain that Q has a fixed point u∈¯U, which is a solution of the problem (1.1).
Next, we establish various of UH–ML stability results for the proposed problem (1.1).
Definition 3.4. The proposed problem (1.1) is called UH-ML stable, if there is a number Cf∈R+ so that for any ϵ>0 and w∈C([a,b],R+) of
|Ca,kDα,ρ;ψw(τ)−f(τ,w(τ))|≤ϵ,τ∈[a,b], | (3.6) |
there exists u∈C([a,b],R+) of the proposed problem (1.1) under the assumption
|w(τ)−u(τ)|≤CfϵEk,α,k(κf(ψ(τ)−ψ(s))αk),κf≥0,τ∈[a,b]. | (3.7) |
Definition 3.5. The proposed problem (1.1) is called generalized UH-ML stable, if there is a function Gf∈C(R+,R+) via Gf(0)=0, such that for any ϵ>0 and w∈C([a,b],R+) of
|Ca,kDα,ρ;ψw(τ)−f(τ,w(τ))|≤Gf(τ),τ∈[a,b], | (3.8) |
there exists u∈C([a,b],R+) of the proposed problem (1.1) under the assumption
|w(τ)−u(τ)|≤Gf(ϵ)Ek,α,k(κf(ψ(τ)−ψ(s))αk),κf≥0,τ∈[a,b]. | (3.9) |
Definition 3.6. The proposed problem (1.1) is called Ulam-Hyers-Rassias-Mittag-Leffler (UHR-ML) stable with respect to another function Φ(τ), if there is a number CfΦ∈R+ so that for any ϵ>0 and w∈C([a,b],R+) of
|Ca,kDα,ρ;ψw(τ)−f(τ,w(τ))|≤ϵΦ(τ),τ∈[a,b], | (3.10) |
there exists u∈C([a,b],R+) of the proposed problem (1.1) under the assumption
|w(τ)−u(τ)|≤CfΦϵΦ(τ)Ek,α,k(κfΦ(ψ(τ)−ψ(s))αk),κfΦ≥0,τ∈[a,b]. | (3.11) |
Definition 3.7. The proposed problem (1.1) is called generalized UHR-ML stable concerning function Φ(τ) so that for any ϵ>0 and w∈C([a,b],R+) of
|Ca,kDα,ρ;ψw(τ)−f(τ,w(τ))|≤Φ(τ),τ∈[a,b], | (3.12) |
there exists u∈C([a,b],R+) of the proposed problem (1.1) under the assumption
|w(τ)−u(τ)|≤CfΦΦ(τ)Ek,α,k(κfΦ(ψ(τ)−ψ(s))αk),κfΦ≥0,τ∈[a.b]. | (3.13) |
Remark 3.8. Assume that w∈C([a,b],R) is the solution of (3.6) if and only if there is uw∈C([a,b],R), which depends on w, so that (i). |uw(τ)|≤ϵ, τ∈[a,b]; (ii). Ca,kDα,ρ;ψw(τ)=f(τ,w(τ))+uw(τ), τ∈[a,b].
Remark 3.9. Assume that w∈C([a,b],R) is the solution of (3.10) if and only if there is vw∈C([a,b],R), that is depends on w, such that (i). |vw(τ)|≤ϵΦ(τ), τ∈[a,b]; (ii). Ca,kDα,ρ;ψz(τ)=f(τ,w(τ))+vw(τ), τ∈[a,b].
Theorem 3.10. Suppose that f∈C([a,b]×R,R), (A1), and (3.2) hold. Then, the proposed problem (1.1) is UH–ML stable and consequently generalized UH-ML stable on [a,b].
Proof. Let ϵ>0 and w∈C([a,b],R) be a solution of (3.6). From (ii) in Remark 3.8, we have
{Ca,kDα,ρ;ψw(τ)=f(τ,w(τ))+uw(τ),τ∈(a,b],w(a)=wa,wa∈R. | (3.14) |
Applying Lemma 3.1, the solution of the problem (3.14) is
w(τ)=waeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,w(s))ds+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)uw(s)ds. | (3.15) |
Let u∈C([a,b],R) be a solution of (1.1). Hence,
u(τ)=uaeρ−1kρ(ψ(τ)−ψ(a))+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)f(s,u(s))ds. | (3.16) |
By applying the property |u−v|≤|u|+|v|, τ∈J with (3.15)-(3.16), one has
|w(τ)−u(τ)|≤1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)|f(s,w(s))−f(s,u(s))|ds+1ραkkΓk(α)∫τaρkΨαk−1ψ(τ,s)ψ′(s)|uw(s)|ds. | (3.17) |
By using the property 0 \leq e^{{\frac{\rho -1}{k\rho}(\psi(\tau)-\psi(a))}} \leq 1 , (\mathcal{H}_{1}) , and (i) in Remark 3.8, the inequality (3.17) can be written
\begin{equation} \left\vert w(\tau) - u(\tau) \right\vert \leq \frac{\left(\psi(b)-\psi(a)\right)^{\frac{\alpha}{k}}}{\rho^{\frac{\alpha}{k}}\Gamma_{k}(\alpha+k)} \epsilon + \frac{\Gamma_{k}(\alpha)}{k}\left(\frac{k\mathcal{L}}{\Gamma_{k}(\alpha)}\right) {_{a,k}^{}}\mathcal{I}^{\alpha,\rho;\psi} \left\vert w(\tau) - u(\tau) \right\vert. \end{equation} | (3.18) |
From Theorem 2.5 and Corollary 2.7, which implies that
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \frac{\left(\psi(b)-\psi(a)\right)^{\frac{\alpha}{k}}}{\rho^{\frac{\alpha}{k}}\Gamma_{k}(\alpha+k)} \epsilon \mathbb{E}_{k,\alpha,k}\left(\mathcal{L}\rho^{-\frac{\alpha}{k}}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
By setting a constant \mathfrak{C}_{f} : = \frac{\left(\psi(b)-\psi(a)\right)^{\frac{\alpha}{k}}}{\rho^{\frac{\alpha}{k}}\Gamma_{k}(\alpha+k)} and \kappa_{f} : = \mathcal{L}\rho^{-\frac{\alpha}{k}} , which yields that
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \mathfrak{C}_{f} \epsilon \mathbb{E}_{k,\alpha,k}\left(\kappa_{f}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
Hence, the proposed problem (1.1) is UH–ML stable. In addition, by setting \mathcal{G}_{f}(\epsilon) = \mathfrak{C}_{f} \epsilon under the condition \mathcal{G}_{f}(0) = 0 , then
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \mathcal{G}_{f}(\epsilon) \mathbb{E}_{k,\alpha,k}\left(\kappa_{f}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
Therefore, the solution of the proposed problem (1.1) is generalized UH-ML stable.
Now, we give the required assumption that is used in Theorem 3.11.
(\mathcal{P}_{1}) Suppose that \Phi \in \mathcal{C}([a, b], \mathbb{R}) is a non-decreasing function. There exists \chi_{\Phi} > 0 so that
\begin{equation*} \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\displaystyle{\int}_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) \Phi(s) ds \leq \chi_{\Phi} \Phi(\tau),\quad \tau \in [a,b]. \end{equation*} |
Theorem 3.11. Suppose that f \in \mathcal{C}([a, b]\times\mathbb{R}, \mathbb{R}) , (\mathcal{A}_{1}) , and (3.2) hold. Then, the proposed problem (1.1) is UHR-ML stable and consequently generalized UHR-ML stable on [a, b] .
Proof. Let \epsilon > 0 and w \in \mathcal{C}([a, b], \mathbb{R}) be the solution of (3.6). From (ii) in Remark 3.9, it follows form
\begin{equation} \left\{ \begin{array}{l} {}_{a,k}^{{C}}\mathfrak{D}^{\alpha,\rho;\psi} w(\tau) = f(\tau,w(\tau)) + v_{w}(\tau), \quad \tau \in (a, b],\\ w(a) = w_{a},\quad w_{a} \in \mathbb{R}. \end{array} \right. \end{equation} | (3.19) |
From Lemma 3.1, the solution of the problem (3.19) is defined by
\begin{eqnarray} w(\tau) & = & w_{a} e^{{\frac{\rho -1}{k\rho}(\psi(\tau)-\psi(a))}} + \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\int_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) f(s,w(s)) ds\\ && + \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\int_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) v_{w}(s) ds. \end{eqnarray} | (3.20) |
Let u \in \mathcal{C}([a, b], \mathbb{R}) be the solution of (1.1). It follows that
\begin{equation} u(\tau) = u_{a} e^{{\frac{\rho -1}{k\rho}(\psi(\tau)-\psi(a))}} + \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\displaystyle{\int}_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) f(s,u(s)) ds. \end{equation} | (3.21) |
By applying the property \vert u - v \vert \leq \vert u \vert + \vert v \vert , for \tau \in [a, b] with (3.20)-(3.21), which yields that
\begin{eqnarray} \left\vert w(\tau) - u(\tau) \right\vert &\leq& \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\int_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) \left\vert f(s,w(s)) - f(s,u(s)) \right\vert ds\\ && + \frac{1}{\rho^{\frac{\alpha}{k}}k\Gamma_{k}(\alpha)}\int_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,s) \psi^{\prime}(s) \left\vert v_{w}(s) \right\vert ds. \end{eqnarray} | (3.22) |
By using 0 \leq e^{{\frac{\rho -1}{k\rho}(\psi(\tau)-\psi(a))}} \leq 1 , (\mathcal{A}_{1}) , (\mathcal{P}_{1}) , and (i) in Remark 3.9, the inequality (3.22) can be written
\begin{equation} \left\vert w(\tau) - u(\tau) \right\vert \leq \epsilon \chi_{\Phi} \Phi(\tau) + \frac{\Gamma_{k}(\alpha)}{k}\left(\frac{k\mathcal{L}}{\Gamma_{k}(\alpha)}\right) {_{a,k}^{}}\mathcal{I}^{\alpha,\rho;\psi} \left\vert w(\tau) - u(\tau) \right\vert. \end{equation} | (3.23) |
Applying Theorem 2.5 and Corollary 2.7, which implies that
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \epsilon \chi_{\Phi} \Phi(\tau) \mathbb{E}_{k,\alpha,k}\left(\mathcal{L}\rho^{-\frac{\alpha}{k}}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
By choosing \mathfrak{C}_{f_{\Phi}} : = \chi_{\Phi} and \kappa_{f} : = \mathcal{L}\rho^{-\frac{\alpha}{k}} , then
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \mathfrak{C}_{f_{\Phi}} \epsilon \Phi(\tau) \mathbb{E}_{k,\alpha,k}\left(\kappa_{f}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
Thus, the proposed problem (1.1) is UHR-ML stable. Additionally, by setting \epsilon = 1 , one has
\begin{equation*} \left\vert w(\tau) - u(\tau) \right\vert \leq \mathfrak{C}_{f_{\Phi}} \Phi(\tau) \mathbb{E}_{k,\alpha,k}\left(\kappa_{f}\left(\psi(\tau)-\psi(s)\right)^{\frac{\alpha}{k}} \right). \end{equation*} |
Then, the solution of the proposed problem (1.1) is generalized UHR-ML stable.
Now, we construct an approximation form for the (k, \psi) -PFDO in Caputo's sense under \alpha\in(0, 1] of a function u(\tau) . We may generate a sequence of N+1 equations with N+1 conditions for a given (k, \psi) proportional fractional Cauchy-type problem under this tool. A sequence (u_{N}) of solutions to such systems eventually leads to the solution of the proposed problem.
Theorem 4.1. Assume that N = 1, 2, \ldots, and u \in \mathcal{AC}^{2}([a, b], \mathbb{R}) . Let
\begin{eqnarray} A_{N} & = & \frac{1}{\rho^{1-\frac{\alpha}{k}} \Gamma_{k}(2k-\alpha)} \sum\limits_{i = 0}^{N}\frac{\Gamma\left(i+\frac{\alpha}{k}-1\right)}{i! \Gamma\left(\frac{\alpha}{k}-1\right)}, \end{eqnarray} | (4.1) |
\begin{eqnarray} B_{N, i} & = & \frac{\Gamma\left(i+\frac{\alpha}{k}-1\right)}{\rho^{1-\frac{\alpha}{k}} (i-1)! \Gamma_{k}(2k-\alpha) \Gamma\left(\frac{\alpha}{k}-1\right)},\quad i = 1, 2, \ldots, N, \end{eqnarray} | (4.2) |
and let \mathcal{V}_{i}:[a, b]\to\mathbb{R} be a function, which is given by
\begin{equation} \mathcal{V}_{i}(\tau) = \displaystyle{\int}_{a}^{\tau}{{\left( \psi (s)-\psi (a) \right)}^{i-1}}{\psi }^{\prime}(s){{e}^{\frac{1-\rho }{k\rho }\psi (s)}} {}_{k}{{\mathfrak{D}}^{\rho ;\psi }}u(s) ds,\quad i = 1, 2, \ldots, N. \end{equation} | (4.3) |
Hence,
\begin{eqnarray} {}_{a,k}^{C}{{\mathfrak{D}}^{\alpha ,\rho ;\psi }} u(\tau) &\approx& {A_{N}}{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}} {}_{k}{{\mathfrak{D}}^{\rho ;\psi }}u(\tau) \\ && -{e^{\frac{\rho -1}{k\rho }\psi (\tau)}}\sum\limits_{i = 1}^{N}{{{B}_{N,i}}{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}-i}}{{\mathcal{V}}_{i}}(\tau)} + E_{tr}(\tau), \end{eqnarray} | (4.4) |
where {}_{k}{{\mathfrak{D}}^{\rho; \psi }}u(\tau) = (1-\rho)u(\tau) + k\rho\frac{u^{\prime}(\tau)}{\psi^{\prime}(\tau)} and \lim\limits_{N\to\infty}{E}_{tr}(\tau) = 0 , \tau \in [a, b] .
Proof. By using Definition 2.3 and {}_{k}{{\mathfrak{D}}^{\rho; \psi }}u(\tau) = (1-\rho)u(\tau) + k\rho\frac{u^{\prime}(\tau)}{\psi^{\prime}(\tau)} , for any \alpha \in (0, 1] , we have
\begin{eqnarray} &&{}_{a,k}^{C}{{\mathfrak{D}}^{\alpha,\rho ;\psi }} u(\tau)\\ & = & \frac{1}{k\rho^{1-\frac{1}{k}}\Gamma_{k}(k-\alpha )}\int_{a}^{\tau} {_{k}^{\rho}}\Psi_{\psi}^{-\frac{\alpha}{k}}(\tau,s) {\psi }^{\prime}(s)\Big({{}_{k}}{{\mathfrak{D}}^{\rho ;\psi }}u(s)\Big)ds\\ & = & \frac{1}{k{{\rho }^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(k-\alpha )}\int_{{{a}^{}}}^{\tau}{{{\left( \psi (\tau)-\psi (s) \right)}^{-\frac{\alpha}{k}}}}{\psi }^{\prime}(s) e^{\frac{\rho -1}{k\rho }\left( \psi (\tau)-\psi (s) \right)} \left[(1-\rho )u(s) + k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right]ds. \end{eqnarray} | (4.5) |
Changing new variables
\begin{equation*} x(s) = {e^{\frac{\rho-1}{k\rho}\left(\psi(\tau)-\psi(s)\right)}}\left[ (1-\rho )u(s)+k\rho\frac{{u}^{\prime}(s)}{{\psi}^{\prime}(s)}\right]\quad\text{and}\quad y^{\prime}(s) = \left(\psi(\tau)-\psi(s)\right)^{-\frac{\alpha}{k}} {\psi}^{\prime}(s), \end{equation*} |
and by helping the integrating by parts technique, the equation (4.5) can be solved as
\begin{eqnarray} &&{}_{a,k}^{C}{{\mathfrak{D}}^{\alpha ,\rho;\psi}} u(\tau)\\ & = & \frac{{_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,a) }{{\rho^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )} \left[ (1-\rho )u(a)+k\rho \frac{{u}^{\prime}(a)}{{\psi }^{\prime}(a)} \right] \\ && +\frac{1}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma}_{k}}(2k-\alpha )}\int_{{{a}^{}}}^{\tau}{\left( \psi (\tau)-\psi (s) \right)^{1-\frac{\alpha }{k}}} \frac{d}{ds}\left( {e^{\frac{\rho -1}{k\rho }\left( \psi (\tau)-\psi (s) \right)}}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right)ds. \end{eqnarray} | (4.6) |
Applying the Newton's generalized binomial theorem, it follows that
\begin{eqnarray} (\psi (\tau)-\psi (s))^{1-\frac{\alpha }{k}} & = & {{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}{{\left( 1-\frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{1-\frac{\alpha }{k}}} \\ & = & (\psi (\tau)-\psi (a))^{1-\frac{\alpha }{k}}\sum\limits_{i = 0}^{N }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i! \Gamma \left( \frac{\alpha }{k}-1 \right)}{{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{i}}} \\ && + (\psi (\tau)-\psi (a))^{1-\frac{\alpha }{k}}\sum\limits_{i = N+1}^{\infty }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i! \Gamma \left( \frac{\alpha }{k}-1 \right)}{{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{i}}}. \end{eqnarray} | (4.7) |
Inserting (4.7) into (4.6), we have
\begin{eqnarray} {}_{a,k}^{C}{{\mathfrak{D}}^{\alpha ,\rho ;\psi }} u(\tau) & = & \frac{{_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}-1}(\tau,a) }{{{\rho}^{1-\frac{\alpha}{k}}}{{\Gamma}_{k}}(2k-\alpha)} \left[ (1-\rho )u(a)+k\rho \frac{u^{\prime}(a)}{{\psi }^{\prime}(a)} \right] \\ && +\frac{1}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )}\int_{{a}}^{\tau}{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}\sum\limits_{i = 0}^{N}\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i!\Gamma \left( \frac{\alpha }{k}-1 \right)}{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)^{i}} \\ && \times \frac{d}{ds}\left( e^{\frac{\rho -1}{k\rho }\left( \psi (\tau)-\psi (s) \right)}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right) ds \\ && + \frac{1}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma}_{k}}(2k-\alpha )}\int_{{{a}^{}}}^{\tau}{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}\sum\limits_{i = N+1}^{\infty }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i! \Gamma \left( \frac{\alpha }{k}-1 \right)}{{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{i}}} \\ && \times \frac{d}{ds} \left( {e^{\frac{\rho -1}{k\rho }\left( \psi (\tau)-\psi (s) \right)}}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right) ds. \end{eqnarray} | (4.8) |
where
\begin{eqnarray} E_{tr}(\tau) & = & \frac{{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma}_{k}}(2k-\alpha )}\int_{{{a}^{}}}^{\tau}R_{N}(s) \frac{d}{ds} \left( {e^{\frac{\rho -1}{k\rho }\left( \psi (\tau)-\psi (s) \right)}}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right) ds, \end{eqnarray} | (4.9) |
\begin{eqnarray} R_{N}(s) & = & \sum\limits_{i = N+1}^{\infty }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i! \Gamma \left( \frac{\alpha }{k}-1 \right)}{{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{i}}}. \end{eqnarray} | (4.10) |
Then, we have
\begin{eqnarray} {}_{a,k}^{C}{{\mathfrak{D}}^{\alpha ,\rho ;\psi }} u(\tau) & = & \frac{(\psi (\tau)-\psi (a))^{1-\frac{\alpha }{k}}}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{{\psi }^{\prime}(\tau)} \right] \\ && + \frac{{{e}^{\frac{\rho -1}{k\rho }\psi (\tau)}}{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}}{{{\rho }^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )}\sum\limits_{i = 1}^{N}{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i!\Gamma \left( \frac{\alpha }{k}-1 \right){{\left( \psi (\tau)-\psi (a) \right)}^{i}}}} \\ && \times \int_{{{a}^{}}}^{\tau}{{\left( \psi (s)-\psi (a) \right)}^{i}} \frac{d}{ds}\left( {e^{\frac{1-\rho }{k\rho }\psi (s)}}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right) ds + {{E}_{tr}(\tau)}. \end{eqnarray} |
By using integrating by parts with changing variable x(s) = {{\left(\psi (s)-\psi (a) \right)}^{i}} and y^{\prime}(s) = \frac{d}{ds}({e^{\frac{1-\rho }{k\rho }\psi (s)}}[ (1-\rho)u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} ]) , which yields that
\begin{eqnarray} {}_{a,k}^{C}{{\mathfrak{D}}^{\alpha ,\rho ;\psi }} u(\tau) & = & \frac{{{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}}}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{{\psi }^{\prime}(\tau)} \right]\left( 1 + \sum\limits_{i = 1}^{N}{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i!\Gamma \left( \frac{\alpha }{k}-1 \right)}} \right) \\ && -\frac{{e^{\frac{\rho -1}{k\rho }\psi (\tau)}}}{{{\rho }^{1-\frac{\alpha }{k}}}{{\Gamma }_{k}}(2k-\alpha )}\sum\limits_{i = 1}^{N}{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right){{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}-i}}}{(i-1)!\Gamma \left( \frac{\alpha }{k}-1 \right)}}\\ && \times \int_{a}^{\tau}{\left( \psi (s)-\psi (a) \right)^{i-1}}{\psi }^{\prime}(s){{e}^{\frac{1-\rho }{k\rho }\psi (s)}}\left[ (1-\rho )u(s)+k\rho \frac{u^{\prime}(s)}{\psi^{\prime}(s)} \right]ds + {E_{tr}(\tau)}, \end{eqnarray} |
Next, we study the term E_{tr}(\tau) , which is the error caused by truncation error. Finally, we show that {E}_{tr}(\tau) \to 0 as {N\to\infty} , \tau \in [a, b] , and to prove this, we provide an upper bound for the error. Since 0 \leq \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \leq 1 for s , \tau \in \mathcal{J} , The equation (4.10) can be computed that
\begin{equation*} {\left\vert R_{N}(s) \right\vert \leq} \sum\limits_{i = N+1}^{\infty }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i!\Gamma \left( \frac{\alpha }{k}-1 \right)}{{\left( \frac{\psi (s)-\psi (a)}{\psi (\tau)-\psi (a)} \right)}^{i}}} \leq \sum\limits_{i = N+1}^{\infty }{\frac{\Gamma \left( i+\frac{\alpha }{k}-1 \right)}{i!\Gamma \left( \frac{\alpha }{k}-1 \right)}} \leq \sum\limits_{i = N+1}^{\infty }{\frac{{{i}^{\frac{\alpha }{k}-2}}}{\Gamma \left( \frac{\alpha }{k}-1 \right)}} \leq \int\limits_{N}^{\infty }{\frac{{{s}^{\frac{\alpha }{k}-2}}}{\Gamma \left( \frac{\alpha }{k}-1 \right)}ds}, \end{equation*} |
which implies that
\begin{equation} {\left\vert R_{N}(s) \right\vert \leq} \frac{1}{{{N}^{1-\frac{\alpha }{k}}}\left( 1-\frac{\alpha }{k} \right)\Gamma \left( \frac{\alpha }{k}-1 \right)}. \end{equation} | (4.11) |
Taking \mathcal{M}(\tau) = \max\limits_{s\in [a, \tau]} \left\vert \frac{d}{ds} \left({e^{\frac{\rho -1}{k\rho }\left(\psi (\tau)-\psi (s) \right)}}\left[ (1-\rho)u(s)+k\rho \frac{u^{\prime}(s)}{{\psi }^{\prime}(s)} \right] \right) \right\vert and (4.11) into (4.9), we obtain the following upper bound
\begin{eqnarray} \left\vert E_{tr}(\tau) \right\vert \leq \frac{\mathcal{M}(\tau){{(\psi (\tau)-\psi (a))}^{1-\frac{\alpha }{k}}} (\tau - a)}{{{\rho}^{1-\frac{\alpha }{k}}}{{\Gamma}_{k}}(2k-\alpha ) {{N}^{1-\frac{\alpha }{k}}}\left( 1-\frac{\alpha }{k} \right)\Gamma \left( \frac{\alpha }{k}-1 \right)}. \end{eqnarray} | (4.12) |
The right-sided of (4.12) tends to zero as N \to \infty , \tau \in [a, b] . The proof is done.
From the property (iv) in Lemma 2.4, the integral equation (3.1) in Lemma 3.1 can be reformed as
\begin{equation} u(\tau) = u_{a} e^{{\frac{\rho -1}{k\rho}(\psi(\tau)-\psi(a))}} + {}_{a,k}^{C}\mathfrak{D}^{k-\alpha,\rho;\psi}\left(\frac{1}{\rho k}\displaystyle{\int}_{a}^{\tau}e^{\frac{\rho-1}{k\rho}(\psi(\tau)-\psi(s))} \psi^{\prime}(s) f(s,u(s)) ds\right). \end{equation} | (4.13) |
By applying Theorem 4.1, the equation (4.13) can be re-written as
\begin{eqnarray} u(s) & = & u_{a} e^{{\frac{\rho -1}{k\rho}(\psi(s)-\psi(a))}} + A^{*}_{N}{{(\psi (s)-\psi (a))}^{\frac{\alpha }{k}}}f(s,u(s)) \\ && -{e^{\frac{\rho -1}{k\rho }\psi (s)}}\sum\limits_{i = 1}^{N}{B^{*}_{N,i}{{(\psi (s)-\psi (a))}^{\frac{\alpha }{k}-i}}{{\mathcal{V}}^{*}_{i}}(s)} +{E^{*}_{tr}(s)}, \end{eqnarray} | (4.14) |
where
\begin{eqnarray*} A^{*}_{N} & = &\frac{1}{{{\rho }^{\frac{\alpha }{k}}}{{\Gamma }_{k}}(k+\alpha )} \sum\limits_{i = 0}^{N}{\frac{\Gamma \left( i-\frac{\alpha }{k} \right)}{i! \Gamma \left( -\frac{\alpha }{k} \right)}} , \notag\\ B^{*}_{N,i} & = &\frac{\Gamma \left( i-\frac{\alpha }{k} \right)}{{{\rho }^{\frac{\alpha }{k}}}{{(i-1)!\Gamma }_{k}}(k+\alpha )\Gamma \left( -\frac{\alpha }{k} \right)},\quad i = 1,2,\ldots,N, \notag\\ \mathcal{V}^{*}_{i}(s) & = & \int_{a}^{s}{{{\left( \psi (r)-\psi (a) \right)}^{i-1}}{\psi }^{\prime}(r){{e}^{\frac{1-\rho }{k\rho }\psi (r)}}f(r,u(r))dr}, \notag\\ \left\vert {E}^{*}_{tr}(s) \right\vert &\leq& \frac{\mathcal{M}(s)(s-a){{(\psi (s)-\psi (a))}^{\frac{\alpha }{k}}}}{{N^{\frac{\alpha }{k}}}\left( \frac{\alpha }{k} \right)\Gamma \left( -\frac{\alpha }{k} \right){{\rho }^{\frac{\alpha }{k}}}{{\Gamma }_{k}}(k+\alpha )} . \end{eqnarray*} |
In order to obtain the formula of an approximated solution, u_{N}(\tau) , we truncate the formula up to order N , getting
\begin{eqnarray} u_{N}(s) & = & u_{a} e^{{\frac{\rho -1}{k\rho}(\psi(s)-\psi(a))}} + {{A}^{*}_{N}}{{(\psi (s)-\psi (a))}^{\frac{\alpha }{k}}}f(s,u_{N}(s)) \\ && -{{e}^{\frac{\rho -1}{k\rho }\psi (s)}}\sum\limits_{i = 1}^{N}{{{B}^{*}_{N,i}}{{(\psi (s)-\psi (a))}^{\frac{\alpha }{k}-i}}{{\mathcal{V}}^{*}_{i,N}}(s)}, \end{eqnarray} | (4.15) |
where
\begin{equation*} \mathcal{V}^{*}_{i,N}(s) = \displaystyle{\int}_{a}^{s}{{{\left( \psi (r)-\psi (a) \right)}^{i-1}}{\psi }^{\prime}(r){{e}^{\frac{1-\rho }{k\rho }\psi (r)}}f(r,u_{N}(r))dr},\quad i = 1,2,\ldots, N. \end{equation*} |
Note that u_{N}(a) = u(a) .
Theorem 4.2. Assume that f \in \mathcal{C}([a, b]\times\mathbb{R}, \mathbb{R}) which verifies (\mathcal{A}_{1}) , and assume that u and u_{N} as in (4.14) and (4.15), respectively, for N \in\mathbb{N} . Also, suppose that b \in \mathbb{R} is a real in the open interval
\begin{equation} \psi(a) < \psi(b) < \psi(a) + \rho k \left( \frac{ \Gamma(\frac{\alpha}{k}+1)}{\mathcal{L}} \right)^{\frac{k}{\alpha}} . \end{equation} | (4.16) |
Then, u_{N}(\tau) \to u(\tau) as N \to \infty , for all \tau \in [a, b] .
Proof. From the equations (4.14) and (4.15), for any \tau\in [a, b] , we have
\begin{eqnarray} \vert x_{N}(\tau)-x(\tau)\vert &\leq& \vert A^{*}_{N} \vert (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}} \vert f(\tau,x_{N}(\tau))-f(\tau,x(\tau)) \vert \\ && + e^{\frac{\rho -1}{k\rho }\psi (\tau)}\sum\limits_{i = 1}^{N} \vert B^{*}_{N,i} \vert (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}-i} \vert \mathcal{V}^{*}_{i,N}(\tau) - \mathcal{V}^{*}_{i}(\tau) \vert + {\vert {E}^{*}_{tr}(\tau) \vert}. \end{eqnarray} | (4.17) |
Now, we define \delta_{u_{N}}: = \max\limits_{\tau\in[a, b]} \vert u_{N}(\tau)-u(\tau)\vert . Then, we have
\begin{eqnarray} \vert f(\tau,u_{N}(\tau))-f(\tau,u(\tau)) \vert &\leq& \mathcal{L} \vert u_{N}(\tau)-u(\tau)\vert \leq \mathcal{L}\delta_{u_{N}}, \end{eqnarray} | (4.18) |
\begin{eqnarray} \vert \mathcal{V}^{*}_{i,N}(\tau) - \mathcal{V}^{*}_{i}(\tau) \vert &\leq& \int_{a}^{\tau} \left(\psi(s)-\psi(a)\right)^{i-1}\psi^{\prime}(s) e^{\frac{1-\rho }{k\rho}\psi (s)} \vert f(s,u_{N}(s))- f(s,u(s)) \vert ds \\ &\leq& \mathcal{L}\delta_{u_{N}} \int_{a}^{\tau} \left(\psi(s)-\psi(a)\right)^{i-1}\psi^{\prime}(s) e^{\frac{1-\rho }{k\rho}\psi (s)} ds \\ &\leq& \frac{\mathcal{L}\delta_{u_{N}}}{i} \left(\psi(\tau)-\psi(a)\right)^{i}. \end{eqnarray} | (4.19) |
By using the formula (3) in [43], we have
\begin{equation} \vert A^{*}_{N} \vert = \frac{1}{{{\rho }^{\frac{\alpha }{k}}}{{\Gamma }_{k}}(k+\alpha )}\left\vert \sum\limits_{i = 0}^{N}{\frac{\Gamma \left( i-\frac{\alpha }{k} \right)}{\Gamma\left(-\frac{\alpha}{k} \right)i!}} \right\vert = \frac{\rho^{-\frac{\alpha }{k}}\Gamma(N+1-\frac{\alpha}{k})}{\Gamma_{k}(k+\alpha) \vert\Gamma\left(-\frac{\alpha}{k} \right)\vert \frac{\alpha}{k}\Gamma(N+1)} \leq \frac{k^{1-\frac{\alpha}{k}}\Gamma(N+1-\frac{\alpha}{k})}{\alpha \rho^{\frac{\alpha}{k}} \pi \Gamma(N+1)}, \end{equation} | (4.20) |
and using (4.19), it follows that
\begin{eqnarray} && \sum\limits_{i = 1}^{N} \vert B^{*}_{N,i} \vert (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}-i} \vert \mathcal{V}^{*}_{i,N}(\tau) - \mathcal{V}^{*}_{i}(\tau) \vert\\ &\leq& \frac{\mathcal{L}\delta_{x_{N}}(\psi (\tau)-\psi (a))^{\frac{\alpha }{k}}} {\rho^{\frac{\alpha }{k}} \Gamma_{k}(k+\alpha) \vert \Gamma \left( -\frac{\alpha }{k} \right)\vert} \left\vert \sum\limits_{i = 1}^{N} \frac{\Gamma\left(i-\frac{\alpha}{k}\right)}{i!} \right\vert \\ &\leq& \frac{\mathcal{L}\delta_{x_{N}}(\psi (\tau)-\psi (a))^{\frac{\alpha }{k}}} {\rho^{\frac{\alpha }{k}} \Gamma_{k}(k+\alpha) \vert \Gamma \left( -\frac{\alpha }{k} \right)\vert} \left( \frac{\Gamma(N+1-\frac{\alpha}{k})}{\frac{\alpha}{k}\Gamma(N+1)} + \left\vert \Gamma\left(-\frac{\alpha}{k}\right)\right\vert \right) \\ &\leq& \frac{\mathcal{L}\delta_{x_{N}}}{\rho^{\frac{\alpha }{k}}} (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{k^{1-\frac{\alpha}{k}}\Gamma(N+1-\frac{\alpha}{k})}{\alpha \pi \Gamma(N+1)} + \frac{k^{1-\frac{\alpha}{k}}}{\alpha \Gamma(\frac{\alpha}{k})} \right). \end{eqnarray} | (4.21) |
Substitution (4.18), (4.20), and (4.21) into the inequality (4.17), for any \tau\in[a, b] , one has
\begin{eqnarray} \vert u_{N}(\tau)-u(\tau)\vert &\leq& \frac{\mathcal{L}\delta_{u_{N}}k^{1-\frac{\alpha}{k}}}{\rho^{\frac{\alpha }{k}} \alpha \pi } (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}} \frac{\Gamma(N+1-\frac{\alpha}{k})}{\Gamma(N+1)} + {\vert {E}^{*}_{tr}(\tau) \vert} \\ && + e^{\frac{\rho -1}{k\rho }\psi (\tau)}\frac{\mathcal{L}\delta_{u_{N}}}{\rho^{\frac{\alpha }{k}}} (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{k^{1-\frac{\alpha}{k}}\Gamma(N+1-\frac{\alpha}{k})}{\alpha \pi \Gamma(N+1)} + \frac{k^{1-\frac{\alpha}{k}}}{\alpha \Gamma(\frac{\alpha}{k})} \right) \\ &\leq& \frac{\mathcal{L}\delta_{u_{N}}k^{1-\frac{\alpha}{k}}}{\rho^{\frac{\alpha }{k}} \alpha } (\psi (\tau)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{2\Gamma(N+1-\frac{\alpha}{k})}{\pi \Gamma(N+1)} + \frac{1}{\Gamma(\frac{\alpha}{k})} \right) + {\vert {E}^{*}_{tr}(\tau) \vert}. \end{eqnarray} | (4.22) |
Applying the maximum into the inequality (4.22) over \tau\in[a, b] , we obtain that
\begin{equation} \delta_{u_{N}} \leq \frac{\mathcal{L}\delta_{u_{N}} k^{1-\frac{\alpha}{k}}}{\rho^{\frac{\alpha }{k}} \alpha } (\psi (b)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{2\Gamma(N+1-\frac{\alpha}{k})}{\pi \Gamma(N+1)} + \frac{1}{\Gamma(\frac{\alpha}{k})} \right) + \max\limits_{\tau\to[a,b]} {\vert {E}^{*}_{tr}(\tau) \vert}. \end{equation} | (4.23) |
Clearly, \lim\limits_{N\to\infty} {\vert {E}^{*}_{tr}(\tau) \vert} = 0 . In addition, by applying the Stirling's formula [44], we have
\begin{equation} \lim\limits_{N\to\infty} \frac{\Gamma(N+1-\frac{\alpha}{k})}{\Gamma(N+1)} = 0. \end{equation} | (4.24) |
Then, taking N \to \infty in (4.23), one has
\begin{equation*} \lim\limits_{N\to\infty} \delta_{u_{N}}\left( 1- \frac{\mathcal{L}(\psi (b)-\psi (a))^{\frac{\alpha }{k}}}{\rho^{\frac{\alpha }{k}} k^{\frac{\alpha}{k}} \Gamma(\frac{\alpha}{k}+1) } \right) \leq 0. \end{equation*} |
By the condition (4.16), we get that \delta_{u_{N}} \to 0 as N \to \infty . From (4.23), which implies that
\begin{equation*} \delta_{u_{N}} \left[1-\frac{\mathcal{L} k^{1-\frac{\alpha}{k}}}{\rho^{\frac{\alpha }{k}} \alpha } (\psi (b)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{2\Gamma(N+1-\frac{\alpha}{k})}{\pi \Gamma(N+1)} + \frac{1}{\Gamma(\frac{\alpha}{k})} \right) \right] \leq \max\limits_{\tau\to\mathcal{J}} \vert {E}^{*}_{tr}(\tau) \vert. \end{equation*} |
Applying the condition (4.16), we obtain the following result
\begin{equation*} 0 < \frac{\mathcal{L}k^{\frac{1-\alpha}{k}}}{\rho^{\frac{\alpha }{k}}\alpha}(\psi (b)-\psi (a))^{\frac{\alpha }{k}} < \Gamma\left(\frac{\alpha}{k}\right). \end{equation*} |
Then,
\begin{equation*} -\frac{2\Gamma(\frac{\alpha}{k})\Gamma(N+1-\frac{\alpha}{k})}{\pi \Gamma(N+1)} < 1-\frac{\mathcal{L} k^{1-\frac{\alpha}{k}}}{\rho^{\frac{\alpha }{k}} \alpha } (\psi (b)-\psi (a))^{\frac{\alpha }{k}} \left( \frac{2\Gamma(N+1-\frac{\alpha}{k})}{\pi \Gamma(N+1)} + \frac{1}{\Gamma(\frac{\alpha}{k})} \right) < 1. \end{equation*} |
For N sufficiently large number, we obtain the property (4.24). Therefore, there exists a function \mathcal{B} (depending on \alpha , \rho , k , \tau but independent of N ), such that
\begin{equation*} \delta_{u_{N}} \leq \mathcal{B}_{\psi}(\alpha, k, \rho, \tau) N^{-\frac{\alpha }{k}}, \end{equation*} |
where
\begin{equation} \mathcal{B}_{\psi}(\alpha, k, \rho, \tau) = \frac{\mathcal{M}(\tau){{(\psi (b)-\psi (a))}^{\frac{\alpha }{k}}} (b-a)}{\left( \frac{\alpha }{k} \right)\Gamma \left( -\frac{\alpha }{k} \right){{\Gamma }_{k}}(k+\alpha )} \left(\frac{\alpha \Gamma(\frac{\alpha}{k})} {\rho^{\frac{\alpha }{k}} \alpha \Gamma(\frac{\alpha}{k})-\mathcal{L} k^{1-\frac{\alpha}{k}}(\psi (b)-\psi (a))^{\frac{\alpha }{k}}} \right). \end{equation} | (4.25) |
Remark 4.3. Under Theorem 4.1 and Theorem 4.2, we have the following results:
(ⅰ) If k = 1 and \psi(\tau) = \tau then, Theorem 4.1 and Theorem 4.2 reduce to the results as in [20].
(ⅱ) If k > 1 , \rho = 1 , and \psi(\tau) = \tau^{\mu} , where \mu > 0 then, Theorem 4.1 and Theorem 4.2 reduce to the results as in [21].
(ⅲ) If k > 1 , \rho = 1 and \psi(\tau) = \tau then, Theorem 4.1 and Theorem 4.2 reduce to the results as in [22,23].
(ⅳ) If k > 1 , \rho = 1 , and \psi(\tau) = \tau^{\mu} , where \mu \to 0^{+} then, Theorem 4.1 and Theorem 4.2 reduce to the results as in [24].
This section provides two different numerical examples of applications to verify the theoretical results of our main results.
Example 5.1. Consider the following Cauchy-type problem under the (k, \psi) -Caputo-PFDO:
\begin{equation} \left\{ \begin{array}{l} {}_{0,k}^{{C}}\mathfrak{D}^{\alpha,\rho;\psi} u(\tau) = \frac{{\rho^{\frac{\alpha}{k}}}\Gamma_{k}\left( \omega + k \right)}{\Gamma_{k}\left( \omega+k-\alpha \right)} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega-\alpha}{k}}(\tau,0) + \frac{1}{2^{\tau+3}} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega}{k}}(\tau,0) - \frac{1}{2^{\tau+3}}u(\tau), \quad \tau\in[0,1], \\ u(0) = 0,\quad \alpha\in(0,1],\quad \rho \in (0,1],\quad k > 0. \end{array} \right. \end{equation} | (5.1) |
The exact solution of the problem (5.1) is provided by u(\tau) = {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega}{k}}(\tau, 0) , where \omega \in \mathbb{R} and \omega > k .
From the problem (5.1), we get
\begin{equation*} f(\tau, u(\tau)) = \frac{{\rho^{\frac{\alpha}{k}}}\Gamma_{k}\left( \omega + k \right)}{\Gamma_{k}\left( \omega+k-\alpha \right)} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega-\alpha}{k}}(\tau,0) + \frac{1}{2^{\tau+3}} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega}{k}}(\tau,0) - \frac{1}{2^{\tau+3}} u(\tau). \end{equation*} |
For every u , v \in \mathbb{R} , \tau \in [0, 1] , we get that
\begin{equation} \left\vert f(\tau, u(\tau)) - f(\tau, v(\tau)) \right\vert \leq \frac{1}{2^{\tau+3}} \left\vert u(\tau) - v(\tau) \right\vert \leq \frac{1}{8} \left\vert u(\tau) - v(\tau) \right\vert. \end{equation} | (5.2) |
It is noticed that the assumption (\mathcal{A}_{1}) is held with \mathcal{L} = 1/8 . By setting \alpha = 22/25 , \rho = 1/4 , k = 5/4 , \omega = 3/2 , and \psi(\tau) = \tau , then (\mathcal{L}\left(\psi(b)-\psi(a)\right)^{\frac{\alpha}{k}})/(\rho^{\frac{\alpha}{k}}\Gamma_{k}(\alpha+k)) \approx 0.311730513 < 1 . Since all assumptions in Theorem 3.2 are satisfied, the proposed problem (5.1) has a unique solution on [0, 1] . From Theorem 3.10, we can compute that
\begin{equation*} \mathfrak{C}_{f} : = \frac{\left(\psi(b)-\psi(a)\right)^{\frac{\alpha}{k}}}{\rho^{\frac{\alpha}{k}}\Gamma_{k}(\alpha+k)} \approx 2.493844103 > 0,\quad \text{and} \quad \kappa_{f} : = \mathcal{L}\rho^{-\frac{\alpha}{k}} \approx 0.331711285 \geq 0. \end{equation*} |
Then, the proposed problem (5.1) is UH-ML stable on [0, 1] . If we set \mathcal{G}_{f}(\epsilon) = \mathfrak{C}_{f} \epsilon under the condition \mathcal{G}_{f}(0) = 0 , then we obtain the proposed problem (5.1) is generalized UH-ML stable on [0, 1] . In addition, by setting a non-decreasing function \Phi(\tau) = {_{k}^{\rho}}\Psi_{\psi}^{\frac{3}{k}}(\tau, a) , we have
\begin{equation*} {}_{a,k}\mathcal{I}_{{{}^{}}}^{\alpha ,\rho ;\psi }\Phi(\tau) = \frac{\Gamma_{k}(3) {_{k}^{\rho}}\Psi_{\psi}^{\frac{3+\alpha}{k}-1}(\tau,a) }{\rho^{\frac{\alpha}{k}}\Gamma_{k}(3+\alpha) } = \frac{\Gamma_{k}(3) (\psi(\tau)-\psi(a))^{\frac{\alpha}{k}}}{\rho^{\frac{\alpha}{k}}\Gamma_{k}(3+\alpha)}\; \Phi(\tau). \end{equation*} |
This yields that \mathfrak{C}_{f_{\Phi}} : = \chi_{\Phi} = (\Gamma_{k}(3))/(\rho^{\frac{\alpha}{k}}\Gamma_{k}(3+\alpha))\; (\psi(\tau)-\psi(a))^{\frac{\alpha}{k}} \approx 1.276974061 > 0 and \kappa_{f_{\Phi}} : = \mathcal{L}\rho^{-\frac{\alpha}{k}} \approx 0.331711285 \geq 0 . Hence, the proposed problem (5.1) is UHR-ML stable on [0, 1] . If we take \epsilon = 1 , then the proposed problem (5.1) is generalized UHR-ML stable on [0, 1] .
To achieve the numerical approximation of the proposed problem (5.1), we apply (4.4) in Theorem 4.1. Then, the operator {}_{0, k}^{C}{\mathfrak{D}^{\alpha, \rho; \psi}}u(\tau) can be replaced as follows:
\begin{eqnarray} {}_{0,k}^{C}{\mathfrak{D}^{\alpha,\rho;\psi}}u(\tau) &\approx& A_{N}{{(\psi (\tau)-\psi (0))}^{1-\frac{\alpha }{k}}}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{\psi^{\prime} (\tau)} \right] \\ && - e^{\frac{\rho -1}{k\rho }\psi (\tau)}\sum\limits_{i = 1}^{N}B_{N,i}(\psi (\tau)-\psi (0))^{1-\frac{\alpha }{k}-i}\mathcal{V}_{i}(\tau), \end{eqnarray} | (5.3) |
where {A}_{N} , {B}_{N, i} , and \mathcal{V}_{i}(\tau) are given by (4.1), (4.2), and (4.3), respectively, in Theorem 4.1. Then
\begin{equation} \left\{ \begin{array}{l} \mathcal{V}^{\prime}_{i}(\tau) = \left( \psi (\tau)-\psi (0) \right)^{i-1}\psi^{\prime}(\tau){{e}^{\frac{1-\rho }{k\rho }\psi (\tau)}}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{\psi^{\prime}(\tau)} \right], \quad i = 1,2,\ldots,N,\\ \mathcal{V}_{i}(0) = 0. \end{array} \right. \end{equation} | (5.4) |
By applying (5.3) and (5.4), then
\begin{equation} \left\{ \begin{array}{l} A_{N}{{(\psi(\tau)-\psi (0))}^{1-\frac{\alpha }{k}}}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{\psi^{\prime} (\tau)} \right] \\ \qquad - e^{\frac{\rho -1}{k\rho }\psi (\tau)}\sum\limits_{i = 1}^{N}{B_{N,i}{{(\psi (\tau)-\psi (0))}^{1-\frac{\alpha }{k}-i}}{{\mathcal{V}}_{i}}(\tau)} \\ \qquad\qquad = \frac{{\rho^{\frac{\alpha }{k}}}{{\Gamma }_{k}}\left( \omega +k \right)}{{{\Gamma }_{k}}\left( \omega +k-\alpha \right)} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\omega-\alpha}{k}}(\tau,0) + \frac{1}{2^{\tau+3}} {_{k}^{\rho}}\Psi_{\psi}^{\frac{\alpha}{k}}(\tau,0) -\frac{1}{2^{\tau+3}}u(\tau), \quad \tau\in[0,1], \\ \mathcal{V}^{\prime}_{i}(\tau) = {{\left( \psi (\tau)-\psi (0) \right)}^{i-1}}{\psi }^{\prime}(\tau){{e}^{\frac{1-\rho }{k\rho }\psi (\tau)}}\left[ (1-\rho )u(\tau)+k\rho \frac{u^{\prime}(\tau)}{\psi^{\prime}(\tau)} \right], \quad i = 1,2,\ldots,N,\\ u(0) = 0, \qquad \mathcal{V}_{i}(0) = 0. \end{array} \right. \end{equation} | (5.5) |
The proposed system (5.5) can solve the numerical solution by applying ode45 in MATLAB software, which applies the explicit Runge-Kutta technique to achieve the approximated solution. The exact solution (continuous line) and numerical approximations (dot lines) of the proposed problem are compared with different parameters \rho = 1/4 , 1/3 , 1/2 , and k = 1/2 , 3/4 .
Figure 1 displays the behavior of the exact and approximate solutions under the fixed value \alpha = 9/10 , and vary values \rho = 1/4 , 1/3 , 1/2 , k = 1/2 , 3/4 , N = 2, 8, 16 for the proposed problem (5.1). While Figure 2 displays the behavior of the absolute error between the exact and approximate solutions under the fixed values \alpha = 9/10 , and vary values \rho = 1/4 , 1/3 , 1/2 , k = 1/2 , 3/4 , N = 2, 8, 16 for the proposed problem (5.1). Figures 1a–1b illustrate that as the value k increases, the approximate solution quickly converges to the exact solution. Moreover, the absolute error values rapidly decrease to zero, as can be seen from Figures 2a–2b. Conversely, Figures 1c–1d indicate that as the value \rho increases, the approximate solution converges to the exact solution more quickly and to greater values. While the absolute error values increase, as can be seen from the Figures 2c–2d.
Now, an application has been developed to support our theoretical results. We study a basic model for determining BALs that is characterized by two differential equations with real-data individuals. The BALs problem is explained by the following Cauchy-type system of the form:
\begin{equation} \left\{ \begin{array}{l} \frac{d\mathcal{S}}{d\tau} = - c_{1} \mathcal{S}(\tau),\qquad\qquad\,\, \mathcal{S}(0) = \mathcal{S}_{0}, \\ \frac{d\mathcal{B}}{d\tau} = c_{1} \mathcal{S}(\tau) - c_{2} \mathcal{B}(\tau),\quad \mathcal{B}(0) = \mathcal{B}_{0} = 0, \end{array} \right. \end{equation} | (5.6) |
where \mathcal{S}(\tau) and \mathcal{B}(\tau) denote the concentrations (mg/L) of alcohol in the stomach and blood of a human body at time \tau (min), respectively. The rate law constants i (min ^{-1} ) are given by c_{i} , i = 1, 2 , and \mathcal{S}_{0} , \mathcal{B}_{0} represent the subject's initial alcohol intake in the stomach and blood, respectively. The exact solution for the system (5.6) is fairly simple and may be determined using the Laplace transform approach
\begin{eqnarray} \mathcal{S}(\tau) & = & \mathcal{S}_{0} \exp(-c_{1} \tau), \end{eqnarray} | (5.7) |
\begin{eqnarray} \mathcal{B}(\tau) & = & \frac{c_{1} \mathcal{S}_{0}}{c_{1} - c_{2}} \left( \exp(-c_{2}\tau) - \exp(-c_{1}\tau) \right),\quad c_{1} \neq c_{2}. \end{eqnarray} | (5.8) |
Fractional-order models have been recognized as being more useful in estimating experimental data than integer-order models. Another reason is that various parameters can be adjusted accordingly. Then, by taking the (k, \psi) –Caputo–PFDO into the left-sided of the problem (5.6), the problem (5.6) can be rewritten as:
\begin{equation} \left\{ \begin{array}{l} {}_{0,k}^{{C}}\mathfrak{D}^{\alpha,\rho;\psi} \mathcal{S}(\tau) = -c_{1} \mathcal{S}(\tau),\qquad\qquad\,\,\alpha\in(0,1],\quad \rho \in (0,1],\quad k > 0, \\ {}_{0,k}^{{C}}\mathfrak{D}^{\alpha,\rho;\psi} \mathcal{B}(\tau) = c_{1} S(\tau) - c_{2} \mathcal{B}(\tau),\quad \alpha\in(0,1],\quad \rho \in (0,1],\quad k > 0,\\ \mathcal{S}(0) = \mathcal{S}_{0},\quad \mathcal{B}(0) = \mathcal{B}_{0} = 0. \end{array} \right. \end{equation} | (5.9) |
Similarly procedure in Example 5.1, we obtain that
\begin{equation} \left\{ \begin{array}{l} A_{N}{{(\psi(\tau)-\psi (0))}^{1-\frac{\alpha }{k}}}\left[ (1-\rho )\mathcal{S}(\tau)+k\rho \frac{\mathcal{S}^{\prime}(\tau)}{\psi^{\prime} (\tau)} \right] \\ \qquad\qquad - e^{\frac{\rho -1}{k\rho }\psi (\tau)}\sum\limits_{i = 1}^{N}{B_{N,i}{{(\psi (\tau)-\psi (0))}^{1-\frac{\alpha }{k}-i}}{{\mathcal{V}}_{i}}(\tau)} = -c_{1} \mathcal{S}(\tau),\\ A_{N}{{(\psi(\tau)-\psi (0))}^{1-\frac{\alpha }{k}}}\left[ (1-\rho )\mathcal{B}(\tau)+k\rho \frac{\mathcal{B}^{\prime}(\tau)}{\psi^{\prime} (\tau)} \right] \\ \qquad\qquad - e^{\frac{\rho -1}{k\rho }\psi (\tau)}\sum\limits_{i = 1}^{N}{B_{N,i}{{(\psi (\tau)-\psi (0))}^{1-\frac{\alpha }{k}-i}}{{\mathcal{U}}_{i}}(\tau)} = c_{1} \mathcal{S}(\tau) - c_{2} \mathcal{B}(\tau), \\ \mathcal{V}^{\prime}_{i}(\tau) = {{\left( \psi (\tau)-\psi (0) \right)}^{i-1}}{\psi }^{\prime}(\tau){{e}^{\frac{1-\rho }{k\rho }\psi (\tau)}}\left[ (1-\rho )\mathcal{S}(\tau)+k\rho \frac{\mathcal{S}^{\prime}(\tau)}{\psi^{\prime}(\tau)} \right], \quad i = 1,2,\ldots,N,\\ \mathcal{U}^{\prime}_{i}(\tau) = {{\left( \psi (\tau)-\psi (0) \right)}^{i-1}}{\psi }^{\prime}(\tau){{e}^{\frac{1-\rho }{k\rho }\psi (\tau)}}\left[ (1-\rho )\mathcal{B}(\tau)+k\rho \frac{\mathcal{B}^{\prime}(\tau)}{\psi^{\prime}(\tau)} \right], \quad i = 1,2,\ldots,N,\\ \mathcal{S}(0) = \mathcal{S}_{0}, \quad \mathcal{B}(0) = \mathcal{B}_{0} = 0, \quad \mathcal{V}_{i}(0) = 0,\quad \mathcal{U}_{i}(0) = 0. \end{array} \right. \end{equation} | (5.10) |
In all graphical simulations, we use the parameters k_{1} = 0.12 , k_{2} = 0.01 , and \psi(\tau) = \tau based on the experimental results for the BALs of a real individual in Table 1. We separated them into three cases, which are shown in Figure 3.
Time (min) | 0 | 10 | 20 | 30 | 45 | 80 | 90 | 110 | 170 |
BAL (mg/L) | 0 | 150 | 200 | 160 | 130 | 170 | 60 | 40 | 20 |
Case (1): We set the initial conditions \mathcal{S}_{0} = 330 , \mathcal{B}_{0} = 0 , \rho = 0.988 , and k = 0.995 with the varied \alpha \in \{0.80, 0.85, 0.90, 0.95, 1.00\} . The graphical simulation of B(\tau) with the common parameters is shown as in Figure 3a. Case (2): We set the initial conditions \mathcal{S}_{0} = 280 , \mathcal{B}_{0} = 0 , \alpha = 0.95 , and k = 0.9 with the varied \rho \in \{0.992, 0.994, 0.996, 0.998, 1.000\} . The graphical simulation of B(\tau) with the common parameters is shown as in Figure 3b. Case (3): We set the initial conditions \mathcal{S}_{0} = 255 , \mathcal{B}_{0} = 0 , \alpha = 0.95 , and \rho = 0.998 with the varied k \in \{0.88, 0.90, 0.92, 0.94, 0.96\} . The graphical simulation of B(\tau) with the common parameters is shown as in Figure 3c. Case (4): We set the initial conditions \mathcal{S}_{0} = 255 , \mathcal{B}_{0} = 0 , \alpha = 0.90 , \rho = 0.25 , k = 0.5 , and N = 2, 5, 20 . The graphical simulation of B(\tau) with the common parameters is shown as in Figure 4.
We calculate the relative absolute error, E_{N} = \vert (x_{e}-x_{a})/x_{e} \vert between the exact solution x_{e}(\tau) and the approximation solution x_{a}(\tau) for the concentration of alcohol in the blood of a human body at time \tau , B(\tau) , as displayed in Table 2. The presented technique is applicable. Also, the results agree with the exact solutions, and the error decreases as N increases.
Time (min) | BAL (mg/L) | \mathcal{S}_2 | \mathcal{S}_5 | \mathcal{S}_{20} | E_2 | E_5 | E_{20} |
0 | 0 | 0.0000 | 0.0000 | 0.0000 | - | - | - |
10 | 150 | 182.4541 | 183.2365 | 183.7853 | 0.2164 | 0.2216 | 0.2252 |
20 | 200 | 198.0376 | 198.9739 | 199.5539 | 0.0098 | 0.0051 | 0.0022 |
30 | 160 | 176.2629 | 176.7983 | 177.0553 | 0.1016 | 0.1050 | 0.1066 |
45 | 130 | 138.3014 | 138.3158 | 138.2298 | 0.0639 | 0.0640 | 0.0633 |
80 | 70 | 75.1240 | 74.6521 | 74.2932 | 0.0732 | 0.0665 | 0.0613 |
90 | 60 | 62.7600 | 62.2160 | 61.8149 | 0.0460 | 0.0369 | 0.0302 |
110 | 40 | 43.4365 | 42.7890 | 42.3281 | 0.0859 | 0.0697 | 0.0582 |
170 | 20 | 13.0767 | 12.3641 | 11.8921 | 0.3462 | 0.3818 | 0.4054 |
In this work, we analyzed an extended Gronwall inequality in the context of the (k, \psi) -PFOs and proved its properties. A Cauchy-type problem under the (k, \psi) -Caputo-PFDO is the highlight of our presentation. First, the existence and uniqueness results of the proposed problem (1.1) were established by applying Banach's contraction mapping principle and Leray-Schauder's nonlinear alternative. Then, UH-ML stability was studied to guarantee the existing results. Moreover, a numerical technique is demonstrated based on a decomposition formula for the (k, \psi) -Caputo-PFDO. Finally, numerical examples are shown to verify the theoretical results. As seen in Example 5.1, we set the specific type of Cauchy problem and showed numeric calculations to confirm the accuracy of Theorem 3.2, which expresses the unique solution. We also showed the calculation satisfied the conditions for various Ulam stability to guarantee the results. For the fixed value of order \alpha and with the different parameters \rho and k , we found the result's behavior through the graphic numerical simulation that the approximate solution of the system converges to the exact solution as \rho and k increased. In contrast, the absolute error values decreased to zero for increasing value k but increased for decreasing value \rho . Furthermore, to strengthen our study in a broader domain, we applied the established numerical technique to the BALs problem as seen in Example 5.2. We expressed a numerical simulation of this problem under different parameters. The graphical results illustrated that various values of \alpha, \rho, and k gave the corresponding behavior of the system in the same trend as the data. Consequently, the fractional-order form may approximate the alcohol concentration in a human's blood under a range of circumstances, as proven with data.
In future works, we can discuss applying (k, \psi) -Caputo-PFDO to real-world problems. This larger paradigm greatly contributes to the development of fractional calculus. It is paving the way for exciting future studies in this dynamic and developing discipline.
Weerawat Sudsutad: Problem statement, conceptualization, methodology, investigation, software, writing the original draft, writing, reviewing, and editing; Jutarat Kongson: Supervision, problem statement, conceptualization, methodology, investigation, writing the original draft, reviewing, and editing; Chatthai Thaiprayoon: Methodology, software, investigation, writing the original draft, reviewing, and editing; Aphirak Aphithana, Weerapan Sae-dan: Investigation, software, writing, reviewing, and editing. All the authors read and approved the final manuscript.
W. Sudsutad (weerawat.s@rumail.ru.ac.th), A. Aphithana (aphirak.apt@gmail.com), W. Sae-dan (weerapan.s@rumail.ru.ac.th) would like to thank Ramkhamhaeng University, J. Kongson and C. Thaiprayoon (chatthai@go.buu.ac.th) thank Burapha University for supporting this work.
The authors declare no conflicts of interest.
[1] |
D. A. Brannan, W. E. Kirwan, On some classes of bounded univalent functions, J. Lond. Math. Soc., s2-1 (1969), 431–443. https://doi.org/10.1112/jlms/s2-1.1.431 doi: 10.1112/jlms/s2-1.1.431
![]() |
[2] | C. Carathéodory, Über den variabilitätsbereich der coeffizienten von potenzreihen, die gegebene werte nicht annehmen, Math. Ann., 64 (1907), 95–115. |
[3] | C. Carathéodory, Über den variabilitätsbereich der fourier'schen konstanten von positiven harmonischen funktionen, Rend. Circ. Mat. Palermo, 32 (1911), 193–217. |
[4] | P. J. Eenigenburg, S. S. Miller, P. T. Mocanu, M. O. Reade, On a Briot-Bouquet differential subordination, General Inequalities 3, International Series of Numerical Mathematics, Basel: Birkh\ddot{a}user, 1983,339–348. |
[5] | I. H. Kim, N. E. Cho, Sufficient conditions for Carathéodory functions, Comput. Math. Appl., 59 (2010), 2067–2073. |
[6] | O. S. Kwon, Y. J. Sim, Sufficient conditions for Carathéodory functions and applications to univalent functions, Math. Slovaca, 69 (2019), 1065–1076. |
[7] | M. S. Liu, Y. C. Zhu, Criteria for strongly starlike and \Phi-like functions, Complex Var. Elliptic Equ., 53 (2008), 485–500. |
[8] |
S. S. Miller, P. T. Mocanu, Univalent solutions of Briot-Bouquet differential equations, J. Differ. Equ., 56 (1985), 297–309. https://doi.org/10.1016/0022-0396(85)90082-8 doi: 10.1016/0022-0396(85)90082-8
![]() |
[9] | S. S. Miller, P. T. Mocanu, Mark-Strohh\ddot{a}cker differential subordination systems, Proc. Amer. Math. Soc., 99 (1987), 527–534. |
[10] | S. S. Miller, P. T. Mocanu, On some classes of first-order differential subordinations, Michigan Math. J., 32 (1985), 185–195. |
[11] | S. S. Miller, P. T. Mocanu, A special differential subordination and its application to univalency conditions, Current Topics in Analytic Function Theory, Singapore, London: World Scientific, 1992,171–185. |
[12] |
S. S. Miller, P. T. Mocanu, Briot-Bouquet differential equations and differential subordinations, Complex Var. Theory Appl.: Int. J., 33 (1997), 217–237. https://doi.org/10.1080/17476939708815024 doi: 10.1080/17476939708815024
![]() |
[13] | S. S. Miller, P. T. Mocanu, Differential Subordinations: Theory and Applications, New York, Basel: Marcel Dekker, 2000. |
[14] | P. T. Mocanu, Alpha-convex integral operators and strongly starlike functions, Studia Univ. Babes-Bolyai Math., 34 (1989), 18–24. |
[15] |
M. Nunokawa, On the order of strongly starlikeness of strongly convex functions, Proc. Japan Acad. Ser. A Math. Sci., 69 (1993), 234–237. https://doi.org/10.3792/pjaa.69.234 doi: 10.3792/pjaa.69.234
![]() |
[16] | M. Nunokawa, O. S. Kwon, Y. J. Sim, N. E. Cho, Sufficient conditions for Carathéodory functions, Filomat, 32 (2018), 1097–1106. |
[17] |
M. S. Robertson, On the theory of univalent functions, Ann. Math., 37 (1936), 374–408. https://doi.org/10.2307/1968451 doi: 10.2307/1968451
![]() |
[18] | Y. J. Sim, O. S. Kwon, N. E. Cho, H. M. Srivastava, Some sets of sufficient conditions for Carathéodory functions, J. Comput. Anal. Appl., 21 (2016), 1243–1254. |
[19] | R. Singh, On Bazilevć functions, Proc. Amer. Math. Soc., 38 (1973), 261–271. |
[20] | J. Stankiewicz, Quelques problèmes extrémaux dans les classes des fonctions \alpha-angularirement étoilées, Ann. Univ. Mariae Curie-Skłodowska Sect. A, 20 (1966), 59–75. |
[21] | J. Stankiewicz, On a family of starlike functions, Ann. Univ. Mariae Curie-Skłodowska Sect. A, 22-24 (1968–1970), 175–181. |
[22] |
L. M. Wang, The tilted Carathéodory class and its applications, J. Korean Math. Soc., 49 (2012), 671–686. http://doi.org/10.4134/JKMS.2012.49.4.671 doi: 10.4134/JKMS.2012.49.4.671
![]() |
[23] |
D. R. Wilken, J. Feng, A remark on convex and starlike functions, J. Lond. Math. Soc., s2-21 (1980), 287–290. https://doi.org/10.1112/jlms/s2-21.2.287 doi: 10.1112/jlms/s2-21.2.287
![]() |
1. | Jutarat Kongson, Chatthai Thaiprayoon, Weerawat Sudsutad, Analysis of a mathematical model for the spreading of the monkeypox virus with constant proportional-Caputo derivative operator, 2025, 10, 2473-6988, 4000, 10.3934/math.2025187 |
Time (min) | 0 | 10 | 20 | 30 | 45 | 80 | 90 | 110 | 170 |
BAL (mg/L) | 0 | 150 | 200 | 160 | 130 | 170 | 60 | 40 | 20 |
Time (min) | BAL (mg/L) | \mathcal{S}_2 | \mathcal{S}_5 | \mathcal{S}_{20} | E_2 | E_5 | E_{20} |
0 | 0 | 0.0000 | 0.0000 | 0.0000 | - | - | - |
10 | 150 | 182.4541 | 183.2365 | 183.7853 | 0.2164 | 0.2216 | 0.2252 |
20 | 200 | 198.0376 | 198.9739 | 199.5539 | 0.0098 | 0.0051 | 0.0022 |
30 | 160 | 176.2629 | 176.7983 | 177.0553 | 0.1016 | 0.1050 | 0.1066 |
45 | 130 | 138.3014 | 138.3158 | 138.2298 | 0.0639 | 0.0640 | 0.0633 |
80 | 70 | 75.1240 | 74.6521 | 74.2932 | 0.0732 | 0.0665 | 0.0613 |
90 | 60 | 62.7600 | 62.2160 | 61.8149 | 0.0460 | 0.0369 | 0.0302 |
110 | 40 | 43.4365 | 42.7890 | 42.3281 | 0.0859 | 0.0697 | 0.0582 |
170 | 20 | 13.0767 | 12.3641 | 11.8921 | 0.3462 | 0.3818 | 0.4054 |
Time (min) | 0 | 10 | 20 | 30 | 45 | 80 | 90 | 110 | 170 |
BAL (mg/L) | 0 | 150 | 200 | 160 | 130 | 170 | 60 | 40 | 20 |
Time (min) | BAL (mg/L) | \mathcal{S}_2 | \mathcal{S}_5 | \mathcal{S}_{20} | E_2 | E_5 | E_{20} |
0 | 0 | 0.0000 | 0.0000 | 0.0000 | - | - | - |
10 | 150 | 182.4541 | 183.2365 | 183.7853 | 0.2164 | 0.2216 | 0.2252 |
20 | 200 | 198.0376 | 198.9739 | 199.5539 | 0.0098 | 0.0051 | 0.0022 |
30 | 160 | 176.2629 | 176.7983 | 177.0553 | 0.1016 | 0.1050 | 0.1066 |
45 | 130 | 138.3014 | 138.3158 | 138.2298 | 0.0639 | 0.0640 | 0.0633 |
80 | 70 | 75.1240 | 74.6521 | 74.2932 | 0.0732 | 0.0665 | 0.0613 |
90 | 60 | 62.7600 | 62.2160 | 61.8149 | 0.0460 | 0.0369 | 0.0302 |
110 | 40 | 43.4365 | 42.7890 | 42.3281 | 0.0859 | 0.0697 | 0.0582 |
170 | 20 | 13.0767 | 12.3641 | 11.8921 | 0.3462 | 0.3818 | 0.4054 |