Next Article in Journal
General Relativity with a Positive Cosmological Constant Λ as a Gauge Theory
Next Article in Special Issue
The 3D Navier–Stokes Equations: Invariants, Local and Global Solutions
Previous Article in Journal / Special Issue
Heteroclinic Solutions for Classical and Singular ϕ-Laplacian Non-Autonomous Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Order Coupled Systems With Functional and Periodic Boundary Conditions: Existence Results and Application to an SIRS Model

1
Department of Mathematics, British University of Vietnam, Ecopark Campus, 160000 Hung Yen, Hanoi, Vietnam
2
Centro de Investigação em Matemática e Aplicações (CIMA), Instituto de Investigação e Formação Avançada, Universidade de Évora, Rua Romão Ramalho, 59, 7000-671 Évora, Portugal
3
Departamento de Matemática, Escola de Ciências e Tecnologia, Universidade de Évora, 7000-812 Évora, Portugal
*
Author to whom correspondence should be addressed.
Axioms 2019, 8(1), 23; https://doi.org/10.3390/axioms8010023
Submission received: 6 January 2019 / Revised: 13 February 2019 / Accepted: 13 February 2019 / Published: 16 February 2019
(This article belongs to the Special Issue New Trends in Differential and Difference Equations and Applications)

Abstract

:
The results presented in this paper deal with the existence of solutions of a first order fully coupled system of three equations, and they are split in two parts: 1. Case with coupled functional boundary conditions, and 2. Case with periodic boundary conditions. Functional boundary conditions, which are becoming increasingly popular in the literature, as they generalize most of the classical cases and in addition can be used to tackle global conditions, such as maximum or minimum conditions. The arguments used are based on the Arzèla Ascoli theorem and Schauder’s fixed point theorem. The existence results are directly applied to an epidemic SIRS (Susceptible-Infectious-Recovered-Susceptible) model, with global boundary conditions.

1. Introduction

In this paper two different problems are analyzed.
Part one is concerned with the study of a fully nonlinear coupled system of equations
u 1 ( t ) = f 1 t , u 1 ( t ) , u 2 ( t ) , u 3 ( t ) u 2 ( t ) = f 2 t , u 1 ( t ) , u 2 ( t ) , u 3 ( t ) u 3 ( t ) = f 3 t , u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ,
f i : [ a , b ] × R 3 R and i = 1 , 2 , 3 are L 1 Carathéodory functions, subject to the nonlinear functional boundary conditions
u 1 ( a ) = L 1 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 2 ( a ) , u 3 ( a ) u 2 ( a ) = L 2 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 3 ( a ) u 3 ( a ) = L 3 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 2 ( a )
where L i : C [ a , b ] 3 × R 5 R , i = 1 , 2 , 3 , are continuous functions with properties later to be defined. The technique used for the functional problem is based on the Arzèla-Ascoli theorem and Schauder’s fixed point theorem.
Part two, more precisely, Section 4, deals with the fully nonlinear couple system of Equation (1) coupled with the periodic boundary conditions
u 1 ( a ) = u 1 ( b ) u 2 ( a ) = u 2 ( b ) u 3 ( a ) = u 3 ( b ) .
Given that the conditions on L i , do not allow the problem (1)–(2) to cover the periodic case, a different approach for the problem (1)–(3) is required. In this case, in order to obtain the existence and location of periodic solutions, the upper and lower solutions method, along with some adequate local monotone assumptions on the nonlinearities, is used.
Mathematical modelling and applications are becoming increasingly popular nowadays. With the sudden outburst of keywords such as big data, data analytics and modelling, the quest for mathematical models is on high demand. In the area of mathematical modelling, systems of differential equations are a must, due to their high applicability in areas such as population dynamics [1,2,3], finance [4], medicine [5], biotechnology [6] and physics [7,8], and also examples treated in [9,10].
Nevertheless, in the literature available, the cases dealing with coupled systems of equations are not abundant. Such systems can be found in [11,12,13], however, in this paper the authors present a problem where both equations are coupled. In addition to that, to the best of our knowledge, it is the first time where coupled systems are considered with coupled functional boundary conditions.
This feature allows to generalize the classical boundary data in the literature, such as two-point or multi-point, nonlinear, nonlocal, integro-differential conditions, among others. Indeed, the functional part can deal with global boundary assumptions, such as minimum or maximum arguments, infinite multi-point data, and integral conditions on the several unknown functions. Functional problems, along with their features, can be seen in [14,15,16,17,18,19] and the references therein.
The methods and techniques applied in this paper can be easily adapted to coupled systems with n equations and variables. However, as the notation and writing appear to be heavy, and may avoid the clarity of the results, we prefer to prove our theoretical part for n = 3 , which is adequate for our application.
This paper has the following structure: Section 2 contains some definitions and generic assumptions on the nonlinearities. Section 3 shows the main result for problems with functional boundary conditions. In Section 4 it is studied the periodic problem via lower and upper solutions technique together with some local growth conditions. The final section presents an application of (1)–(2) to an epidemic SIRS model to illustrate the applicability of the problem discussed and to show the potentialities of the functional boundary conditions, exploring global initial boundary conditions on the system.

