Next Article in Journal / Special Issue
Letnikov vs. Marchaud: A Survey on Two Prominent Constructions of Fractional Derivatives
Previous Article in Journal
A Numerical Solution of Fractional Lienard’s Equation by Using the Residual Power Series Method
Previous Article in Special Issue
On Some New Properties of the Fundamental Solution to the Multi-Dimensional Space- and Time-Fractional Diffusion-Wave Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Iterative Method for Solving a Class of Fractional Functional Differential Equations with “Maxima”

1
Department of Mathematics, Faculty of Exact Sciences, Echahide Hamma Lakhdar University, B.P. 789, El Oued 39000, Algeria
2
Laboratory of Applied Mathematics, Badji Mokhtar University, B.P. 12, Annaba 23000, Algeria
*
Author to whom correspondence should be addressed.
Mathematics 2018, 6(1), 2; https://doi.org/10.3390/math6010002
Submission received: 14 November 2017 / Revised: 16 December 2017 / Accepted: 19 December 2017 / Published: 22 December 2017
(This article belongs to the Special Issue Fractional Calculus: Theory and Applications)

Abstract

:
In the present work, we deal with nonlinear fractional differential equations with “maxima” and deviating arguments. The nonlinear part of the problem under consideration depends on the maximum values of the unknown function taken in time-dependent intervals. Proceeding by an iterative approach, we obtain the existence and uniqueness of the solution, in a context that does not fit within the framework of fixed point theory methods for the self-mappings, frequently used in the study of such problems. An example illustrating our main result is also given.

1. Introduction

One of the most interesting kinds of nonlinear functional differential equations is the case when the nonlinear part depends on the maximum values of the unknown function. These equations, called functional differential equations with “maxima”, arise in many technological processes. For instance, in the automatic control theory of various technical systems, it occurs that the law of regulation depends on the maximal deviation of the regulated quantity (see [1,2]). Such problems are often modeled by differential equations that contain the maximum values of the unknown function (see [3,4,5]). Recently, ordinary differential equations with “maxima” have received wide attention and have been investigated in diverse directions (see, for example, [4,6,7,8,9,10,11] and the references therein). As far as we know, in the fractional case, these equations are not yet sufficiently discussed in the existing literature, and thus form a natural subject for further investigation. Motivated by the previous fact and inspired by [11], in this work, we focus on the existence and uniqueness of the solution for similar systems in a fractional context, and in more general terms. We consider the following nonlinear fractional differential equation with “maxima” and deviating arguments:
C D α u t = f t , max σ [ a ( t ) , b ( t ) ] u σ , u t τ 1 t , , u t τ N t , t > 0 ,
with the initial condition function
u t = ϕ t , t 0 ,
where C D α denotes the Caputo fractional derivative operator of order α 0 , 1 , N is a positive integer, a, b and τ i (with 1 i N ) are real continuous functions defined on R + = [ 0 , + ] subject to conditions that will be specified later, ϕ : [ , 0 ] R is a continuous function such that ϕ 0 = ϕ 0 > 0 , and f : R + × R 1 + N R is a nonlinear continuous function.
Our aim is to give sufficient assumptions leading to an iterative process that converges to the unique continuous solutions of Equations (1) and (2). These being under weaker conditions compared to the usual contractions (see Remark 3), and in a setting for which the standard process of Picard’s iterations fails to be well defined.
It should be pointed out here, that the maximums in Equation (1) are taken on time-dependent intervals and not on a fixed one as is the case of the example given in [11].
Moreover, the Equation (1) will be supposedly of mixed type, namely with both retarded and advanced deviations τ i , while, in [11], only the delays are considered. It is also important to note that, in the Lipschitz condition of the nonlinear function f, we take into account the direction of maximums too, which is not the case of the corresponding assumption in [11].
Due to all of these generalizations, our work attempts to extend the application of [11] (Theorem 3) to the fractional case by a constructive approach.
To our knowledge, the studies devoted to the question of the existence and uniqueness of the solutions for fractional differential equations are based on different variants from the fixed point theory for self-mappings, or on the upper and lower solutions method (see, e.g., [12,13,14,15,16,17,18] and the references therein). We emphasize here that our result answers this question for a class of problems of the forms Equations (1) and (2), even when the previous versions of the theory fail to do so directly. That is, when the integral operator associated with Equations (1) and (2) is allowed to be a non self-mapping (see Remark 1).
The rest of the paper is organized as follows. In the next section, we introduce some basic definitions from the fractional calculus as well as preliminary lemmas. In Section 3, under some sufficient conditions allowing the integral operator associated with Equations (1) and (2) to be non self, we prove an existence–uniqueness result by means of an iterative process. The applicability of our theoretical result is illustrated in Section 4.

2. Preliminaries

