In this paper, we investigate the existence and uniqueness of solutions to a nonlinear coupled systems of (k,φ)-Hilfer fractional differential equations supplemented with nonlocal integro-multi-point boundary conditions. We make use of the Banach contraction mapping principle to obtain the uniqueness result, while the existence results are proved with the aid of Krasnosel'ski⌣i's fixed point theorem and Leray-Schauder alternative for the given problem. Examples demonstrating the application of the abstract results are also presented. Our results are of quite general nature and specialize in several new results for appropriate values of the parameters β1, β2, and the function φ involved in the problem at hand.
Citation: Nattapong Kamsrisuk, Sotiris K. Ntouyas, Bashir Ahmad, Ayub Samadi, Jessada Tariboon. Existence results for a coupled system of (k,φ)-Hilfer fractional differential equations with nonlocal integro-multi-point boundary conditions[J]. AIMS Mathematics, 2023, 8(2): 4079-4097. doi: 10.3934/math.2023203
[1] | Ibtesam Alshammari, Islam M. Taha . On fuzzy soft β-continuity and β-irresoluteness: some new results. AIMS Mathematics, 2024, 9(5): 11304-11319. doi: 10.3934/math.2024554 |
[2] | R. Mareay, Radwan Abu-Gdairi, M. Badr . Soft rough fuzzy sets based on covering. AIMS Mathematics, 2024, 9(5): 11180-11193. doi: 10.3934/math.2024548 |
[3] | Fahad Alsharari, Ahmed O. M. Abubaker, Islam M. Taha . On r-fuzzy soft γ-open sets and fuzzy soft γ-continuous functions with some applications. AIMS Mathematics, 2025, 10(3): 5285-5306. doi: 10.3934/math.2025244 |
[4] | Abdelghani Taouti, Waheed Ahmad Khan . Fuzzy subnear-semirings and fuzzy soft subnear-semirings. AIMS Mathematics, 2021, 6(3): 2268-2286. doi: 10.3934/math.2021137 |
[5] | Arife Atay . Disjoint union of fuzzy soft topological spaces. AIMS Mathematics, 2023, 8(5): 10547-10557. doi: 10.3934/math.2023535 |
[6] | Samirah Alzahrani, A. A. Nasef, N. Youns, A. I. EL-Maghrabi, M. S. Badr . Soft topological approaches via soft γ-open sets. AIMS Mathematics, 2022, 7(7): 12144-12153. doi: 10.3934/math.2022675 |
[7] | Rui Gao, Jianrong Wu . Filter with its applications in fuzzy soft topological spaces. AIMS Mathematics, 2021, 6(3): 2359-2368. doi: 10.3934/math.2021143 |
[8] | Warud Nakkhasen, Teerapan Jodnok, Ronnason Chinram . Intra-regular semihypergroups characterized by Fermatean fuzzy bi-hyperideals. AIMS Mathematics, 2024, 9(12): 35800-35822. doi: 10.3934/math.20241698 |
[9] | Fenhong Li, Liang Kong, Chao Li . Non-global nonlinear mixed skew Jordan Lie triple derivations on prime ∗-rings. AIMS Mathematics, 2025, 10(4): 7795-7812. doi: 10.3934/math.2025357 |
[10] | Umar Ishtiaq, Fahad Jahangeer, Doha A. Kattan, Manuel De la Sen . Generalized common best proximity point results in fuzzy multiplicative metric spaces. AIMS Mathematics, 2023, 8(11): 25454-25476. doi: 10.3934/math.20231299 |
In this paper, we investigate the existence and uniqueness of solutions to a nonlinear coupled systems of (k,φ)-Hilfer fractional differential equations supplemented with nonlocal integro-multi-point boundary conditions. We make use of the Banach contraction mapping principle to obtain the uniqueness result, while the existence results are proved with the aid of Krasnosel'ski⌣i's fixed point theorem and Leray-Schauder alternative for the given problem. Examples demonstrating the application of the abstract results are also presented. Our results are of quite general nature and specialize in several new results for appropriate values of the parameters β1, β2, and the function φ involved in the problem at hand.
In this paper, we consider the existence of solutions and a generalized Lyapunov-type inequality to the following boundary value problem for differential equation of variable order
{Dq(t)0+x(t)+f(t,x)=0, 0<t<T,x(0)=0,x(T)=0, | (1.1) |
where 0<T<+∞, Dq(t)0+ denotes derivative of variable order([1,2,3,4]) defined by
Dq(t)0+x(t)=d2dt2∫t0(t−s)1−q(s)Γ(2−q(s))x(s)ds,t>0, | (1.2) |
and
I2−q(t)0+x(t)=∫t0(t−s)1−q(s)Γ(2−q(s))x(s)ds,t>0, | (1.3) |
denotes integral of variable order 2−q(t), 1<q(t)≤2, 0≤t≤T. f:(0,T]×R→R is given continuous function satisfying some assumption conditions.
The operators of variable order, which fall into a more complex operator category, are the derivatives and integrals whose order is the function of certain variables. The variable order fractional derivative is an extension of constant order fractional derivative. In recent years, the operator and differential equations of variable order have been applied in engineering more and more frequently, for the examples and details, see [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17].
The subject of fractional calculus has gained considerable popularity and importance due to its frequent appearance in different research areas and engineering, such as physics, chemistry, control of dynamical systems etc. Recently, many people paid attention to the existence and uniqueness of solutions to boundary value problems for fractional differential equations. Although the existing literature on solutions of boundary value problems of fractional order (constant order) is quite wide, few papers deal with the existence of solutions to boundary value problems of variable order. According to (1.2) and (1.3), it is obviously that when q(t) is a constant function, i.e. q(t)≡q (q is a finite positive constant), then Iq(t)0+,Dq(t)0+ are the usual Riemann-Liouville fractional integral and derivative [18].
The following properties of fractional calculus operators Dq0+, Iq0+ play an important part in discussing the existence of solutions of fractional differential equations.
Proposition 1.1. [18] The equality Iγ0+Iδ0+f(t)=Iγ+δ0+f(t), γ>0,δ>0 holds for f∈L(0,b),0<b<+∞.
Proposition 1.2. [18] The equality Dγ0+Iγ0+f(t)=f(t), γ>0 holds for f∈L(0,b),0<b<+∞.
Proposition 1.3. [18] Let 1<α≤2. Then the differential equation
Dα0+f=0 |
has solutions
f(t)=c1tα−1+c2tα−2,c1,c2∈R. |
Proposition 1.4. [18] Let 1<α≤2, f(t)∈L(0,b), Dα0+f∈L(0,b). Then the following equality holds
Iα0+Dα0+f(t)=f(t)+c1tα−1+c2tα−2,c1,c2∈R. |
These properties play a very important role in considering the existence of the solutions of differential equations for the Riemnn-Liouville fractional derivative, for details, please refer to [18]. However, from [1,2,16], for general functions h(t),g(t), we notice that the semigroup property doesn't hold, i.e., Ih(t)a+Ig(t)a+≠Ih(t)+g(t)a+. Thus, it brings us extreme difficulties, we can't get these properties like Propositions 1.1-1.4 for the variable order fractional operators (integral and derivative). Without these properties for variable order fractional derivative and integral, we can hardly consider the existence of solutions of differential equations for variable order derivative by means of nonlinear functional analysis (for instance, some fixed point theorems).
Let's take Proposition 1.1 for example. To begin with the simplest case,
Example 1.5. Let p(t)=t, q(t)=1, f(t)=1,0≤t≤3. Now, we calculate Ip(t)0+Iq(t)0+f(t)|t=1 and Ip(t)+q(t)0+f(t)|t=1 defined in (1.3).
Ip(t)0+Iq(t)0+f(t)|t=1=∫10(1−s)s−1Γ(s)∫s0(s−τ)1−1Γ(1)dτds=∫10(1−s)s−1sΓ(s)ds≈0.472. |
and
Ip(t)+q(t)0+f(t)|t=1=∫10(1−s)sΓ(s+1)ds=∫10(1−s)ssΓ(s)ds≈0.686. |
Therefore,
Ip(t)0+Iq(t)0+f(t)|t=1≠Ip(t)+q(t)0+f(t)|t=1. |
As a result, the Propositions 1.2 and Propositions 1.4 do not hold for Dp(t)0+ and Ip(t)0+, such as, for function f∈L(0,T),0<p(t)<1,0≤t≤T, we get
Dp(t)0+Ip(t)0+f(t)=D1(I1−p(t)0+Ip(t)0+f(t))≠D1I1−p(t)+p(t)0+f(t)=f(t),t∈(0,T], |
since we know that I1−p(t)0+Ip(t)0+f(t)≠I1−p(t)+p(t)0+f(t) for general function f.
Now, we can conclude that Propositions 1.1–1.4 do not hold for Dq(t)0+ and Iq(t)0+.
So, one can not transform a differential equation of variable order into an equivalent interval equation without the Propositions 1.1–1.4. It is a difficulty for us in dealing with the boundary value problems of differential equations of variable order. Since the equations described by the variable order derivatives are highly complex, difficult to handle analytically, it is necessary and significant to investigate their solutions.
In [16], by means of Banach Contraction Principle, Zhang considered the uniqueness result of solutions to initial value problem of differential equation of variable order
{Dp(t)0+x(t)=f(t,x),0<t≤T,x(0)=0, | (1.4) |
where 0<T<+∞, Dp(t)0+ denotes derivative of variable order p(t) ([1,2,3,4]) defined by
Dp(t)0+x(t)=ddt∫t0(t−s)−p(t)Γ(1−p(t))x(s)ds,t>0. | (1.5) |
and 1Γ(1−p(t))∫t0(t−s)−p(t)x(s)ds is integral of variable order 1−p(t) for function x(t). And p:[0,T]→(0,1] is a piecewise constant function with partition P={[0,T1],(T1,T2],(T2,T3],⋯, (TN∗−1,T]} (N∗ is a given natural number) of the finite interval [0,T], i.e.
p(t)=N∗∑k=1qkIk(t),t∈[0,T], |
where 0<qk≤1,k=1,2,⋯,N∗ are constants, and Ik is the indicator of the interval [Tk−1,Tk],k=1,2,⋯,N∗(here T0=0,TN∗=T), that is Ik=1 for t∈[Tk−1,Tk], Ik=0 for elsewhere.
In [17], the authors studied the Cauchy problem for variable order differential equations with a piecewise constant order function[19]. Inspired by these works, we will study the boundary value problem (1.1) for variable order differential equation with a piecewise constant order function q(t) in this paper.
Lyapunov's inequality is an outstanding result in mathematics with many different applications, see [20,21,22,23,24,25] and references therein. The result, as proved by Lyapunov[20] in 1907, asserts that if h:[a,b]→R is a continuous function, then a necessary condition for the boundary value problem
{y″(t)+h(t)y(t)=0,a<t<b,y(a)=y(b)=0, | (1.6) |
to have a nontrivial solution is given by
∫ba|h(s)|ds>4b−a, | (1.7) |
where −∞<a<b<+∞.
Lyapunov's inequality has taken many forms, including versions in the context of fractional (noninteger order) calculus, where the second-order derivative in (1.6) is substituted by a fractional operator of order α,
{Dαa+y(t)+h(t)y(t)=0,a<t<b,y(a)=y(b)=0, | (1.8) |
where Dαa+ is the Riemann-Liouville derivative of order α∈(1,2] and h:[a,b]→R is a continuous function. If (1.8) has a nontrivial solution, then
∫ba|h(s)|ds>Γ(α)(4b−a)α−1. |
A Lyapunov fractional inequality can also be obtained by considering the fractional derivative in in the sense of Caputo instead of Riemann-Liouville [22]. More recently, there are some results of Lyapunov type inequalities for fractional boundary value problems. see [23,24]. In [25], authors obtained a generalization of inequality to boundary value problem as following
{Dαa+y(t)+h(t)f(y)=0,a<t<b,y(a)=y(b)=0, | (1.9) |
where Dαa+ is the Riemann-Liouville derivative, 1<α≤2, and h:[a,b]→R is a Lebesgue integrable function. Under some assumptions on the nonlinear term f, authors got a generalization of inequality to the boundary value problem (1.9).
∫ba|h(s)|ds>4α−1Γ(α)η(b−a)α−1f(η), | (1.10) |
where η is maximum value of nontrivial solution to the boundary value problem (1.9).
Motivated by [21,22,23,24,25] and the above results, we focus on a generalized Lyapunov-type inequality to the boundary value problem (1.1) under certain assumptions of nonlinear term.
The paper is organized as following. In Section 2, we provide some necessary definitions associated with the boundary value problem (1.1). In Section 3, we establish the existence of solutions for the boundary value problem (1.1) by using the Schauder fixed point theorem. In Section 4, we investigative the generalized Lyapunov-type inequalities to the boundary value problem (1.1). In section 5, we give some examples are presented to illustrate the main results.
For the convenience of the reader, we present here some necessary definitions that will be used to prove our main results.
Definition 2.1. A generalized interval is a subset I of R which is either an interval (i.e. a set of the form [a,b],(a,b),[a,b) or (a,b]); a point {a}; or the empty set ∅.
Definition 2.2. If I is a generalized interval. A partition of I is a finite set P of generalized intervals contained in I, such that every x in I lies in exactly one of the generalized intervals J in P.
Example 2.3. The set P={{1},(1,6),[6,7),{7},(7,8]} of generalized intervals is a partition of [1,8].
Definition 2.4. Let I be a generalized interval, let f:I→R be a function, and let P a partition of I. f is said to be piecewise constant with respect to P if for every J∈P, f is constant on J.
Example 2.5. The function f:[1,6]→R defined by
f(x)={3, 1≤x<3,4, x=3,5, 3<x<6,2, x=6, |
is piecewise constant with respect to the partition {[1,3],{3},(3,6),{6}} of [1,6].
The following example illustrates that the semigroup property of the variable order fractional integral doesn't holds for the piecewise constant functions p(t) and q(t) defined in the same partition of finite interval [a,b].
Example 2.6. Let p(t)={3, 0≤t≤1,2, 1<t≤4, q(t)={2, 0≤t≤1,3, 1<t≤4, and f(t)=1,0≤t≤4. We'll verify Ip(t)0+Iq(t)0+f(t)|t=2≠Ip(t)+q(t)0+f(t)|t=2, here, the variable order fractional integral is defined in (1.3). For 1≤t≤4, we have
Ip(t)0+Iq(t)0+f(t)=∫10(t−s)p(s)−1Γ(p(s))∫s0(s−τ)q(τ)−1Γ(q(τ))dτds+∫t1(t−s)p(s)−1Γ(p(s))∫s0(s−τ)q(τ)−1Γ(q(τ))dτds=∫10(t−s)2Γ(3)∫s0(s−τ)Γ(2)dτds+∫t1(t−s)Γ(2)[∫10(s−τ)Γ(2)dτ+∫s1(s−τ)2Γ(3)dτ]ds=∫10(t−s)2s22Γ(3)ds+∫t1(t−s)[s22−(s−1)22+(s−1)36]ds=∫10(t−s)2s22Γ(3)ds+16∫t1(t−s)(s3−3s2+9s−4)ds, |
thus, we have
Ip(t)0+Iq(t)0+f(t)|t=2=∫10(2−s)2s22Γ(3)ds+16∫21(2−s)(s3−3s2+9s−4)ds=215+1760=512 |
Ip(t)+q(t)0+f(t)|t=2=∫20(2−s)p(s)+q(s)−1Γ(p(s)+q(s))ds=∫10(2−s)3+2−1Γ(3+2)ds+∫21(2−s)2+3−1Γ(2+3)ds=31120+1120=415. |
Therefore, we obtain
Ip(t)0+Iq(t)0+f(t)|t=2≠Ip(t)+q(t)0+f(t)|t=2, |
which implies that the semigroup property of the variable order fractional integral doesn't hold for the piecewise constant functions p(t) and q(t) defined in the same partition [0,1],(1,4] of finite interval [0,4].
We need the following assumptions:
(H1) Let n∗∈N be an integer, P={[0,T1],(T1,T2],(T2,T3],⋯,(Tn∗−1,T]} be a partition of the interval [0,T], and let q(t):[0,T]→(1,2] be a piecewise constant function with respect to P, i.e.,
q(t)=N∑k=1qkIk(t)={q1, 0≤t≤T1,q2, T1<t≤T2,⋯, ⋯,qn∗, Tn∗−1<t≤Tn∗=T, | (3.1) |
where 1<qk≤2(k=1,2,⋯,n∗) are constants, and Ik is the indicator of the interval [Tk−1,Tk], k=1,2,⋯,n∗ (here T0=0,Tn∗=T), that is, Ik(t)=1 for t∈[Tk−1,Tk] and Ik(t)=0 for elsewhere.
(H2) Let trf:[0,T]×R→R be a continuous function (0≤r<1), there exist constants c1>0,c2>0, 0<γ<1 such that
tr|f(t,x(t))|≤c1+c2|x(t)|γ,0≤t≤T,x(t)∈R. |
In order to obtain our main results, we firstly carry on essential analysis to the boundary value problem (1.1).
By (1.2), the equation of the boundary value problem (1.1) can be written as
d2dt2∫t0(t−s)1−q(s)Γ(2−q(s))x(s)ds+f(t,x)=0, 0<t<T, | (3.2) |
According to (H1), Eq (3.2) in the interval (0,T1] can be written as
Dq10+x(t)+f(t,x)=0, 0<t≤T1. | (3.3) |
Equation (3.2) in the interval (T1,T2] can be written by
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x(s)ds+∫tT1(t−s)1−q2Γ(2−q2)x(s)ds)+f(t,x)=0, | (3.4) |
and Eq (3.2) in the interval (T2,T3] can be written by
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x(s)ds+∫T2T1(t−s)1−q2Γ(2−q2)x(s)ds+∫tT2(t−s)1−q3Γ(2−q3)x(s)ds)+f(t,x)=0. | (3.5) |
In the same way, Eq (3.2) in the interval (Ti−1,Ti],i=4,5,⋯,n∗−1 can be written by
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x(s)ds+⋯+∫tTi−1(t−s)1−qiΓ(2−qi)x(s)ds)+f(t,x)=0. | (3.6) |
As for the last interval (Tn∗−1,T), similar to above argument, Eq (3.2) can be written by
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x(s)ds+⋯+∫tTn∗−1(t−s)1−qn∗Γ(2−qn∗)x(s)ds)+f(t,x)=0. | (3.7) |
Remark 3.1. From the arguments above, we find that, according to condition (H1), in the different interval, the equation of the boundary value problem (1.1) must be represented by different expression. For instance, in the interval (0,T1], the equation of the boundary value problem (1.1) is represented by (3.3); in the interval (T1,T2], the equation of the boundary value problem (1.1) is represented by (3.4); in the interval (T2,T3], the equation of the boundary value problem (1.1) is represented by (3.5), etc. But, as far as we know, in the different intervals, the equation of integer order or constant fractional order problems may be represented by the same expression. Based these facts, different than integer order or constant fractional order problems, in order to consider the existence results of solution to the boundary value problem (1.1), we need consider the relevant problem defined in the different interval, respectively.
Now, based on arguments previous, we present definition of solution to the boundary value problem (1.1), which is fundamental in our work.
Definition 3.2. We say the boundary value problem (1.1) has a solution, if there exist functions xi(t),i=1,2,⋯,n∗ such that x1∈C[0,T1] satisfying equation (3.3) and x1(0)=0=x1(T1); x2∈C[0,T2] satisfying equation (3.4) and x2(0)=0=x2(T2); x3∈C[0,T3] satisfying equation (3.5) and x3(0)=0=x3(T3); xi∈C[0,Ti] satisfying equation (3.6) and xi(0)=0=xi(Ti)(i=4,5,⋯,n∗−1); xn∗∈C[0,T] satisfying equation (3.7) and xn∗(0)=xn∗(T)=0.
Theorem 3.3. Assume that conditions (H1) and (H2) hold, then the boundary value problem (1.1) has one solution.
Proof. According the above analysis, the equation of the boundary value problem (1.1) can be written as Eq (3.2). Equation (3.2) in the interval (0,T1] can be written as
Dq10+x(t)+f(t,x)=0, 0<t≤T1. |
Now, we consider the following two-point boundary value problem
{Dq10+x(t)+f(t,x)=0, 0<t<T1,x(0)=0,x(T1)=0. | (3.8) |
Let x∈C[0,T1] be solution of the boundary value problem (3.8). Now, applying the operator Iq10+ to both sides of the above equation. By Propositions 1.4, we have
x(t)=d1tq1−1+d2tq1−2−1Γ(q1)∫t0(t−s)q1−1f(s,x(s))ds,0<t≤T1. |
By x(0)=0 and the assumption of function f, we could get d2=0. Let x(t) satisfying x(T1)=0, thus we can get d1=Iq10+f(T1,x)T1−q11. Then, we have
x(t)=Iq10+f(T1,x)T1−q11tq1−1−Iq10+f(t,x),0≤t≤T1 | (3.9) |
Conversely, let x∈C[0,T1] be solution of integral Eq (3.9), then, by the continuity of function trf and Proposition 1.2, we can easily get that x is the solution of boundary value problem (3.8).
Define operator T:C[0,T1]→C[0,T1] by
Tx(t)=Iq10+f(T1,x)T1−q11tq1−1−Iq10+f(t,x(t)),0≤t≤T1. |
It follows from the properties of fractional integrals and assumptions on function f that the operator T:C[0,T1]→C[0,T1] defined above is well defined. By the standard arguments, we could verify that T:C[0,T1]→C[0,T1] is a completely continuous operator.
In the next analysis, we take
M(r,q)=max{2(1−r)Γ(q1),2(1−r)Γ(q2),⋯,2(1−r)Γ(qn∗)}. |
Let Ω={x∈C[0,T1]:‖x‖≤R} be a bounded closed convex subset of C[0,T1], where
R=max{2c1M(r,q)(T+1)2,(2c2M(r,q)(T+1)2)11−γ}. |
For x∈Ω and by (H2), we have
|Tx(t)|≤T1−q11tq1−1Γ(q1)∫T10(T1−s)q1−1|f(s,x(s))|ds+1Γ(q1)∫t0(t−s)q1−1|f(s,x(s))|ds≤2Γ(q1)∫T10(T1−s)q1−1|f(s,x(s))|ds≤2Γ(q1)∫T10(T1−s)q1−1s−r(c1+c2|x(s)|γ)ds≤2Tq1−11Γ(q1)∫T10s−r(c1+c2Rγ)ds≤2Tq1−11T1−r1(1−r)Γ(q1)(c1+c2Rγ)≤M(r,q)Tq1−r1(c1+c2Rγ)≤M(r,q)(T+1)2(c1+c2RRγ−1)≤R2+R2=R, |
which means that T(Ω)⊆Ω. Then the Schauder fixed point theorem assures that the operator T has one fixed point x1∈Ω, which is a solution of the boundary value problem (3.8).
Also, we have obtained that Eq (3.2) in the interval (T1,T2] can be written by (3.4). In order to consider the existence result of solution to (3.4), we rewrite (3.4) as following
d2dt2∫T10(t−s)1−q2Γ(2−q2)x(s)ds+d2dt2∫tT1(t−s)1−q2Γ(2−q2)x(s)ds=f(t,x). T1<t≤T2, |
For 0≤s≤T1, we take x(s)≡0, then, by the above equation, we get
Dq2T1x(t)+f(t,x)=0,T1<t<T2. |
Now, we consider the following boundary value problem
{Dq2T1+x(t)+f(t,x)=0, T1<t<T2,x(T1)=0,x(T2)=0, | (3.10) |
Let x∈C[T1,T2] be solution of the boundary value problem (3.10). Now, applying operator Iq2T1+ on both sides of equation to boundary value problem (3.10) and by Propositions 1.4, we have
x(t)=d1(t−T1)q2−1+d2(t−T1)q2−2−1Γ(q2)∫tT1(t−s)q2−1f(s,x(s))ds,T1<t≤T2. |
By x(T1)=0,x(T2)=0, we have d2=0 and d1=Iq2T1+f(T2,x)(T2−T1)1−q2. Then, we have
x(t)=Iq20+f(T2,x)(T2−T1)1−q2(t−T1)q2−1−1Γ(q2)∫tT1(t−s)q2−1f(s,x(s))ds,T1≤t≤T2. |
Conversely, let x∈C[T1,T2] be solution of integral equation above, then, by the continuity assumption of function trf and Proposition (1.2), we can get that x is solution solution of the boundary value problem (3.10).
Define operator T:C[T1,T2]→C[T1,T2] by
Tx(t)=Iq20+f(T2,x)(T2−T1)1−q2(t−T1)q2−1−1Γ(q2)∫tT1(t−s)q2−1f(s,x(s))ds. |
It follows from the continuity of function trf that operator T:C[T1,T2]→C[T1,T2] is well defined. By the standard arguments, we know that T:C[T1,T2]→C[T1,T2] is a completely continuous operator.
For x∈Ω and by (H2), we get
|Tx(t)|≤(T2−T1)1−q2(t−T1)q2−1Γ(q2)∫T2T1(T2−s)q2−1|f(s,x(s))|ds+1Γ(q2)∫tT1(t−s)q2−1|f(s,x(s))|ds≤2Γ(q2)∫T2T1(T2−s)q2−1|f(s,x(s))|ds≤2Γ(q2)∫T2T1(T2−s)q2−1s−r(c1+c2|x(s)|γ)ds≤2Tq2−12Γ(q2)∫T2T1s−r(c1+c2Rγ)ds=2Tq2−12(T1−r2−T1−r1)(1−r)Γ(q2)(c1+c2Rγ)≤2Tq2−r2(1−r)Γ(q2)(c1+c2Rγ)≤M(r,q)(T+1)2(c1+c2RRγ−1)≤R2+R2=R, |
which means that T(Ω)⊆Ω. Then the Schauder fixed point theorem assures that operator T has one fixed point ˜x2∈Ω, which is one solution of the following integral equation, that is,
˜x2(t)=Iq20+f(T2,˜x2)(T2−T1)1−q2(t−T1)q2−1−1Γ(q2)∫tT1(t−s)q2−1f(s,˜x2(s))ds,T1≤t≤T2. | (3.11) |
Applying operator Dq2T1+ on both sides of (3.11), by Proposition 1.2, we can obtain that
Dq2T1+˜x2(t)+f(t,˜x2)=0,T1<t≤T2, |
that is, ˜x2(t) satisfies the following equation
d2dt21Γ(2−q2)∫tT1(t−s)1−q2˜x2(s)ds+f(t,˜x2)=0,T1<t≤T2. | (3.12) |
We let
x2(t)={0, 0≤t≤T1,˜x2(t),T1<t≤T2 | (3.13) |
hence, from (3.12), we know that x2∈C[0,T2] defined by (3.13) satisfies equation
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x2(s)ds+∫tT1(t−s)1−q2Γ(2−q2)x2(s)ds)+f(t,x2)=0, |
which means that x2∈C[0,T2] is one solution of (3.4) with x2(0)=0,x2(T2)=˜x2(T2)=0.
Again, we have known that Eq (3.2) in the interval (T2,T3] can be written by (3.5). In order to consider the existence result of solution to Eq (3.5), for 0≤s≤T2, we take x(s)≡0, then, by (3.5), we get
Dq3T2x(t)+f(t,x)=0,T2<t<T3. |
Now, we consider the following boundary value problem
{Dq3T2+x(t)+f(t,x)=0, T2<t<T3,x(T2)=0,x(T3)=0. | (3.14) |
By the standard way, we know that the boundary value problem (3.14) exists one solution ˜x3∈Ω. Since ˜x3 satisfies equation
Dq3T2+˜x3(t)+f(t,˜x3)=0,T2<t≤T3, |
that is, ˜x3(t) satisfies the following equation
d2dt21Γ(2−q3)∫tT2(t−s)1−q3˜x3(s)ds+f(t,˜x3)=0,T2<t≤T3. | (3.15) |
We let
x3(t)={0, 0≤t≤T2,˜x3(t),T2<t≤T3, | (3.16) |
hence, from (3.15), we know that x3∈C[0,T3] defined by (3.16) satisfies equation
d2dt2(∫T10(t−s)1−q1Γ(2−q1)x3(s)ds+∫T2T1(t−s)1−q2Γ(2−q2)x3(s)ds+∫tT2(t−s)1−q3Γ(2−q3)x3(s)ds)+f(t,x3)=0, |
which means that x3∈C[0,T3] is one solution of (3.5) with x3(0)=0,x(T3)=˜x3(T3)=0.
By the similar way, in order to consider the existence of solution to Eq (3.6) defined on [Ti−1,Ti] of (3.2), we can investigate the following two-point boundary value problem
{DqiTi−1+x(t)+f(t,x)=0, Ti−1<t<Ti,x(Ti−1)=0,x(Ti)=0. | (3.17) |
By the same arguments previous, we obtain that the Eq (3.6) defined on [Ti−1,Ti] of (3.2) has solution
xi(t)={0, 0≤t≤Ti−1,˜xi(t),Ti−1<t≤Ti, | (3.18) |
where ˜xi∈Ω with ˜xi(Ti−1)=0=˜xi(Ti), i=4,5,⋯,n∗−1.
Similar to the above argument, in order to consider the existence result of solution to Eq (3.7), we may consider the following boundary value problem
{Dqn∗Tn∗−1+x(t)+f(t,x)=0, Tn∗−1<t<Tn∗=T,x(Tn∗−1)=0,x(T)=0. | (3.19) |
So by the same considering, for Tn∗−1≤t≤T we get
x(t)=(T−Tn∗−1)1−qn∗(t−Tn∗−1)qn∗−1Iqn∗Tn∗−1+f(T,x)−Iqn∗Tn∗−1+f(t,x). |
Define operator T:C[Tn∗−1,T]→C[Tn∗−1,T] by
Tx(t)=(T−Tn∗−1)1−qn∗(t−Tn∗−1)qn∗−1Iqn∗Tn∗−1+f(T,x)−1Γ(qn∗)∫tTn∗−1(t−s)qn∗−1f(s,x(s))ds, |
Tn∗−1≤t≤T. It follows from the continuity assumption of function trf that operator T:C[Tn∗−1,T]→C[Tn∗−1,T] is well defined. By the standard arguments, we note that T:C[Tn∗−1,T]→C[Tn∗−1,T] is a completely continuous operator.
For x∈Ω and by (H2), we get
|Tx(t)|≤(T−Tn∗−1)1−qn∗(t−Tn∗−1)qn∗−1Γ(qn∗)∫TTn∗−1(T−s)qn∗−1|f(s,x(s))|ds+1Γ(qn∗)∫tTn∗−1(t−s)qn∗−1|f(s,x(s))|ds≤2Γ(qn∗)∫TTn∗−1(T−s)qn∗−1|f(s,x(s))|ds≤2Γ(qn∗)∫TTn∗−1(T−s)qn∗−1s−r(c1+c2|x(s)|γ)ds≤2Tqn∗−1Γ(qn∗)∫TTn∗−1s−r(c1+c2Rγ)ds≤2Tqn∗−1(T1−r−T1−rn∗−1)(1−r)Γ(qn∗)(c1+c2Rγ)≤2(T+1)2(1−r)Γ(qn∗)(c1+c2Rγ)≤M(r,q)(T+1)2(c1+c2RRγ−1)≤R2+R2=R, |
which means that T(Ω)⊆Ω. Then the Schauder fixed point theorem assures that operator T has one fixed point ˜xn∗∈Ω, which is one solution of the following integral equation, that is,
˜xn∗(t)=(T−Tn∗−1)1−qn∗(t−Tn∗−1)qn∗−1Iqn∗Tn∗−1+f(T,˜xn∗)−1Γ(qn∗)∫tTn∗−1(t−s)qn∗−1f(s,˜xn∗(s))ds,Tn∗−1≤t≤T. | (3.20) |
Applying operator Dqn∗Tn∗−1+ on both sides of (3.20), by Proposition 1.2, we can obtain that
Dqn∗Tn∗−1+˜xn∗(t)+f(t,˜xn∗)=0,Tn∗−1<t≤T, |
that is, ˜xT(t) satisfies the following equation
d2dt21Γ(2−qn∗)∫tTn∗−1(t−s)1−qn∗˜xn∗(s)ds+f(t,˜xn∗)=0,Tn∗−1<t≤T. | (3.21) |
We let
xn∗(t)={0, 0≤t≤Tn∗−1,˜xn∗(t),Tn∗−1<t≤T, | (3.22) |
hence, from (3.21), we know that xn∗∈C[0,T] defined by (3.22) satisfies equation
d2dt2(∫T10(t−s)1−q1Γ(2−q1)xn∗(s)ds+⋯+∫tTn∗−1(t−s)1−qn∗Γ(2−qn∗)xn∗(s)ds)+f(t,xn∗)=0. |
for Tn∗−1<t<T, which means that xn∗∈C[0,T] is one solution of (3.7) with xn∗(0)=0,xn∗(T)=˜xn∗(T)=0.
As a result, we know that the boundary value problem (1.1) has a solution. Thus we complete the proof.
Remark 3.4. For condition (H2), if γ≥1, then using similar way, we can obtain the existence result of solution to the boundary value problem (1.1) provided that we impose some additional conditions on c1,c2.
In this section, we investigate the generalized Lyapunov-type inequalities for the boundary value problem (1.1).
Now, we explore characters of Green functions to the boundary value problems (3.8), (3.10), (3.14), ⋯, (3.17) and (3.19).
Proposition 4.1. Assume that trf:[0,T]×R→R, (0≤r<1) is continuous function, q(t):[0,T]→(1,2] satisfies (H1), then the Green functions
Gi(t,s)={1Γ(qi)[(Ti−Ti−1)1−qi(t−Ti−1)qi−1(Ti−s)qi−1−(t−s)qi−1],Ti−1≤s≤t≤Ti,1Γ(qi)(Ti−Ti−1)1−qi(t−Ti−1)qi−1(Ti−s)qi−1,Ti−1≤t≤s≤Ti, | (4.1) |
of the boundary value problems (3.8), (3.10), (3.14), ⋯, (3.17) and (3.19) satisfy the following properties:
(1) Gi(t,s)≥0 for all Ti−1≤t,s≤Ti;
(2) maxt∈[Ti−1,Ti]Gi(t,s)=Gi(s,s), s∈[Ti−1,Ti];
(3) Gi(s,s) has one unique maximum given by
maxs∈[Ti−1,Ti]Gi(s,s)=1Γ(qi)(Ti−Ti−14)qi−1, |
where i=1,2,⋯,n∗, T0=0,Tn∗=T.
Proof. From the proof of Theorem 3.1, we know that Green functions of the boundary value problems (3.8), (3.10), (3.14), ⋯, (3.17) and (3.19) are given by (4.1).
Using a similar way, we can verify these three results. In fact, let
g(t,s)=(Ti−Ti−1)1−qi(t−Ti−1)qi−1(Ti−s)qi−1−(t−s)qi−1,Ti−1≤s≤t≤Ti. |
We see that
gt(t,s)=(qi−1)[(Ti−Ti−1)1−qi(t−Ti−1)qi−2(Ti−s)qi−1−(t−s)qi−2]≤(qi−1)[(Ti−Ti−1)1−qi(t−s)qi−2(Ti−Ti−1)qi−1−(t−s)qi−2]=0, |
which means that g(t,s) is nonincreasing with respect to t, so g(t,s)≥g(Ti,s)=0 for Ti−1≤s≤t≤Ti. Thus, together this with the expression of Gi(t,s), we get that Gi(t,s)≥0 for all Ti−1≤t,s≤Ti, i=1,2,⋯,n∗, T0=0,Tn∗=T.
Since g(t,s) is nonincreasing with respect to t, it holds that g(t,s)≤g(s,s) for Ti−1≤s≤t≤Ti. On the other hand, for Ti−1≤t≤s≤Ti, we have
(Ti−Ti−1)1−qi(t−Ti−1)qi−1(Ti−s)qi−1≤(Ti−Ti−1)1−qi(s−Ti−1)qi−1(Ti−s)qi−1. |
These assure that maxt∈[Ti−1,Ti]Gi(t,s)=Gi(s,s), s∈[Ti−1,Ti], i=1,2,⋯,n∗, T0=0,Tn∗=T.
Next, we verify (3) of Proposition 4.1. Obviously, the maximum points of Gi(s,s) are not Ti−1 and Ti, i=1,2,⋯,n∗, T0=0,Tn∗=T. For s∈(Ti−1,Ti), i=1,2,⋯,n∗, T0=0,Tn∗=T, we have that
dGi(s,s)ds=1Γ(qi)(Ti−Ti−1)1−qi(qi−1)[(s−Ti−1)qi−2(Ti−s)qi−1−(s−Ti−1)qi−1(Ti−s)qi−2]=1Γ(qi)(Ti−Ti−1)1−qi(qi−1)(s−Ti−1)qi−2(Ti−s)qi−2[Ti+Ti−1−2s], |
which implies that the maximum points of Gi(s,s) is s=Ti−1+Ti2, i=1,2,⋯,n∗, T0=0,Tn∗=T. Hence, for i=1,2,⋯,n∗,T0=0,Tn∗=T,
maxs∈[Ti−1,Ti]Gi(s,s)=G(Ti−1+Ti2,Ti−1+Ti2)=1Γ(qi)(Ti−Ti−14)qi−1. |
Thus, we complete this proof.
Theorem 4.2. Let (H1) holds and trf:[0,T]×R→R, (0≤r<1) be a continuous function. Assume that there exists nonnegative continuous function h(t) defined on [0,T] such that
tr|f(t,x)|≤h(t)|x(t)|,0≤t≤T,x(t)∈R |
If the boundary value problem (1.1) has a nontrivial solution x, then
∫T0s−rh(s)ds>n∗∑i=1Γ(qi)(4Ti−Ti−1)qi−1. | (4.2) |
Proof. Let x be a nontrivial solution of the boundary value problem (1.1). Using Definition 3.2 and the proof of Theorem 3.3, we know that
x(t)={x1(t),0≤t≤T1x2(t)={0,0≤t≤T1,˜x1(t),T1<t≤T2,x3(t)={0,0≤t≤T2,˜x2(t),T2<t≤T3,⋮xi(t)={0,0≤t≤Ti−1,˜xi(t),Ti−1<t≤Ti,⋮xn∗(t)={0,0≤t≤Tn∗−1,˜xn∗(t),Tn∗−1<t≤T, | (4.3) |
where x1∈C[0,T1] is nontrivial solution of the boundary value problem (3.8) with a1=0, ˜x2∈C[T1,T2] is nontrivial solution of the boundary value problem (3.10) with a2=0, ˜x3∈C[T2,T3] is nontrivial solution of the boundary value problem (3.14) with a3=0, ˜xi∈C[Ti−1,Ti] is nontrivial solution of the boundary value problem (3.17) with ai=0, ˜xn∗∈C[Tn∗−1,T] is nontrivial solution of the boundary value problem (3.19). From (4.3) and Proposition 4.1, we have
‖x1‖C[0,T1]=max0≤t≤T1|x1(t)|≤max0≤t≤T1∫T10G1(t,s)|f(s,x1(s))|ds≤∫T10G1(s,s)s−rh(s)|x1(s)|ds<‖x1‖C[0,T1]Γ(q1)(T14)q1−1∫T10s−rh(s)ds, |
which implies that
∫T10s−rh(s)ds>Γ(q1)(4T1)q1−1. | (4.4) |
‖x2‖C[0,T2]=maxT1≤t≤T2|˜x2(t)|=maxT1≤t≤T2|∫T2T1G2(t,s)f(s,˜x2(s))ds|≤∫T2T1G2(s,s)s−rh(s)|˜x2(s)|ds<‖˜x2‖C[T1,T2]Γ(q2)(T2−T14)q2−1∫T2T1s−rh(s)ds,=‖x2‖C[0,T2]Γ(q2)(T2−T14)q2−1∫T2T1s−rh(s)ds, |
which implies that
∫T2T1s−rh(s)ds>Γ(q2)(4T2−T1)q2−1. | (4.5) |
Similar, for i=3,4,⋯,n∗ (Tn∗=T), we have
‖xi‖C[0,Ti]=maxTi−1≤t≤Ti|˜xi(t)|=maxTi−1≤t≤Ti|∫TiTi−1Gi(t,s)f(s,˜xi(s))ds|≤∫TiTi−1Gi(s,s)s−rh(s)|˜xi(s)|ds<‖˜xi‖C[Ti−1,Ti]Γ(qi)(Ti−Ti−14)qi−1∫TiTi−1s−rh(s)ds,=‖xi‖C[0,Ti]Γ(qi)(Ti−Ti−14)qi−1∫TiTi−1s−rh(s)ds, |
which implies that
∫TiTi−1s−rh(s)ds>Γ(qi)(4Ti−Ti−1)qi−1. |
So, we get
∫T0s−rh(s)ds>n∗∑i=1Γ(qi)(4Ti−Ti−1)qi−1. |
We complete the proof.
Remark 4.3. We notice that if r=0 and q(t)=q, q is a constant, i.e., BVP (1.1) is a fractional differential equation with constant order, then by similar arguments as done in [22], we get
∫T0h(s)ds>Γ(q)(4T)q−1. |
So, the inequalities (4.2) is a generalized Lyapunov-type inequality for the boundary value problem (1.1).
Example 5.1. Let us consider the following nonlinear boundary value problem
{Dq(t)0+x(t)+t−0.5|x|121+x2=0,0<t<2,u(0)=u(2)=0, | (5.1) |
where
q(t)={1.2, 0≤t≤1,1.6, 1<t≤2. |
We see that q(t) satisfies condition (H1); t0.5f(t,x(t))=|x(t)|121+x(t)2:[0,2]×R→R is continuous. Moreover, we have
t0.5|f(t,x(t))|=|x(t)|121+x(t)2≤|x(t)|12. |
Let r=0.5, c1=c2=1 and γ=12. We could verify that f(t,x)=t−0.5|x|121+x2 satisfies condition (H2). This suggests that the boundary value problem (5.1) has a solution by the conclusion of Theorem 3.3.
Example 5.2. Let us consider the following linear boundary value problem
{Dq(t)0+x(t)+t0.4=0,0<t<3,u(0)=0,u(3)=0, | (5.2) |
where
q(t)={1.2, 0≤t≤1,1.5, 1<t≤2,1.8, 2<t≤3. |
We see that q(t) satisfies condition (H1); f(t,x(t))=t0.4:[0,3]×R→R is continuous. Moreover, |f(t,x(t))|=t0.4≤30.4, thus we could take suitable constants to verify f(t,x)=t0.4 satisfies condition (H2). Then Theorem 3.3 assures the boundary value problem (5.2) has a solution.
In fact, we know that equation of (5.1) can been divided into three expressions as following
D1.20+x(t)+t0.4=0, 0<t≤1. | (5.3) |
For 1<t≤2,
d2dt2(∫10(t−s)−0.2Γ(0.8)x(s)ds+∫t1(t−s)−0.5Γ(0.5)x(s)ds)+t0.4=0. | (5.4) |
For 2<t≤3,
d2dt2(∫10(t−s)−0.2Γ(0.8)x(s)ds+∫21(t−s)−0.5Γ(0.5)x(s)ds+∫t2(t−s)−0.8Γ(0.2)x(s)ds)+t0.4=0. | (5.5) |
By [18], we can easily obtain that the following boundary value problems
{D1.20+x(t)+t0.4=0, 0<t≤1,x(0)=0,x(1)=0 |
{D1.51+x(t)=d2dt2∫t1(t−s)−0.5Γ(0.5)x(s)ds+t0.4=0, 1<t<2,x(1)=0,x(2)=0 |
{D1.82+x(t)=d2dt2∫t2(t−s)−0.8Γ(0.2)x(s)ds+t0.4=0, 2<t<3,x(2)=0,x(3)=0 |
respectively have solutions
x1(t)=Γ(1.4)Γ(2.6)(t0.2−t1.6)∈C[0,1];˜x2(t)=Γ(1.4)Γ(2.9)((t−1)0.5−(t−1)1.9)∈C[1,2];˜x3(t)=Γ(1.4)Γ(3.2)((t−2)0.8−(t−2)2.2)∈C[2,3]. |
It is known by calculation that
x1(t),0≤t≤1,x2(t)={0,0≤t≤1,˜x2(t),1<t≤2,x3(t)={0,0≤t≤2,˜x3(t),2<t≤3, | (5.6) |
are the solutions of (5.3)–(5.5), respectively. By Definition 3.2 and (5.6), we know that
x(t)={x1(t)=Γ(1.4)Γ(2.6)(t0.2−t1.6),0≤t≤1,x2(t)={0,0≤t≤1,Γ(1.4)Γ(2.9)((t−1)0.5−(t−1)1.9),1<t≤2,x3(t)={0,0≤t≤2,Γ(1.4)Γ(3.2)((t−2)0.8−(t−2)2.2),2<t≤3 |
is one solution of the boundary value problem (5.2).
In this paper, we consider a two-points boundary value problem of differential equations of variable order, which is a piecewise constant function. Based the essential difference about the variable order fractional calculus (derivative and integral) and the integer order and the constant fractional order calculus (derivative and integral), we carry on essential analysis to the boundary value problem (1.1). According to our analysis, we give the definition of solution to the boundary value problem (1.1). The existence result of solution to the boundary value problem (1.1) is derived. We present a Lyapunov-type inequality for the boundary value problems (1.1). Since the variable order fractional calculus (derivative and integral) and the integer order and the constant fractional order calculus (derivative and integral) has the essential difference, it is interesting and challenging about the existence, uniqueness of solutions, Lyapunov-type inequality, etc, to the boundary value problems of differential equations of variable order.
This research is supported by the Natural Science Foundation of China (11671181). The authors are thankful to the referees for their careful reading of the manuscript and insightful comments.
The author declares no conflicts of interest in this paper.
[1] |
M. Faieghi, S. Kuntanapreeda, H. Delavari, D. Baleanu, LMI-based stabilization of a class of fractional-order chaotic systems, Nonlinear Dyn., 72 (2013), 301–309. https://doi.org/10.1007/s11071-012-0714-6 doi: 10.1007/s11071-012-0714-6
![]() |
[2] |
F. Zhang, G. Chen, C. Li, J. Kurths, Chaos synchronization in fractional differential systems, Philos. T. R. Soc. A, 371 (2013), 20120155. https://doi.org/10.1098/rsta.2012.0155 doi: 10.1098/rsta.2012.0155
![]() |
[3] |
Y. Xu, W. Li, Finite-time synchronization of fractional-order complex-valued coupled systems, Phys. A, 549 (2020), 123903. https://doi.org/10.1016/j.physa.2019.123903 doi: 10.1016/j.physa.2019.123903
![]() |
[4] |
M. S. Ali, G. Narayanan, V. Shekher, A. Alsaedi, B. Ahmad, Global Mittag-Leffler stability analysis of impulsive fractional-order complex-valued BAM neural networks with time varying delays, Commun. Nonlinear Sci., 83 (2020), 105088. https://doi.org/10.1016/j.cnsns.2019.105088 doi: 10.1016/j.cnsns.2019.105088
![]() |
[5] |
Y. Ding, Z. Wang, H. Ye, Optimal control of a fractional-order HIV-immune system with memory, IEEE T. Contr. Syst. T., 20 (2012), 763–769. https://doi.org/10.1109/tcst.2011.2153203 doi: 10.1109/tcst.2011.2153203
![]() |
[6] |
R. Metzler, J. Klafter, The random walk's guide to anomalous diffusion: A fractional dynamics approach, Phys. Rep., 339 (2000), 1–77. https://doi.org/10.1016/s0370-1573(00)00070-3 doi: 10.1016/s0370-1573(00)00070-3
![]() |
[7] |
M. Javidi, B. Ahmad, Dynamic analysis of time fractional order phytoplankton-toxic phytoplankton-zooplankton system, Ecol. Model., 318 (2015), 8–18. https://doi.org/10.1016/j.ecolmodel.2015.06.016 doi: 10.1016/j.ecolmodel.2015.06.016
![]() |
[8] |
A. Carvalho, C. M. A. Pinto, A delay fractional order model for the co-infection of malaria and HIV/AIDS, Int. J. Dyn. Control, 5 (2017), 168–186. https://doi.org/10.1007/s40435-016-0224-3 doi: 10.1007/s40435-016-0224-3
![]() |
[9] |
J. Henderson, R. Luca, A. Tudorache, On a system of fractional differential equations with coupled integral boundary conditions, Fract. Calc. Appl. Anal., 18 (2015), 361–386. https://doi.org/10.1515/fca-2015-0024 doi: 10.1515/fca-2015-0024
![]() |
[10] |
J. R. Wang, Y. Zhang, Analysis of fractional order differential coupled systems, Math. Method. Appl. Sci., 38 (2015), 3322–3338. https://doi.org/10.1002/mma.3298 doi: 10.1002/mma.3298
![]() |
[11] |
L. Zhang, B. Ahmad, G. Wang, Monotone iterative method for a class of nonlinear fractional differential equations on unbounded domains in Banach spaces, Filomat, 31 (2017), 1331–1338. https://doi.org/10.2298/fil1705331z doi: 10.2298/fil1705331z
![]() |
[12] | B. Ahmad, A. Alsaedi, S. Aljoudi, S. K. Ntouyas, On a coupled system of sequential fractional differential equations with variable coeffcients and coupled integral boundary conditions, Bull. Math. Soc. Sci. Math. Roumanie, 60 (2017), 3–18. |
[13] |
M. S. Abdo, K. Shah, S. K. Panchal, H. A. Wahash, Existence and Ulam stability results of a coupled system for terminal value problems involving ψ-Hilfer fractional operator, Adv. Differ. Equ., 2020 (2020), 316. https://doi.org/10.1186/s13662-020-02775-x doi: 10.1186/s13662-020-02775-x
![]() |
[14] |
A. M. Saeed, M. S. Abdo, M. B. Jeelani, Existence and Ulam-Hyers stability of a fractional-order coupled system in the frame of generalized Hilfer derivatives, Mathematics, 9 (2021), 2543. https://doi.org/10.3390/math9202543 doi: 10.3390/math9202543
![]() |
[15] |
B. Ahmad, R. Luca, Existence of solutions for a sequential fractional integro-differential system with coupled integral boundary conditions, Chaos Soliton. Fract., 104 (2017), 378–388. https://doi.org/10.1016/j.chaos.2017.08.035 doi: 10.1016/j.chaos.2017.08.035
![]() |
[16] | R. S. Adigüzel, U. Aksoy, E. Karapinar, I. M. Erhan, On the solutions of fractional differential equations via geraghty type hybrid contractions, Appl. Comput. Math., 20 (2021), 313–333. |
[17] |
A. Salim, M. Benchohra, E. Karapinar, J. E. Lazreg, Existence and Ulam stability for impulsive generalized Hilfer-type fractional differential equations, Adv. Differ. Equ., 2020 (2020), 601. https://doi.org/10.1186/s13662-020-03063-4 doi: 10.1186/s13662-020-03063-4
![]() |
[18] | R. S. Adigüzel, U. Aksoy, E. Karapinar, I. M. Erhan, On the solution of a boundary value problem associated with a fractional differential equation, Math. Method. Appl. Sci., 2020. https://doi.org/10.1002/mma.6652 |
[19] | K. Diethelm, The analysis of fractional differential equations, Springer, 2010. https://doi.org/10.1007/978-3-642-14574-2 |
[20] | A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Preface, North-Holland Math. Stud., 204 (2006), vii-x. https://doi.org/10.1016/s0304-0208(06)80001-0 |
[21] | V. Lakshmikantham, S. Leela, J. V. Devi, Theory of fractional dynamic systems, Cambridge Scientific Publishers, 2009. |
[22] | K. S. Miller, B. Ross, An introduction to the fractional calculus and differential equations, New York: John Wiley, 1993. |
[23] | I. Podlubny, Fractional differential equations, New York: Academic Press, 1999. |
[24] | S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional integrals and derivatives, Yverdon: Gordon and Breach Science, 1993. |
[25] | B. Ahmad, A. Alsaedi, S. K. Ntouyas, J. Tariboon, Hadamard-type fractional differential equations, inclusions and inequalities, Switzerland: Springer, 2017. https://doi.org/10.1007/978-3-319-52141-1 |
[26] | B. Ahmad, S. K. Ntouyas, Nonlocal nonlinear fractional-order boundary value problems, Singapore: World Scientific, 2021. https://doi.org/10.1142/12102 |
[27] | Y. Zhou, Basic theory of fractional differential equations, Singapore: World Scientific, 2014. https://doi.org/10.1142/9069 |
[28] |
K. D. Kucche, A. D. Mali, On the nonlinear (k,φ)-Hilfer fractional differential equations, Chaos Soliton. Fract., 152 (2021), 111335. https://doi.org/10.1016/j.chaos.2021.111335 doi: 10.1016/j.chaos.2021.111335
![]() |
[29] |
S. K. Ntouyas, B. Ahmad, J. Tariboon, M. S. Alhodaly, Nonlocal integro-multi-point (k,ψ)-Hilfer type fractional boundary value problems, Mathematics, 10 (2022), 2357. https://doi.org/10.3390/math10132357 doi: 10.3390/math10132357
![]() |
[30] |
A. Samadi, S. K. Ntouyas, J. Tariboon, Nonlocal coupled system for (k,φ)-Hilfer fractional differential equations, Fractal Fract., 6 (2022), 234. https://doi.org/10.3390/fractalfract6050234 doi: 10.3390/fractalfract6050234
![]() |
[31] |
T. Li, A class of nonlocal boundary value problems for partial differential equations and its applications in numerical analysis, J. Comput. Appl. Math., 28 (1989), 49–62. https://doi.org/10.1016/0377-0427(89)90320-8 doi: 10.1016/0377-0427(89)90320-8
![]() |
[32] |
B. Ahmad, J. J. Nieto, Existence of solutions for nonlocal boundary value problems of higher-order nonlinear fractional differential equations, Abstr. Appl. Anal., 2009 (2009), 494720. https://doi.org/10.1155/2009/494720 doi: 10.1155/2009/494720
![]() |
[33] |
A. Alsaedi, B. Ahmad, S. Aljoudi, S. K. Ntouyas, A study of a fully coupled two-parameter system of sequential fractional integro-differential equations with nonlocal integro-multipoint boundary conditions, Acta Math. Sci., 39 (2019), 927–944. https://doi.org/10.1007/s10473-019-0402-4 doi: 10.1007/s10473-019-0402-4
![]() |
[34] |
R. S. Adigüzel, U. Aksoy, E. Karapinar, I. M. Erhan, Uniqueness of solution for higher-order nonlinear fractional differential equations with multi-point and integral boundary conditions, RACSAM, 115 (2021), 155. https://doi.org/10.1007/s13398-021-01095-3 doi: 10.1007/s13398-021-01095-3
![]() |
[35] |
Y. Alruwaily, B. Ahmad, S. K. Ntouyas, A. S. M. Alzaidi, Existence results for coupled nonlinear sequential fractional differential equations with coupled Riemann–Stieltjes integro-multipoint boundary conditions, Fractal Fract., 6 (2022), 123. https://doi.org/10.3390/fractalfract6020123 doi: 10.3390/fractalfract6020123
![]() |
[36] |
J. Tariboon, A. Samadi, S. K. Ntouyas, Multi-point boundary value problems for (k,ψ)-Hilfer fractional differential equations and inclusions, Axioms, 11 (2022), 110. https://doi.org/10.3390/axioms11030110 doi: 10.3390/axioms11030110
![]() |
[37] |
Y. C. Kwun, G. Farid, W. Nazeer, S. Ullah, S. M. Kang, Generalized Riemann-Liouville k-fractional integrals associated with Ostrowski type inequalities and error bounds of Hadamard inequalities, IEEE Access, 6 (2018), 64946–64953. https://doi.org/10.1109/access.2018.2878266 doi: 10.1109/access.2018.2878266
![]() |
[38] | K. Deimling, Nonlinear functional analysis, New York: Springer, 1985. https://doi.org/10.1007/978-3-662-00547-7 |
[39] | A. Granas, J. Dugundji, Fixed point theory, New York: Springer-Verlag, 2005. https://doi.org/10.1007/978-0-387-21593-8 |
[40] | M. A. Krasnoselski, Two remarks on the method of successive approximations, Uspekhi Mat. Nauk, 10 (1955), 123–127. |
1. | Ibtesam Alshammari, Islam M. Taha, On fuzzy soft β-continuity and β-irresoluteness: some new results, 2024, 9, 2473-6988, 11304, 10.3934/math.2024554 | |
2. | D. I. Taher, R. Abu-Gdairi, M. K. El-Bably, M. A. El-Gayar, Decision-making in diagnosing heart failure problems using basic rough sets, 2024, 9, 2473-6988, 21816, 10.3934/math.20241061 | |
3. | Fahad Alsharari, Ahmed O. M. Abubaker, Islam M. Taha, On r-fuzzy soft γ-open sets and fuzzy soft γ-continuous functions with some applications, 2025, 10, 2473-6988, 5285, 10.3934/math.2025244 | |
4. | Fahad Alsharari, Hind Y. Saleh, Islam M. Taha, Some Characterizations of k-Fuzzy γ-Open Sets and Fuzzy γ-Continuity with Further Selected Topics, 2025, 17, 2073-8994, 678, 10.3390/sym17050678 |