2. Definitions and Assumptions

Throughout this work we consider the space of continuous functions in [ a , b ] , on the Banach space E : = C [ a , b ] 3 , equipped with the norm
u E = max u i , i = 1 , 2 , 3
where u i = max t a , b u i t .
The functional boundary functions verify the assumption:
Hypothesis 1.
L i : C [ a , b ] 3 × R 5 R are continuous functions. Moreover, each of the functions, L i η , x , y are uniformly continuous when η , x , y is bounded.
The admissible nonlinearities will be L 1 Carathéodory functions, according the following definition:
Definition 1.
The functions f i : [ a , b ] × R 3 R , i = 1 , , n , are L 1 Carathéodory if they verify:
(i
for each y 1 , y 2 , y 3 R 3 , t f i ( t , y 1 , y 2 , y 3 ) are measurable on [ a , b ] , for i = 1 , 2 , 3 ;
(ii
for almost every t [ a , b ] , y 1 , y 2 , y 3 f i ( t , y 1 , y 2 , y 3 ) are continuous on R 3 , for i = 1 , 2 , 3 ;
(iii
for each L > 0 , there exists a positive function ψ i L L 1 a , b , i = 1 , 2 , 3 , such that, for max y i , i = 1 , 2 , 3 < L ,
f i ( t , y 1 ( t ) , y 2 ( t ) , y 3 ( t ) ) ψ i L ( t ) , a . e . t [ a , b ] , i = 1 , 2 , 3 .
To demonstrate the final result, Schauder’s fixed point theorem will be an important tool to guarantee the existence of fixed points for the operator to be defined:
Theorem 1.
([20]) Let Y be a nonempty, closed, bounded and convex subset of a Banach space X, and suppose that P : Y Y is a completely continuous operator. Then P has at least one fixed point in Y.

3. Main Result for Functional Problems

In this section, we present and prove the main existence result for (1)–(2), given by the following theorem:
Theorem 2.
If f i are L 1 Carathéodory functions, for i = 1 , 2 , 3 , and the continuous functions L i , i = 1 , 2 , 3 , verify (H1) and
Hypothesis 2.
there exists R > 0 such that
max k i + a t ψ i R s d s , i = 1 , 2 , 3 R ,
with
k 1 : = max L 1 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 2 ( a ) , u 3 ( a ) , k 2 : = max L 2 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 3 ( a ) , k 3 : = max L 3 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 2 ( a ) ,
then the problem (1)–(2) has at least one solution u = ( u 1 , u 2 , u 3 ) C a , b 3 .
Proof. 
Let us consider the integral system given by
u 1 ( t ) = L 1 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 2 ( a ) , u 3 ( a ) + a t f 1 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s , u 2 ( t ) = L 2 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 3 ( a ) + a t f 2 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s , u 3 ( t ) = L 3 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 2 ( a ) + a t f 3 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s ,
and the operator
T : C a , b 3 C a , b 3
defined by
T u 1 , u 2 , u 3 = T 1 u 1 , u 2 , u 3 , T 2 u 1 , u 2 , u 3 , T 3 u 1 , u 2 , u 3 ,
where T i : C a , b 3 C a , b , i = 1 , 2 , 3 , given by
T 1 u 1 , u 2 , u 3 = L 1 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 2 ( a ) , u 3 ( a ) + a t f 1 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s , T 2 u 1 , u 2 , u 3 = L 2 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 3 ( a ) + a t f 2 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s , T 3 u 1 , u 2 , u 3 = L 3 u 1 , u 2 , u 3 , u 1 ( b ) , u 2 ( b ) , u 3 ( b ) , u 1 ( a ) , u 2 ( a ) + a t f 3 s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s ,
As the fixed points of T are fixed points of T i , i = 1 , 2 , 3 , and vice-versa, which are solutions of (1)–(2), the goal of this architecture will be to use the Arzèla-Ascoli theorem and Schauder’s fixed point theorem to prove that the problem (1)–(2) has at least one solution.
For clarity, we consider several claims:
Claim 1.
T D D , for some D C a , b 3 a bounded, closed and convex subset.
Consider
D = u 1 , u 2 , u 3 C a , b 3 : u 1 , u 2 , u 3 E k ,
with k > 0 such that
max k i + a t ψ i k s d s , i = 1 , 2 , 3 k .
Given u 1 , u 2 , u 3 D , by Definition 1 (iii) and (Hypothesis 2), we have that for u 1 , u 2 , u 3 D , i = 1 , 2 , 3 ,
f i s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) ψ i k s , a . e . s a , b .
Therefore, T u 1 , u 2 , u 3 E k and Claim 1 is proved.
Claim 2.
The operator T is completely continuous.
To prove that the operator T is completely continuous it is sufficient to show that T is uniformly bounded and T is equicontinuous. Using the above arguments, it can be proved that T i are uniformly bounded, for i = 1 , 2 , 3 , and therefore T is uniformly bounded.
In order to show that the operator T is equicontinuous, let us consider t 1 , t 2 a , b , such that, without any loss of generality, t 1 < t 2 .
Then for T i , i = 1 , 2 , 3 , we have
T i u 1 t 2 , u 2 t 2 , u 3 t 2 T i u 1 t 1 , u 2 t 1 , u 3 t 1 t 1 t 2 f i s , u 1 ( s ) , u 2 ( s ) , u 3 ( s ) d s t 1 t 2 ψ i L s d s 0 ,
as t 1 t 2 . So, each operator T i is equicontinuous and, hence, the operator T is equicontinuous.
Therefore, by Arzèla-Ascoli’s theorem, the operator T is compact and using Schauder’s fixed point theorem, we obtain that T has a fixed point, that is, the problem (1)–(2) has at least a solution u C a , b 3 .  □
Remark that (Hypothesis 2) implies that the periodic case is not covered by (2).

4. Existence and Localization Result for the Periodic Case

Consider now the system (1) with the periodic boundary conditions (3). As it is well known, in nonlinear differential equations, the periodic case is more delicate and requires a different approach for general nonlinearities.
The method to be used will apply lower and upper solutions technique, based on the definition:
Definition 2.
Consider the C 1 -functions α i , β i : a , b R , i = 1 , 2 , 3 . The triple α 1 , α 2 , α 3 is a lower solution of the periodic problem (1), (3) if
α 1 ( t ) f 1 t , α 1 ( t ) , α 2 ( t ) , α 3 ( t ) α 2 ( t ) f 2 t , α 1 ( t ) , α 2 ( t ) , α 3 ( t ) α 3 ( t ) f 3 t , α 1 ( t ) , α 2 ( t ) , α 3 ( t )
and
α 1 ( a ) α 1 ( b ) α 2 ( a ) α 2 ( b ) α 3 ( a ) α 3 ( b ) .
The triple β 1 , β 2 , β 3 is an upper solution of the periodic problem (1), (3) if the reversed inequalities hold.
This method allows to obtain an existence and localization theorem:
Theorem 3.
Let α 1 , α 2 , α 3 and β 1 , β 2 , β 3 be lower and upper solutions of (1), (3), respectively, such that α i ( t ) β i ( t ) , t a , b and for i = 1 , 2 , 3 .
Define the set
A = t , u 1 , u 2 , u 3 a , b × R 3 : α i ( t ) u i β i ( t ) , i = 1 , 2 , 3
and assume that f i are L 1 Carathéodory functions on A, for i = 1 , 2 , 3 verifying
f 1 t , x , α 2 ( t ) , α 3 ( t ) f 1 t , x , y , z f 1 t , x , β 2 ( t ) , β 3 ( t ) ,
for t a , b , α 2 ( t ) y β 2 ( t ) , α 3 ( t ) z β 3 ( t ) ,
f 2 t , α 1 ( t ) , y , α 3 ( t ) f 2 t , x , y , z f 2 t , β 1 ( t ) , y , β 3 ( t ) ,
for t a , b , α 1 ( t ) x β 1 ( t ) , α 3 ( t ) z β 3 ( t ) ,
f 3 t , α 1 ( t ) , α 2 ( t ) , z f 3 t , x , y , z f 3 t , β 1 ( t ) , β 2 ( t ) , z ,
for t a , b , α 1 ( t ) x β 1 ( t ) , α 2 ( t ) y β 2 ( t ) .
Then the problem (1), (3) has, at least, a solution u = ( u 1 , u 2 , u 3 ) C a , b 3 such that
α i ( t ) u i ( t ) β i ( t ) , i = 1 , 2 , 3 , for all t a , b .
Proof. 
For i = 1 , 2 , 3 , define the truncature functions δ i given by
δ i ( t , u i ) = α i ( t ) if u i < α i ( t ) u i if α i ( t ) u i β i ( t ) β i ( t ) if u i > β i ( t )
and consider the modified problem composed by the truncated and perturbed differential equations
u 1 ( t ) + u 1 ( t ) = f 1 t , δ 1 ( t , u 1 ) , δ 2 ( t , u 2 ) , δ 3 ( t , u 3 ) + δ 1 ( t , u 1 ) u 2 ( t ) + u 2 ( t ) = f 2 t , δ 1 ( t , u 1 ) , δ 2 ( t , u 2 ) , δ 3 ( t , u 3 ) + δ 2 ( t , u 2 ) u 3 ( t ) + u 3 ( t ) = f 3 t , δ 1 ( t , u 1 ) , δ 2 ( t , u 2 ) , δ 3 ( t , u 3 ) + δ 3 ( t , u 3 ) ,
together with the boundary conditions (3).
As the linear and homogeneous problem associated to (10), (3) has only the null solution, then we can write (10), (3) in the integral form
u 1 ( t ) = a b G 1 ( t , s ) f 1 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 1 ( s , u 1 ( s ) ) d s u 2 ( t ) = a b G 2 ( t , s ) f 2 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 2 ( s , u 2 ( s ) ) d s u 3 ( t ) = a b G 3 ( t , s ) f 3 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 3 ( s , u 3 ( s ) ) d s ,
where G i ( t , s ) are the Green functions corresponding to the problem
u i ( t ) + u i ( t ) = f i t u i ( a ) = u i ( b ) ,
for i = 1 , 2 , 3 .
Then the operator
T : C a , b 3 C a , b 3
given by
T u 1 , u 2 , u 3 = T 1 u 1 , u 2 , u 3 , T 2 u 1 , u 2 , u 3 , T 3 u 1 , u 2 , u 3 ,
with T i : C a , b 3 C a , b , i = 1 , 2 , 3 , defined as
T 1 u 1 , u 2 , u 3 = a b G 1 ( t , s ) f 1 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 1 ( s , u 1 ( s ) ) d s T 2 u 1 , u 2 , u 3 = a b G 2 ( t , s ) f 2 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 2 ( s , u 2 ( s ) ) d s T 3 u 1 , u 2 , u 3 = a b G 3 ( t , s ) f 3 s , δ 1 ( s , u 1 ( s ) ) , δ 2 ( s , u 2 ( s ) ) , δ 3 ( s , u 3 ( s ) ) + δ 3 ( s , u 3 ( s ) ) d s
is completely continuous in C a , b 3 .
By Schauder’s fixed point theorem, we obtain that T has a fixed point, that is, the problem (10), (3) has at least a solution u : = u 1 , u 2 , u 3 C a , b 3 .
To prove that this function u is a solution of the initial problem (1), (3), it will be enough to show that
α i ( t ) u i ( t ) β i ( t ) , i = 1 , 2 , 3 , for all t a , b .
Suppose that the first inequality does not hold for i = 1 and some t 1 a , b , that is
u 1 ( t 1 ) < α 1 ( t 1 ) .
Extend the function u 1 α 1 by periodicity and consider
t 0 : = inf t t 1 b + a , t 1 : s ] t , t 1 ] , u 1 ( s ) < α 1 ( s ) .
Therefore, for t [ t 0 , t 1 ] we have, by Definition 2 and (9),
u 1 ( t ) α 1 ( t ) f 1 t , α 1 ( t ) , δ 2 ( t , u 2 ) , δ 3 ( t , u 3 ) + α 1 ( t ) u 1 ( t ) f 1 t , α 1 ( t ) , α 2 ( t ) , α 3 ( t ) α 1 ( t ) u 1 ( t ) > 0 .
So, α 1 u 1 is increasing on t 0 , t 1 and
0 > u 1 ( t 1 ) α 1 ( t 1 ) > u 1 ( t 0 ) α 1 ( t 0 ) ,
where t 0 = t 1 b + a , which is in contradiction of the periodicity of the extension of u 1 α 1 . So α 1 ( t ) u 1 ( t ) , for all t a , b .
Applying similar arguments it can be proved that u 1 ( t ) β 1 ( t ) , for t a , b , and
α i ( t ) u i ( t ) β i ( t ) , i = 2 , 3 , for all t a , b .
 □

5. An Epidemic Model of an SIRS System With Nonlinear Incidence Rate and Interaction from Infectious to Susceptible Subjects

The existent literature has innumerous examples of applications of SIR models, namely in [21,22,23] where the population is divided into Susceptible S , Infectious I and Recovery R . However, in SIR models, recovered individuals are assumed to develop lifelong immunity, which for some diseases such as seasonal flu, influenza or venereal diseases is not necessarily true. A recovered individual becomes susceptible and possibly infected after some time. In this case, SIRS models, where recovered individuals lose immunity and become susceptible again, are far more adequate. Examples can be found in [24,25,26,27,28,29]. In [30], the authors develop and explore a mathematical model of an SIRS epidemic model where a transfer between infectious and susceptible rate is included, as shown in Figure 1.
This approach was designed to model the cases where recovery cannot generate immunity for a long time. Infected individuals may recover after some treatments and go back directly to the susceptible class. In addition, a nonlinear incidence with the average number of new cases of a disease per unit time, S g ( I ) , is included, as suggested in [31]. This nonlinear incidence replicates a more complex dynamic than the ones presented on bilinear or standard incidence models, and it seems to perform better when modelling more complex cholera cases, as shown in [31].
The model presented is then
d S d t = 1 100 Λ μ S S g ( I ) + γ 1 I + δ R , d I d t = 1 100 S g ( I ) μ + γ 1 + γ 2 + α I , d R d t = 1 100 γ 2 I μ + δ R
where:
  • Λ represents the recruitment rate of susceptible individuals;
  • μ is the natural death rate;
  • γ 1 is the transfer rate from the infected class to the susceptible class;
  • γ 2 is the transfer rate from the infected class to the recovered class;
  • α is the disease-induced death rate;
  • δ the immunity loss rate.
In this model Λ and μ are assumed to be positive and δ , γ 1 , γ 2 and α are assumed to be nonnegative. As per the nonlinear incidence, S g ( I ) , g is a real locally Lipschitz function on R + with the following conditions, as presented in [32]:
  • g 0 = 0 and g I > 0 for I > 0 ,
  • g I / I is continuous and monotonously increasing for I > 0 and lim I 0 + g ( I ) I exists as β > 0 .
Unlike the previous models presented, in this paper we couple the system (12) with the boundary conditions, for t 0 , T and T > 0 ,
S ( 0 ) = max t 0 , T I ( t ) , I 0 = 0 T S ( T ) S ( s ) d s , R 0 = 0 T S ( T ) I ( s ) d s .
These functional boundary conditions fully reflect the considerations and conclusions presented in [30]. These types of boundary conditions have the following practical meaning:
  • S ( 0 ) , the initial number of susceptible subjects, is equal to the maximum of the infected subjects;
  • I 0 , the initial number of infected subjects, is a weighted average of the susceptible individuals, weighted by the final value of the susceptible S, at time T ;
  • R 0 , the initial number of individuals who recovered, is equal to a weighted average of the infected individuals, weighted by the final value of the susceptible, S, at time T.
Functional boundary conditions allow these assumptions to be examined in full, as operators can be considered as boundary conditions. These types of global conditions, that can include somewhat more abstract conditions, can only be contemplated via functional boundary conditions.
This model can therefore be presented in the form of (1)–(2), where u 1 = S , u 2 = I , u 3 = R , a = 0 and b = T ,
f 1 t , u 1 , u 2 , u 3 = Λ μ u 1 S g ( u 2 ) + γ 1 u 2 + δ u 3 , f 2 t , u 1 , u 2 , u 3 = u 1 g ( I ) μ + γ 1 + γ 2 + α u 2 , f 3 t , u 1 , u 2 , u 3 = γ 2 u 1 μ + δ u 3 ,
L 1 u 1 , u 2 , u 3 , u 1 T , u 2 T , u 3 T , u 2 0 , u 3 0 = max t 0 , T u 2 : = k 1 , L 2 u 1 , u 2 , u 3 , u 1 T , u 2 T , u 3 T , u 1 0 , u 3 0 = 0 T u 1 ( T ) u 1 ( s ) d s : = k 2 , L 3 u 1 , u 2 , u 3 , u 1 T , u 2 T , u 3 T , u 1 0 , u 2 0 = 0 T u 1 ( T ) u 2 ( s ) d s : = k 3 .
It is clear that the boundary conditions L 1 , L 2 and L 3 satisfy (H1).
Moreover, f 1 , f 2 and f 3 are L 1 Carathéodory functions such that, for max u 1 , u 2 , u 3 < k ,
f 1 t , u 1 , u 2 , u 3 Λ + μ k + k β + γ 1 k + δ k : = ψ 1 k ( t ) , f 2 t , u 1 , u 2 , u 3 k β + μ + γ 1 + γ 2 + α k : = ψ 2 k ( t ) , f 3 t , u 1 , u 2 , u 3 γ 2 k + μ + δ k : = ψ 3 k ( t ) .
To satisfy (Hypothesis 2), one must verify that there exists k > 0 , such that,
k max k 1 + 0 T Λ + k + μ k + β + γ 1 + δ d s , k 2 + 0 T k β + μ + γ 1 + γ 2 + α d s , k 3 + 0 T k γ 2 + μ + δ d s , .
Let T = 1 , Λ = 2 , μ = 12.6 , γ 1 = 2.7 , γ 2 = 0.28 , α = 6.32 , δ = 0.38 , k 1 = 2 , k 2 = 0.27 , k 3 = 5 , β = 12 .
For k 0.06 , condition (Hypothesis 2) is verified and therefore, by Theorem 2, the system (12)–(13) has at least one solution S , I , R C 0 , 1 3 , for the values considered.

6. Conclusions

In this paper, the authors show the existence of solution for a first order fully coupled system of three equations, involving two different cases. The first case, with coupled functional boundary conditions, is an existence result. The second case, with periodic boundary conditions, which is not covered by the first result, is an existence and location result. The extra information obtained in this second case is associated with the technique used, as it relies on the upper and lower solution method.
The application to an SIRS model, with global boundary conditions, underlines the key advantage and flexibility of the functional boundary conditions. The example shown illustrates not only the theorem proved in this paper, but it also provides guidance on how to arrange global conditions, in order to conform with the layout in (2).
As a matter of fact, a similar approach can be taken in several other models, allowing global conditions to be considered as boundary conditions, highly increasing the level of applicability of these models.

Author Contributions

The authors equally contributed to this work.

Funding

This project was supported by Fundação para a Ciência e a Tecnologia (FCT) via project UID/MAT/04674/2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Al-Moqbali, M.K.A.; Al-Salti, N.S.; Elmojtaba, I.M. Prey–Predator Models with Variable Carrying Capacity. Mathematics 2018, 6, 102. [Google Scholar] [CrossRef]
  2. Gumus, O.A.; Kose, H. On the Stability of Delay Population Dynamics Related with Allee Effects. Math. Comput. Appl. 2012, 17, 56–67. [Google Scholar] [CrossRef] [Green Version]
  3. Song, H.S.; Cannon, W.R.; Beliaev, A.S.; Konopka, A. Mathematical Modeling of Microbial Community Dynamics: A Methodological Review. Processes 2014, 2, 711–752. [Google Scholar] [CrossRef] [Green Version]
  4. Osakwe, C.J.U. Incentive Compatible Decision Making: Real Options with Adverse Incentives. Axioms 2018, 7, 9. [Google Scholar] [CrossRef]
  5. Deuflhard, P. Differential equations in technology and medicine: Computational concepts, adaptive algorithms, and virtual labs. In Computational Mathematics Driven by Industrial Problems; Springer: Berlin/Heidelberg, Germany, 2000; pp. 69–125. [Google Scholar]
  6. Jang, S.S.; de la Hoz, H.; Ben-zvi, A.; McCaffrey, W.C.; Gopaluni, R.B. Parameter estimation in models with hidden variables: An application to a biotech process. Can. J. Chem. Eng. 2012, 90, 690–702. [Google Scholar] [CrossRef]
  7. Akarsu, M.; Özbaş, Ö. Monte Carlo Simulation for Electron Dynamics in Semiconductor Devices. Math. Comput. Appl. 2005, 10, 19–26. [Google Scholar] [CrossRef] [Green Version]
  8. Malinzi, J.; Quaye, P.A. Exact Solutions of Non-Linear Evolution Models in Physics and Biosciences Using the Hyperbolic Tangent Method. Math. Comput. Appl. 2018, 23, 35. [Google Scholar] [CrossRef]
  9. Nieto, J.J. Periodic boundary value problems for first-order impulsive ordinary differential equations. Nonlinear Anal. 2002, 51, 1223–1232. [Google Scholar] [CrossRef]
  10. Zhang, W.; Fan, M. Periodicity in a generalized ecological competition system governed by impulsive differential equations with delays. Math. Comput. Model. 2004, 39, 479–493. [Google Scholar] [CrossRef]
  11. Agarwal, R.P.; O’Reagan, D. A coupled system of boundary value problems. Appl. Anal. 1998, 69, 381–385. [Google Scholar] [CrossRef]
  12. Asif, N.A.; Talib, I.; Tunc, C. Existence of solutions for first-order coupled system with nonlinear coupled boundary conditions. Bound. Val. Prob. 2015, 2015, 134. [Google Scholar] [CrossRef]
  13. Asif, N.A.; Khan, R.A. Positive solutions to singular system with four-point coupled boundary conditions. J. Math. Anal. Appl. 2012, 386, 848–861. [Google Scholar] [CrossRef] [Green Version]
  14. Cabada, A.; Fialho, J.; Minhós, F. Extremal solutions to fourth order discontinuous functional boundary value problems. Math. Nachr. 2013, 286, 1744–1751. [Google Scholar] [CrossRef] [Green Version]
  15. Cabada, A.; Pouso, R.; Minhós, F. Extremal solutions to fourth-order functional boundary value problems including multipoint condition. Nonlinear Anal. Real World Appl. 2009, 10, 2157–2170. [Google Scholar] [CrossRef]
  16. Fialho, J.; Minhós, F. Higher order functional boundary value problems without monotone assumptions. Bound. Val. Prob. 2013, 2013, 81. [Google Scholar] [CrossRef] [Green Version]
  17. Fialho, J.; Minhós, F. Multiplicity and location results for second order functional boundary value problems. Dyn. Syst. Appl. 2014, 23, 453–464. [Google Scholar]
  18. Graef, J.; Kong, L.; Minhós, F. Higher order boundary value problems with ϕ -Laplacian and functional boundary conditions. Comput. Math. Appl. 2011, 61, 236–249. [Google Scholar] [CrossRef]
  19. Graef, J.; Kong, L.; Minhós, F.; Fialho, J. On the lower and upper solution method for higher order functional boundary value problems. Appl. Anal. Discret. Math. 2011, 5, 133–146. [Google Scholar] [CrossRef]
  20. Zeidler, E. Nonlinear Functional Analysis and Its Applications, I: Fixed-Point Theorems; Springer: New York, NY, USA, 1986. [Google Scholar]
  21. Angstmann, C.N.; Henry, B.I.; McGann, A.V. A Fractional-Order Infectivity and Recovery SIR Model. Fract. Fract. 2017, 1, 11. [Google Scholar] [CrossRef]
  22. Cui, Q.; Qiu, Z.; Liu, W.; Hu, Z. Complex Dynamics of an SIR Epidemic Model with Nonlinear Saturate Incidence and Recovery Rate. Entropy 2017, 19, 305. [Google Scholar] [CrossRef]
  23. Secer, A.; Ozdemir, N.; Bayram, M. A Hermite Polynomial Approach for Solving the SIR Model of Epidemics. Mathematics 2018, 6, 305. [Google Scholar] [CrossRef]
  24. Alexander, M.E.; Moghadas, S.M. Bifurcation analysis of an SIRS epidemic model with generalized incidence. SIAM J. Appl. Math. 2005, 65, 1794–1816. [Google Scholar] [CrossRef]
  25. Chen, J. An SIRS epidemic model. Appl. Math. J. Chin. Univ. 2004, 19, 101–108. [Google Scholar] [CrossRef]
  26. Hu, Z.; Bi, P.; Ma, W.; Ruan, S. Bifurcations of an SIRS epidemic model with nonlinear incidence rate. Discret. Contin. Dyn. Syst. Ser. B 2011, 15, 93–112. [Google Scholar] [CrossRef]
  27. Liu, J.; Zhou, Y. Global stability of an SIRS epidemic model with transport-related infection. Chaos Solitons Fract. 2009, 40, 145–158. [Google Scholar] [CrossRef]
  28. Teng, Z.; Liu, Y.; Zhang, L. Persistence and extinction of disease in non-autonomous SIRS epidemic models with disease-induced mortality. Nonlinear Anal. Theory Methods Appl. 2008, 69, 2599–2614. [Google Scholar] [CrossRef]
  29. Jin, Y.; Wang, W.; Xiao, S. An SIRS model with a nonlinear incidence rate. Chaos Solitons Fract. 2007, 34, 1482–1497. [Google Scholar] [CrossRef]
  30. Li, T.; Zhang, F.; Liu, H.; Chen, Y. Threshold dynamics of an SIRS model with nonlinear incidence rate and transfer from infectious to susceptible. Appl. Math. Lett. 2017, 70, 52–57. [Google Scholar] [CrossRef]
  31. Capasso, V.; Serio, G. A generalization of the Kermack-McKendrick deterministic epidemic model. Math. Biosci. 1978, 42, 43–61. [Google Scholar] [CrossRef]
  32. Li, J.; Yang, Y.; Xiao, Y.; Liu, S. A class of Lyapunov functions and the global stability of some epidemic models with nonlinear incidence. J. Appl. Anal. Comput. 2016, 6, 38–46. [Google Scholar] [CrossRef]
Figure 1. SIRS model diagram.
Figure 1. SIRS model diagram.
Axioms 08 00023 g001

Share and Cite

MDPI and ACS Style

Fialho, J.; Minhós, F. First Order Coupled Systems With Functional and Periodic Boundary Conditions: Existence Results and Application to an SIRS Model. Axioms 2019, 8, 23. https://doi.org/10.3390/axioms8010023

AMA Style

Fialho J, Minhós F. First Order Coupled Systems With Functional and Periodic Boundary Conditions: Existence Results and Application to an SIRS Model. Axioms. 2019; 8(1):23. https://doi.org/10.3390/axioms8010023

Chicago/Turabian Style

Fialho, João, and Feliz Minhós. 2019. "First Order Coupled Systems With Functional and Periodic Boundary Conditions: Existence Results and Application to an SIRS Model" Axioms 8, no. 1: 23. https://doi.org/10.3390/axioms8010023

APA Style

Fialho, J., & Minhós, F. (2019). First Order Coupled Systems With Functional and Periodic Boundary Conditions: Existence Results and Application to an SIRS Model. Axioms, 8(1), 23. https://doi.org/10.3390/axioms8010023

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