We start by recalling the definitions of the Riemann–Liouville fractional integrals and the Caputo fractional derivatives on the half real axis. For further details on the historical account and essential properties about the fractional calculus, we refer to [19,20,21,22].
Definition 1.
The Riemann–Liouville fractional integral of a function u : R + R of order α R + is defined by
I α u t : = 1 Γ α 0 t t s α 1 u s d s , t > 0 ,
where Γ . is the Gamma function, provided that the right side is pointwise defined on 0 , .
In the following definition, n denotes the positive integer such that n 1 < α n and d n / d t n is the classical derivative operator of order n. For simplicity, we set d u / d t = u t .
Definition 2.
The Caputo fractional derivative of a function u : R + R of order α R + is defined by
C D α u t : = I n α d n d t n u t : = 1 Γ n α 0 t t s n α 1 d n d s n u s d s , t > 0 ,
provided that the right-hand side exists pointwise on 0 , .
In particular, when 0 < α < 1 ,
C D α u t : = I 1 α u t : = 1 Γ 1 α 0 t u s t s α d s , t > 0 .
Let us denote by C R the set of all real continuous functions on R . Applying the Riemann–Liouville fractional integral operator I α of order α to both sides of Equation (1) and using its properties (see [19,21]), together with the initial condition Equation (2), we easily get the following lemma.
Lemma 1.
If f , a , b and τ i (with 1 i N ) are continuous functions, then u v C R s . t . v t = ϕ t f o r t 0 is a solution of Equations (1) and (2) if and only if u t = F u t , where
F u t = ϕ 0 + 0 t t s α 1 Γ α f s , max σ [ a ( s ) , b ( s ) ] u σ , u s τ 1 s , , u s τ N s d s , t > 0 ,
F u t = ϕ t , t 0 .
Proof. 
Let u C R . The functions a and b are continuous, so, according to the Remark in [7] (page 8) , see also [4] (Remark 3.1.1, page 62), max σ [ a ( t ) , b ( t ) ] u σ is continuous too. Moreover, since τ i are continuous, then f t , max σ [ a ( t ) , b ( t ) ] u σ , u t τ 1 t , , u t τ N t as a composition of continuous functions, it is also continuous. Now, we are able to follow the usual approach to show this type of result (see [15,19,21,23,24]). Note first that the Caputo fractional derivative of order α 0 , 1 can be expressed by means of the Riemann–Liouville fractional derivative denoted by D α , as follows (see [21] (2.4.4) or [19] (Definition 3.2)):
C D α u t = D α u t u 0 : = d d t I 1 α u t u 0 .
Let now u v C R s . t . v t = ϕ t for t 0 be a solution of Equations (1) and (2). Thus, in view of the first equality in Equation (6), Equation (1) can be rewritten as
D α u t u 0 = f t , max σ [ a ( t ) , b ( t ) ] u σ , u t τ 1 t , , u t τ N t , t > 0 .
Since the right-hand side of Equation (7) is continuous, then according to the definition of the Riemann–Liouville fractional derivative given by the second equality in Equation (6), we have
I 1 α u t u 0 C 1 R + .
Thus, using [21] (Lemma 2.9, (d) with γ = 0 ), we have
I α D α u t u 0 = u t u 0 1 Γ α I 1 α U 0 t α 1 ,
where
U t : = u t u 0 .
Since U is continuous, for every T > 0 , there exists L > 0 such that U t L : for all t 0 , T . Thus, the following inequality holds true for every t > 0 , sufficiently small
I 1 α U t 1 Γ 1 α 0 t U s t s α d s L Γ 2 α t 1 α .
Hence, the fact that 1 α > 0 , together with the continuity of I 1 α U resulting from Equation (8), imply that I 1 α U 0 = 0 . Consequently, Equation (9) becomes
I α D α u t u 0 = u t u 0 .
Now, returning to Equation (7), applying the Riemann–Liouville fractional integral to both sides, and then using Equation (11) together with Equation (2), we obtain Equation (4).
Suppose now that u v C R s . t . v t = ϕ t for t 0 is a solution of Equations (4) and (5). Then, in view of Definition 1, we can rewrite Equation (4) as
u t = ϕ 0 + I α f t , max σ [ a ( t ) , b ( t ) ] u σ , u t τ 1 t , , u t τ N t .
Since u is continuous, then the right-hand side above is continuous too. By applying the Caputo fractional derivative operator C D α to both sides, then using its linearity (see [19] (Theorem 3.16)), as well as the fact that the derivative of a constant (in the sense of Caputo) is equal to zero [21] (Property 2.16), together with [21] (Lemma 2.21), we get Equation (1). ☐
In the present work, the state space will be regarded as a complete Hausdorff locally convex space. For further details on these spaces, we refer to [25]. In the sequel of this paper, we make use of the following lemma, which can be found in [26] ([Lemma 2).
Lemma 2.
Let X be a complete Hausdorff locally convex space, E a closed subset of X and u , v X . If u E and v E , then there exists β 0 , 1 such that w β : = ( 1 β ) u + β v E , where E denotes the boundary of E . Furthermore, if u E , then β 0 , 1 .

3. The Main Results

In this section, we not only prove the existence–uniqueness result for Equations (1) and (2), but we also give this solution as a limit of an iterative process.
First, let us set the following hypotheses:
H 1
t 0 : 0 a a ( t ) b ( t ) b , with a : = inf t [ 0 , + ] a ( t ) , and b : = sup t [ 0 , + ] b ( t ) . Furthermore, for all t 0 , b , we assume that a ( t ) = a and b ( t ) = b . In other words, the functions a and b are constant on the interval 0 , b .
H 2
τ > 0 , such that, for i = 1 , , N : τ i t > t τ , t > 0 .
H 3
For i = 1 , , N , t i > 0 : τ i t t , t [ 0 , t i ] , and τ i t < t , t t i , + .
H 4
There exist positive constants l 1 and l 2 , such that f satisfies the Lipschitz condition
f t , ξ , x 1 , , x N f t , η , y 1 , , y N l 1 ξ η + l 2 i = 1 N x i y i .
H 5
There exists a positive constant M > ϕ 0 such that
1 Γ α α f t , M , x 1 , , x N M ϕ 0 b α , t , x 1 , , x N 0 , b × R N .
H 6
f is a non negative function, and, moreover, h [ ϕ 0 , M ] such that t 0 , b
1 Γ α α f t , h , x 1 , , x N > M ϕ 0 | b max 1 i N t i | α , x 1 , , x N [ ϕ 0 , ϕ 0 + h ϕ 0 b τ ] N .
Let X = C R be the locally convex sequentially complete Hausdorff space of all real valued continuous functions defined on R , and { P K : K K } be the saturated family of semi-norms, generating the topology of X, defined by
P K u = sup t K e λ t u t ,
where K runs over the set of all compact subsets of R denoted by K , and λ is a positive real number to be specified later.
We denote by E ϕ , M , the subset of X defined by
E ϕ , M = u X : u t = ϕ t for t 0 , and u t M for t [ a , b ] ,
where a , b and M are the constants given by H 1 and H 5 . It can be easily seen that E ϕ , M is a closed subset of X and its boundary is
E ϕ , M = u X : u t = ϕ t for t 0 and max t [ a , b ] u t = M .
Throughout the remaining of this paper, F denotes the operator defined on E ϕ , M by Equations (4) and (5). Thus, according to Lemma 1, F maps E ϕ , M into X and the fixed points of F are continuous solutions of problems Equations (1) and (2).
Remark 1.
It should be pointed out that under hypotheses H 1 H 3 , H 6 with the additional condition max 1 i N t i < b , F is a non-self mapping on E ϕ , M . Indeed, as is noted in the proof of [11] (Theorem 3) , for any function u E ϕ , M defined by u t = ϕ 0 + h ϕ 0 t / b , where t 0 , b and h is the constant given by H 6 , it can be easily seen that F u E ϕ , M . This will be checked by the example of the last section.
The introduction of a self-mapping of the index set in uniform spaces is motivated by applications in the theory of neutral functional differential equations [11,27,28]. Following this idea, let us define a map j : K K by
j ( K ) : = K , if K + = , 0 , max { K m , τ , b } , if K + ,
where K + : = K 0 , + , K m = sup K , τ and b are the positive constants given in H 1 H 2 . For n N , j n ( K ) is the compact set defined inductively by j n ( K ) = j j n 1 ( K ) and j 0 ( K ) = K .
Remark 2.
Note that, for every K K and every integer n greater than 1, we have j n ( K ) = j ( K ) .
In the next proposition, we show that F satisfies Equation (14), which is a weakened version of the usual contraction when L λ < 1 (see Remark 3).
Proposition 1.
Under hypotheses H 1 H 4 , the operator F : E ϕ , M X satisfies for each u , v E ϕ , M and every K K
P K F x F y L λ P j ( K ) x y
with
L λ = l 1 λ α Γ α Γ α 2 1 1 + α α 1 + α α 1 + α e λ b + N l 2 e λ τ λ α .
Proof. 
Note that it suffices to consider K + , since otherwise P K F u F v = 0 . Letting t K + , we obtain by means of hypotheses H 3 and H 4
F u t F v t 0 t t s α 1 Γ α l 1 | max σ [ a ( s ) , b ( s ) ] u σ max σ [ a ( s ) , b ( s ) ] v σ | d s
+ 0 t t s α 1 Γ α + l 2 i = 1 N | u s τ i s v s τ i s | d s
l 1 0 t t s α 1 Γ α max σ [ a ( s ) , b ( s ) ] | u σ v σ | d s + l 2 i = 1 N 0 min { t i , t } t s α 1 Γ α ϕ r i s ϕ r i s d s
+ l 2 i { 1 , , N : t i t } t i t t s α 1 Γ α u r i s v r i s d s
= l 1 0 t t s α 1 Γ α max σ [ a ( s ) , b ( s ) ] | u σ v σ | d s + l 2 i { 1 , , N : t i t } t i t t s α 1 Γ α u r i s v r i s d s
l 1 0 t t s α 1 Γ α e λ b ( s ) max σ [ a ( s ) , b ( s ) ] e λ σ | u σ v σ | d s
+ l 2 i { 1 , , N : t i t } t i t t s α 1 Γ α e λ r i ( s ) e λ r i ( s ) u r i s v r i s d s
l 1 max σ [ a , b ] e λ σ | u σ v σ | 0 t t s α 1 Γ α e λ b ( s ) d s
+ l 2 i { 1 , , N : t i t } t i t t s α 1 Γ α e λ r i ( s ) e λ r i ( s ) u r i s v r i s d s ,
where r i s = s τ i s . Note that, due to the definition Equation (13) and under hypothesis H 1 , it is clear that, for every K K with K + , we have a , b j ( K ) and further H 2 H 3 lead to r i s j ( K ) when t i s t . Hence,
F x t F y t l 1 P j ( K ) u v 0 t t s α 1 Γ α e λ b ( s ) d s
+ l 2 i { 1 , , N : t i t } t i t t s α 1 Γ α e λ r i ( s ) max ξ j ( K ) e λ ξ | u ξ v ξ | d s
= l 1 P j ( K ) u v 0 t t s α 1 Γ α e λ b ( s ) d s + l 2 P j ( K ) u v i { 1 , , N : t i t } t i t t s α 1 Γ α e λ r i ( s ) d s .
Now, multiplying the both sides of the above inequality by e λ t , then performing the change of variable u = λ t s , we get
e λ t F u t F v t l 1 P j ( K ) u v 0 t t s α 1 Γ α e λ ( t b ( s ) ) d s
+ l 2 P j ( K ) u v i { 1 , , N : t i t } t i t t s α 1 Γ α e λ t r i s d s
= l 1 λ α Γ α p j ( K ) u v 0 λ t x α 1 e λ ( t b ( t x λ ) ) d x
+ l 2 λ α P j ( K ) u v i { 1 , , N : t i t } 0 λ t t i x α 1 Γ α e x e λ τ i t x λ d x
l 1 λ α Γ α P j ( K ) u v 0 λ t x α 1 e x e λ ( ( t x λ ) b ( t x λ ) ) d x
+ l 2 λ α P j ( K ) u v i { 1 , , N : t i t } 0 λ t t i x α 1 Γ α e x e λ τ i t x λ d x .
Let μ : = 1 + α and ν : = 1 + 1 / α . Taking into account H 3 , Hölder’s inequality gives
e λ t F u t F v t ,
l 1 λ α Γ α 0 λ t x μ α 1 e μ x d x 1 μ 0 λ t e ν λ ( ( t x λ ) b ( t x λ ) ) d x 1 ν + N l 2 e λ τ λ α P j ( K ) u v
= l 1 λ α Γ α 0 λ t x μ α 1 e μ x d x 1 μ λ 0 t e ν λ ( s b ( s ) ) d s 1 ν + N l 2 e λ τ λ α P j ( K ) u v
l 1 λ α Γ α 0 λ t x μ α 1 e μ x d x 1 μ λ 0 t e ν λ ( s b ) d s 1 ν + N l 2 e λ τ λ α P j ( K ) u v
l 1 λ α Γ α Γ α 2 1 μ 1 ν 1 ν e λ b + N l 2 e λ τ λ α P j ( K ) u v .
Thus, the result is obtained by taking the supremum on K. ☐
Remark 3.
Since K + j ( K ) , if P K F u F v L λ P K u v is satisfied, then Equation (14) holds true. Therefore, due to the choice of j, in the present context, the usual contraction is a particular case of Equation (14) when L λ < 1 .
To reach our aim, we proceed by adapting the proof of [11], [Theorem 1] with some completeness, for the construction of an iterative process converging to the unique continuous solutions of Equations (1) and (2).
According to Remark 1, the standard process of Picard’s iterations fails to be well defined. To overcome this fact, we make use of Lemma 2 to construct a sequence of elements of E ϕ , M as follows: starting from an arbitrary point u 0 E ϕ , M , we define the terms of a sequence { u n } n N in E ϕ , M iteratively as follows:
u n = F u n 1 , if F u n 1 E ϕ , M , u n = ( 1 β n ) u n 1 + β n F u n 1 E ϕ , M with β n 0 , 1 , if F u n 1 E ϕ , M .
Note that the terms of the sequence { u n } n N belong to A B E ϕ , M , with B E ϕ , M , where
A : = { u i { u n } n N : u i = F u i 1 } and B : = { u i { u n } n N : u i F u i 1 } .
Furthermore, if u n B , a straightforward computation leads to
P K ( u n 1 u n ) + P K ( u n F u n 1 ) = P K ( u n 1 F u n 1 ) K K .
Proposition 2.
Let u 0 E ϕ , M , and { u n } n N be the sequence defined iteratively by Equation (16). Then, under hypotheses H 1 H 5 , for each K K and every integer m greater than or equal to 1, the following estimation holds true:
max { P K ( u 2 m u 2 m + 1 ) , P K ( u 2 m + 1 u 2 m + 2 ) } 2 2 m 1 L λ m C K ,
where L λ is given by Equation (15) and
C K = max { P j ( K ) ( u 0 u 1 ) , P j ( K ) ( u 1 u 2 ) , P j ( K ) ( u 2 u 3 ) } .
Proof. 
If max t [ a , b ] u t = M , then, for a t b , hypothesis H 5 gives
F u t = ϕ 0 + 1 Γ α 0 t t s α 1 f s , M , u s τ 1 s , , u s τ N s d s ϕ 0 + α ( M ϕ 0 ) b α 0 t t s α 1 d s M ,
which means that F ( E ϕ , M ) E ϕ , M . Consequently, two consecutive terms of the sequence { u n } n N can not belong to B (recall that B E ϕ , M ). Thus, it suffices to consider the three cases below.
Case 1.
u n , u n + 1 A . From Equation (14), we have
P K ( u n u n + 1 ) = P K ( F u n 1 F u n ) L λ · P j ( K ) ( u n 1 u n ) .
Case 2.
u n A , u n + 1 B . From the condition Equation (14) together with Equation (17) (for u n + 1 instead of u n ), we get
P K ( u n u n + 1 ) = P K ( u n F u n ) P K ( u n + 1 F u n ) P K ( F u n 1 F u n ) L λ · P j ( K ) ( u n 1 u n ) .
Case 3.
u n B , u n + 1 A . Then, β n 0 , 1 : u n = ( 1 β n ) u n 1 + β n F u n 1 , which implies that
P K ( u n u n + 1 ) max { P K ( u n 1 u n + 1 ) , P K ( F u n 1 u n + 1 ) } .
Thus, by Equation (14), for every integer number n 2 , we obtain either
P K ( u n u n + 1 ) L λ · p j ( K ) ( u n 1 u n ) , or P K ( u n u n + 1 ) L λ · P j ( K ) ( u n 2 u n ) .
Moreover,
P j ( K ) ( u n 2 u n ) P j ( K ) ( u n 2 u n 1 ) + P j ( K ) ( u n 1 u n ) 2 max { P j ( K ) ( u n 2 u n 1 ) , P j ( K ) ( u n 1 u n ) }
In summary, the following inequality is true in all cases
P K ( u n u n + 1 ) 2 L λ max { P j ( K ) ( u n 2 u n 1 ) , P j ( K ) ( u n 1 u n ) } .
We now prove Equation (18) by induction. Using Equation (19), we have either
P K ( u 2 u 3 ) 2 L λ · P j ( K ) ( u 0 u 1 ) 2 L λ · C K ,
or
P K ( u 2 u 3 ) 2 L λ · P j ( K ) ( u 1 u 2 ) 2 L λ · C K ,
and similarly we obtain
P K ( u 3 u 4 ) 2 L λ · C K .
Consequently, Equation (18) is satisfied for m = 1 . Assume now that Equation (18) holds true for some m > 1 . Using Equation (19), we get either
P K ( u 2 m + 2 u 2 m + 3 ) 2 L λ · P j ( K ) ( u 2 m u 2 m + 1 ) ,
or
P K ( u 2 m + 2 u 2 m + 3 ) 2 L λ · P j ( K ) ( u 2 m + 1 u 2 m + 2 ) .
Thus, the fact that C j ( K ) = C K , which follows from Remark 2, leads to
P K ( u 2 m + 2 u 2 m + 3 ) 2 L λ 2 2 m 1 L λ m C j ( K ) = 2 2 m L λ m + 1 C K 2 2 ( m + 1 ) 1 L λ m + 1 C K .
In the same way, we get
P K ( u 2 m + 3 u 2 m + 4 ) 2 2 ( m + 1 ) 1 L λ m + 1 C K ,
which means that Equation (18) holds for m + 1 , and this completes the proof. ☐
We are now ready to prove our main result.
Theorem 1.
Let u 0 E ϕ , M , then under hypotheses H 1 H 6 with
max 1 i N t i < b ,
the sequence { u n } n N defined iteratively by Equation (16), converges in E ϕ , M to the unique continuous solution of Equations (1) and (2) provided that
l 1 Γ α Γ α 2 1 1 + α α 1 + α α 1 + α + N l 2 e max { τ , b } α < 1 4 .
Proof. 
Let us put λ = 1 / max { τ , b } in Equation (12). Thus, according to Proposition 1, for every K K and u , v E ϕ , M , Equation (14) holds true with L λ < 1 / 4 . Therefore, for an arbitrary fixed K K , and, for each ε > 0 , there exists a positive integer s satisfying
m = s 2 2 m L λ m < ε C K .
Hence, for n 2 s , q 1 and a sufficiently large l, we get, by means of Equations (18) and (22),
P K ( u n u n + q ) P K ( u n u n + 1 ) + P K ( u n + 1 u n + 2 ) + + P K ( u n + q 1 u n + q ) m = s l { P K ( u 2 m u 2 m + 1 ) + P K ( u 2 m + 1 u 2 m + 2 ) } m = s l 2 2 m L λ m · C K C K · m = s 2 2 m L λ m < ε .
Consequently, { u n } n N is a Cauchy sequence in the closed subset E ϕ , M of the complete locally convex space X, and so it converges to a point u E ϕ , M . Let { u n k } k 1 be a sub-sequence of { u n } n 1 in A , which is u n k + 1 = F u n k for every positive integer k. Then, for each compact K K , we have
P K ( u F u ) P K ( u u n k ) + P K ( u n k F u ) = P K ( u u n k ) + P K ( F u n k 1 F u ) P K ( u u n k ) + L λ · P j ( K ) ( u n k 1 u ) k + 0 .
Therefore, u = F u and so, according to Lemma 1, u is a solution of Equations (1) and (2). For the uniqueness, assume that there exists another solution v E ϕ , M such that u v . Since X is Hausdorff, then P K 0 ( u v ) 0 for some compact K 0 K . Using Equation (14) and Remark 2, we get for every positive integer n
0 < P K 0 ( u v ) = P K 0 ( F u F v ) L λ · P j ( K 0 ) ( u v ) = L λ · P j ( K 0 ) ( F u F v ) L λ 2 · P j ( K 0 ) ( u v ) L λ n P j ( K 0 ) ( u v ) ,
which contradicts the fact that L λ < 1 4 . This completes the proof. ☐

4. Example

The following example illustrates the applicability of our theoretical result. Let us consider the following equation
C D 0.5 u t = 1.132 2.45 + 10 2 | max t [ 10 1 , 2 ] u t | + 10 4 | u 0.994 1 + t t 0.004 | , t > 0 ,
subject to the initial condition function
u t = 2 t 2 + 1.168 , t 0 ,
Problems in Equations (23) and (24) are identified to Equations (1) and (2) with α = 0.5 , N = 1 , a ( t ) = 10 1 , b ( t ) = 2 , τ 1 t = t 0.994 1 + t t + 0.004 , ϕ t = 2 t 2 + 1.168 and
f t , ξ , η = 1.132 2.45 + 10 2 ξ + 10 4 η .
It can be easily seen that hypotheses H 1 H 4 are satisfied with a = 10 1 , b = 2 , τ = 0.994 , t 1 = 4 × 10 3 / 1 10 2 , l 1 = 1132 × 10 5 and l 2 = 1132 × 10 7 .
In addition, there exists M = 1.9 ( M > ϕ 0 = 1.168 ), such that, for every η R , we have
f t , M , η 1.132 2.45 + 10 2 × M 0.45845216686918 ,
and
0.5 Γ 0.5 M ϕ 0 b 0.5 0.458712974257473 .
Thus, H 5 holds true. To check H 6 , let h = 1.168001 ( ϕ 0 < h < M ), and then we have for every η [ ϕ 0 , ϕ 0 + h ϕ 0 b τ [ = 1.168 , 1.168 + 9.97 × 10 5
f t , h , η > 1.132 2.45 + h × 10 2 + ϕ 0 + h ϕ 0 b τ × 10 4 0.461798645149363 ,
and
0.5 Γ 0.5 M ϕ 0 b t 1 0.5 0.459177023920154 .
Since, moreover, f is clearly non negative, hypothesis H 6 is satisfied too.
Furthermore, we have t 1 < b and
l 1 Γ 0.5 Γ 0.25 1 1.5 0.5 1.5 0.5 1.5 + N l 2 e max { τ , b } 0.5 0.041231690678434 ,
which is all conditions of Theorem 1 being fulfilled. Now, let u 0 be the function defined by
u 0 t = 1.168 + ( 5 × 10 7 ) t : t 0 , b , 1.168001 : t b , 2 t 2 + 1.168 : t 0 .
It is clear that u 0 E ϕ , 1.9 , where
E ϕ , 1.9 = u C R : u t = 2 t 2 + 1.168 for t 0 and u t 1.9 for t 10 1 , 2 .
Note that, for 0 < t τ 1 t < τ , we have
ϕ 0 = 1.168 < u 0 t τ 1 t < 1.168 + ( 5 × 10 7 ) τ = ϕ 0 + h ϕ 0 b τ .
Then, hypothesis H 6 yields to
F u 0 2 = F u 0 b = ϕ 0 + 0 b b s α 1 Γ α f s , h , u 0 s τ 1 s d s ϕ 0 + t 1 b b s α 1 Γ α f s , h , u 0 s τ 1 s d s > ϕ 0 + α M ϕ 0 b t 1 α t 1 b b s α 1 d s = M = 1.9 ,
which means that F u 0 E ϕ , 1.9 . Thus, in this framework, the iterative processes usually used in the self-mapping context can not be applied, while, according to Theorem 1, the process defined by Equation (16), converges in E ϕ , 1.9 , to the unique continuous solutions of Equations (23) and (24). The first term is approximately given by
u 1 t 1.168 + ( 12115 × 10 13 ) t + 0.5176017180 t : 0 < t t 1 , 1.168 + ( 12115 × 10 13 ) t + ( 509761847 × 10 18 ) 27225 t 110 + 0.5176017180 t : t 1 < t b , 1.168000002 + ( 509761847 × 10 18 ) 27225 t 110 + 0.517601718 t : t b , 2 t 2 + 1.168 : t 0 .
For t > 0 , the other terms can be computed using the following formulas:
u 2 t = ϕ 0 + 1.132 Γ α 0 t t s α 1 2.45 + 10 2 | max t [ 10 1 , 2 ] u 1 s | + 10 4 | u 1 0.994 1 + s s 0.004 | d s .
By successive iterations up to the order n 1 , the term u n is given by
u n t = 0 t t s α 1 2.45 + 10 2 | max t [ 10 1 , 2 ] u n 1 s | + 10 4 | u n 1 0.994 1 + s s 0.004 | d s ,
if the right-hand side belongs to E ϕ , 1.9 . If not, we have
u n t = ( 1 β n ) u n 1 + β n ϕ 0 + 1.132 Γ α 0 t t s α 1 2.45 + 10 2 | max t [ 10 1 , 2 ] u n 1 s | + 10 4 | u n 1 0.994 1 + s s 0.004 | d s ,
with β n 0 , 1 , such that the right-hand side belongs to E ϕ , 1.9 .

5. Conclusions

In this contribution, the investigated question concerns the existence and uniqueness of the solution for a class of nonlinear functional differential equations of fractional order. The considered problems in Equations (1) and (2) are distinguished by the fact that the nonlinear part depends on maximum values of the unknown function, which is not frequently discussed in the existing literature. These maximums are taken on time-dependent intervals and, moreover, the equation is of mixed type, i.e., with both retarded and advanced deviations. It should be noted that, if the hypotheses H 6 and Equation (20) are omitted, the operator F can be a self mapping, and thus, by the usual contraction methods, it can be shown that the result of Theorem 1 remains valid with the bound in Equation (21) weakened to 1. When additional conditions are necessary to meet the physical or mechanical requirements of the phenomenon governed by Equations (1) and (2), we leave the previous usual framework of study. In this case, our main result of Theorem 1 shows that the condition in Equation (21) is sufficient for the existence and uniqueness of the solution.

Acknowledgments

The authors thank the anonymous reviewers for their relevant comments and valuable suggestions, which significantly contributed to improving the quality of the paper. This research was supported by the Ministry of Higher Education and Scientific Research of Algeria, under the program of National Commission for the Evaluation of University Research Projects, Grant Agreement No. B01120140060.

Author Contributions

K. Nisse and L. Nisse contributed equally in this work. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Popov, E.P. Automatic Regulation and Control; Nauka: Moscow, Russia, 1966. [Google Scholar]
  2. Popov, E.P. The Dynamics of Automatic Control Systems; Pergamon Press: Oxford, UK, 1962. [Google Scholar]
  3. Agarwal, R.P.; Hristova, S. Strict stability in terms of tow measures for impulsive differential equations with supremum. Appl. Anal. 2012, 91, 1379–1392. [Google Scholar] [CrossRef]
  4. Bainov, D.D.; Hristova, S.G. Differential Equations with Maxima; Chapman and Hall/CRC: Boca Raton, FL, USA, 2011. [Google Scholar]
  5. Stepanov, E. On solvability of some boundary value problems for differential equations with maxima. Topol. Methods Nonlinear Anal. 1996, 8, 315–326. [Google Scholar] [CrossRef]
  6. Agarwal, R.P.; Hristova, S. Quasilinearization for initial value problems involving differential equations with “maxima”. Math. Comput. Model. 2012, 55, 2096–2105. [Google Scholar] [CrossRef]
  7. Angelov, V.G.; Bainov, D.D. On functional differential equations with Maximums. Period. Math. Hung. 1987, 18, 7–15. [Google Scholar] [CrossRef]
  8. Bohner, M.J.; Georgieva, A.T.; Hristova, S.G. Nonlinear differential equations with maxima: Parametric stability in terms of two measures. Appl. Math. Inf. Sci. 2013, 7, 41–48. [Google Scholar] [CrossRef]
  9. Golev, A.; Hristova, S.; Rahnev, A. An algorithm for approximate solving of differential equations with “maxima”. Comput. Math. Appl. 2010, 60, 2771–2778. [Google Scholar] [CrossRef]
  10. Otrocol, D. Systems of functional differential equations with maxima of mixed type, Volume 191 of Graduate Texts in Mathematics. Electron. J. Qual. Theory Differ. Equ. 2014, 5, 1–9. [Google Scholar]
  11. Tsachev, T.; Angelov, V.G. Fixed points of nonself-mappings and applications. Nonlinear Anal. Theory Methods Appl. 1993, 21, 1–12. [Google Scholar] [CrossRef]
  12. Araci, A.; Sen, E.; Acikgoz, M.; Srivastava, H.M. Existence and uniqueness of positive and nondecreasing solutions for a class of fractional boundary value problems involving the p-Laplacian operator. Adv. Differ. Equ. 2015, 40. [Google Scholar] [CrossRef]
  13. Jankowski, T. Fractional equations of Volterra type involving a Riemann-Liouville derivative. Appl. Math. Lett. 2013, 26, 344–350. [Google Scholar] [CrossRef]
  14. Mishra, L.N.; Srivastava, H.M.; Sen, M. Existence results for some nonlinear functional-integral equations in Banach algebra with applications. Int. J. Anal. Appl. 2016, 11, 1–10. [Google Scholar]
  15. Nisse, L.; Bouaziz, A. Existence and stability of the solutions for systems of nonlinear fractional differential equations with deviating arguments. Adv. Differ. Equ. 2014, 275. [Google Scholar] [CrossRef]
  16. Wang, G.; Agarwal, R.P.; Cabada, A. Existence results and the monotone iterative technique for systems of nonlinear fractional differential equations. Appl. Math. Lett. 2012, 25, 1019–1024. [Google Scholar] [CrossRef]
  17. Ahmad, B.; Ntouyas, S.K. Initial value problems of fractional order Hadamard-type functional differential equations. Electron. J. Differ. Equ. 2015, 77, 1–9. [Google Scholar]
  18. Tisdell, C.C. Basic existence and a priori bound results for solutions to systems of boundary value problems for fractional differential equations. Electron. J. Differ. Equ. 2016, 84, 1–9. [Google Scholar]
  19. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin, Germany, 2004. [Google Scholar]
  20. Gorenflo, R.; Mainardi, F. Fractional calculus: Integral and differential equations of fractional order. In Fractional Calculus in Continuum and Statistical Mechanics (Udine, 1996); CISM Courses and Lectures; Springer: Viena, Austria, 1997; pp. 223–276. [Google Scholar]
  21. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  22. Mainardi, F. Fractional calculus: Some basic problems in Continuum and Statistical Mechanics. In Fractal and Fractional Calculus in Continuum Mechanics (Udine, 1996); CIAM Courses and Lectures; Springer: Wien, Austria; New York, NY, USA, 1997; pp. 291–348. [Google Scholar]
  23. Diethelm, K. Analysis of Fractional Differential Equations. J. Math. Anal. Appl. 2002, 265, 229–248. [Google Scholar] [CrossRef]
  24. Jalilian, Y.; Jalilian, R. Existence of solution for delay fractional differential equations. Mediterr. J. Math. 2013, 10, 1731–1747. [Google Scholar] [CrossRef]
  25. Koethe, G. Topological Vector Spaces I; Springer: Berlin, Germany, 1969. [Google Scholar]
  26. Sehgal, V.M.; Singh, S.P. On a fixed point theorem of Krasnoselskii for locally convex spaces. Pac. J. Math. 1976, 62, 561–567. [Google Scholar] [CrossRef]
  27. Angelov, V.G. A Converse to a Contraction Mapping Theorem in Uniform Spaces. Nonlinear Anal. Theory Methods Appl. 1988, 12, 989–996. [Google Scholar] [CrossRef]
  28. Angelov, V.G. Fixed point theorems in uniform spaces and applications. Czechoslov. Math. J. 1987, 37, 19–33. [Google Scholar]

Share and Cite

MDPI and ACS Style

Nisse, K.; Nisse, L. An Iterative Method for Solving a Class of Fractional Functional Differential Equations with “Maxima”. Mathematics 2018, 6, 2. https://doi.org/10.3390/math6010002

AMA Style

Nisse K, Nisse L. An Iterative Method for Solving a Class of Fractional Functional Differential Equations with “Maxima”. Mathematics. 2018; 6(1):2. https://doi.org/10.3390/math6010002

Chicago/Turabian Style

Nisse, Khadidja, and Lamine Nisse. 2018. "An Iterative Method for Solving a Class of Fractional Functional Differential Equations with “Maxima”" Mathematics 6, no. 1: 2. https://doi.org/10.3390/math6010002

APA Style

Nisse, K., & Nisse, L. (2018). An Iterative Method for Solving a Class of Fractional Functional Differential Equations with “Maxima”. Mathematics, 6(1), 2. https://doi.org/10.3390/math6010002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop