![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
This paper investigates a kind of stochastic neutral integro-differential equations with infinite delay and Poisson jumps in the concrete-fading memory-phase space . We suppose that the linear part has a resolvent operator and the nonlinear terms are globally Lipschitzian. We introduce sufficient conditions that ensure the existence and uniqueness of mild solutions by using successive approximation. Moreover, we target exponential stability, including moment exponential stability in
-th (
) and almost surely exponential stability of solutions and their maps. An example illustrates the potential of the main result.
1. Introduction
Research on stochastic differential equations with delay has received attention over the last few decades because of their appropriateness to describe physical systems subject to delays, such as the ones found in biology, medicine, epidemiology, chemistry, physics, and economics (see Helge et al., Citation2010; Intissar., Citation2020; Kostikov & Romanenkov, Citation2020; Mao, Citation2007; Trung, Citation2020 for a brief overview). The qualitative and quantitative properties of solutions of stochastic differential equations with delay, such as the existence, uniqueness, controllability, and stability, have been considered by several authors (see Bouzahir et al., Citation2017; Dieye et al., Citation2017; Diop et al., Citation2014; Taniguchi et al., Citation2002; Zouine et al., Citation2020). One point, in particular, has received a lot of attention: the study of existence and asymptotic behavior of mild solutions of some stochastic differential equations on Hilbert spaces, such as the semigroup approach (Taniguchi et al., Citation2002), comparison theorem (Govindan, Citation2003), Razumikhin-type theorem (Kai & Yufeng, Citation2006), analytic technique (Taniguchi, Citation1998), and Banach fixed-point principle (Diop et al., Citation2014).
The literature shows many dynamical systems modeled by neutral stochastic partial differential Equationequations 2(2)
(2) (Chen et al., Citation2014; Cui et al., Citation2011) [Equation9]
(9)
(9) . For these equations, some contain the derivatives of delayed states, which differ from stochastic partial differential equations with delays that depend on the present and past states only (for more details on this theory and its applications, see Mao et al., Citation2017; Yue, Citation2014). Stochastic integro-differential equations have been intensively studied, with special attention paid to qualitative properties, such as stability, regularity, periodicity, control problems, and optimality conditions (see Dieye et al., Citation2019; Diop et al., Citation2014). Due to the existence of an integral term in the equations, we use here the theory of the resolvent operator instead of the strongly continuous semigroups operator (see Grimmer, Citation1982 for further details).
Yet, most of the researchers dealing with exponential stability have limited their research to finite delay (see Dieye et al., Citation2017; Diop et al., Citation2014). Regarding the infinite delay, most investigations have been done for the case of continuous dependence of solutions on the initial value, considering exponential and asymptotic estimates (see, for instance, the papers Cui & Yan, Citation2012; Mao et al., Citation2017; Ren & Xia, Citation2009; Yue, Citation2014 for an account on phase spaces). It is noteworthy that few contributions exist for characterizing the exponential stability of stochastic equations with infinite delays (see Jiang et al., Citation2016; Wu et al., Citation2017). Jiang et al. (Citation2016) showed the exponential stability for a class of second-order neutral stochastic partial differential equations with infinite delays and impulses using the integral inequality technique. Wu et al. (Citation2017) showed boundedness in the mean square and convergence for both solutions and their maps in the phase space by using the Itô formula.
Motivated by the above discussion, consider the following neutral stochastic integro-differential equations with infinite delay and Poisson jumps given on the complete probability space :
System (1) holds with , where
is
measurable, and the definition of the concrete fading memory-phase space
is detailed in the next section.
is a closed linear operator on a separable Hilbert space
and
represents a closed linear operator such that
has
as a domain, where
, for all
.
,
,
are appropriate functions, the history
, such that
belongs to the phase space
. The process
represents a Wiener process on a separable Hilbert space
and
is a compensated Poisson random measure.
To the best of the authors’ knowledge, this paper is the first to present a study of the existence and exponential stability of neutral stochastic integro-differential equations with infinite delay and Poisson jumps. The main contribution of this paper is to find conditions to ensure existence, uniqueness, exponential stability in -th moment for
and almost surely exponential stability of solutions and their maps of (1). We show the result using stochastic techniques and the resolvent operator theory, as defined in Grimmer (Citation1982). It is worth mentioning that Diop et al. (Citation2014) studied the system (1) with finite delay. They focused only on the existence of mild solutions and their exponential stability in mean square (Diop et al., Citation2014). For this reason, our approach can be seen as an extension of the result of Diop et al. (Citation2014) for the infinite delay case.
The organization of this paper is as follows. Some notations and preliminary results are presented in Section 2. The existence and uniqueness of mild solutions for neutral stochastic integro-differential equations with infinite delay are shown in Section 3. Conditions assuring moment exponential stability in the -th (
) and almost surely exponential stability of the solution
, and the solution maps
,
are shown in Section 4. Finally, an example that illustrates our results is presented in Section 5.
2. Notations and preliminary results
Let and
be two real separable Hilbert spaces (c.f. Taniguchi et al., Citation2002). We let also
be the space of bounded linear operators from
into
associated with
to represent the norm operator in
,
, and
. We assume that System (1) is equipped with a normal filtration
.
Denote by , the Poisson random measure induced by the
finite stationary
adapted Poisson point process
taking values in a measurable space
, and define the compensated Poisson random measure
as
, where
for
and
is the characteristic measure of
.
Let represent the space of all continuous functions from
into
equipped with the norm defined by
. For a given
, consider the fading memory
defined as follows:
as the chosen phase space in this paper. It is a Banach space with the norm
. To know more details, see Hino et al. (Citation1991) and Appendix.
Supposed that represents a
-valued Wiener process which is independent of the Poisson point process on the probability space
with a positive self-adjoint covariance operator
. In addition, we suppose that there exists a complete orthonormal system
in
, a bounded sequence of positive real numbers
such that
, and a sequence
of independent standard Brownian motions such that
for
and
is the
-algebra generated by
(see (Taniguchi et al., Citation2002)). We consider the subspace
of
, it is a Hilbert space equipped with the inner product
. Let
be the space of all Hilbert-Schmidt operators from
to
.
is a separable Hilbert space endowed with the norm
for any
.
Hereafter, and
are closed linear operators on a Banach space denoted by
, and
is the Banach space
endowed with the graph norm
for
.
The notations ,
and
represent the space of continuous functions from
into
, the space of continuously differentiable functions from
into
and the set of bounded linear operators from
into
, respectively.
2.1. Preliminaries on partial integro-differential equations
We now consider the problem
with .
Definition 2.1. (Grimmer, Citation1982). A bounded linear operator-valued function ,
is called a resolvent operator for (2) if the next two conditions are satisfied.
(i) and there exist two constants
and
such that
for all
.
(ii) For each element in
, the function
is strongly continuous for each
and for
in
,
and satisfies
Remark 1. The resolvent operator is said to be exponentially stable when Definition 2.1(i) holds with .
The following two conditions, borrowed from Grimmer (Citation1982), are sufficient to assure the existence of solutions for Equationequation (2)(2)
(2) .
(A1) The operator is an infinitesimal generator of a
-semigroup on
.
(A2) For all ,
denotes a closed, continuous linear operator from
to
and
belongs to
. For any
, the map
is bounded, differentiable, and its derivative
is bounded and uniformly continuous on
.
Lemma 2.2. (Grimmer, Citation1982). Under Assumptions and
, the existence of a resolvent operator for (2) is guaranteed, and it is unique.
We now recall conditions that assure existence of solutions for the deterministic, integro-differential equation
with and
is a continuous function.
Lemma 2.3. (Grimmer, Citation1982). Suppose that Assumptions and
hold, if
is a strict solution of (3) (i.e.,
satisfies (3), for
and it belongs to
), then
Let us give the concept of solutions for the stochastic system in the next definition.
Definition 2.4. A mild solution of (1) is an -valued process
(with
) which satisfies the next two conditions.
(i) is
-adapted and
almost surely.
(ii) is continuous for any
and equals
with , where
in (5) represents the resolvent operator of (2).
To achieve our main goal, we assume the following three assumptions:
(A3) The resolvent operator satisfying Lemma 2.3 is exponentially stable, i.e., there exist two constants
and
such that
, for all
.
(A4) Let be an integer, there exists a real number
, such that
for all , and all
,
(A5) Let be an integer, there exists a real number
, such that
for all and all
,
In particular, there holds .
Remark 2. We consider the assumption for all
, to guarantee that there exists a zero equilibrium solution to the stochastic Equationequation (1)
(1)
(1) . If this assumption does not hold, the equilibrium solution for Equationequation (1)
(1)
(1) can always be transformed into the zero equilibrium solution of another equation.
3. Existence and uniqueness
In this section, we present sufficient conditions to guarantee the existence and uniqueness of mild solutions of the equation in (1). To do so, we use the method of successive approximations and some stochastic analysis techniques. Still, we have to develop some new techniques to deal with infinite delay. Hereafter, we replace by the Hilbert space
in
and
. Now, we present the following main result.
Theorem 3.1. Assumed that Assumptions ,
,
and
hold, and
. Then, Equationequation (1)
(1)
(1) has a unique mild solution.
Proof. The proof of this theorem uses the following sequence of successive approximations that is defined for by
for any
and for
by
for any and
when
. Take
, from the uniform boundedness
. The remaining arguments are divided into three main steps.
Step 1: First, let us show that is a bounded sequence. From (6), for any
, we have
. For any
, we obtain that
.
If , consider the five terms
Considering the definition of the sequence , together with the elementary inequality
, we have
From Assumption , we obtain
and
In addition, combining Assumption and Holder inequality yields
Now, using Lemma 6.2 (see Appendix), and by a similar reasoning on , we show that
On the other hand, combining Assumption , Lemma 6.3 (see Appendix) and Holder inequality on
, we find that
Substituting (8)–(12) into (7), we obtain
where By using the definition of the norm
(see Appendix), we can write
Define
We can write
Therefore, for any , one can conclude that
where . Besides, we have
It follows that
where The Gronwall inequality gives
Since is arbitrary, we have
which assures that the sequence is bounded.
Step 2: Now we show that is a Cauchy sequence. From the construction of successive approximations, we have
on
, for
. For
, we can prove that
. Therefore, observe from (6) that
where . Using
, we obtain
To show the result for , we combine
with the inequality
to obtain
Regarding , combining
and the Holder inequality produces
Similarly to Step , from Lemma 6.2, we have
Finally, by employing and Lemma 6.3, we can proceed similarly to obtain
Substituting (16)–(20) into (15) results
where
On the other hand, note that
and recalling that , we can deduce
By similar arguments as above, we get
Therefore
where
We can also prove that
Indeed, by repeating the iteration as in, for all , we obtain
Therefore, for any , we obtain
This argument proves that is a Cauchy sequence in
Step 3: Now we prove the existence and uniqueness of the solution of Equationequation (1)(1)
(1) . One has that
as
in
. The Borel–Cantelli lemma gives us
uniformly converge to
as
, for
. Using Assumption
and
, for all
, we can prove the next inequality holds:
Note that
Similarly, we have
and
Therefore, we take the limits on both sides of (6) with respect to to obtain
We can check the uniqueness of the solution by employing the Gronwall lemma, together with a similar argument as that used in the proof of Step . This argument completes the proof.
Remark 3. We point out that the local solution exists and it is unique on for each real number
, then, existence of the solution to Equationequation (1)
(1)
(1) is global, that is,
is defined in
.
4. Exponential stability
Here, we use the Gronwall lemma and the properties of the concrete-phase space to obtain the exponential stability for the solutions of the stochastic Equationequation (1)
(1)
(1) and their maps. Other researchers have studied the stability as well (Dieye et al., Citation2017, Citation2019), but their results are based on the stochastic convolution, an approach completely detached from ours.
Definition 4.1. The mild solution of (1) is said to be -th moment exponentially stable when
if, for any initial value
,
measurable, there exist two positive real numbers
and
such that
, for all
For the sake of notational simplicity, we define the function
Now, we can introduce the main result of this paper in the following theorem.
Theorem 4.1 Suppose that and all conditions of Theorem 3.1 hold. Suppose in addition that the next two inequalities hold:
and
Then, the mild solution and the solution maps
to Equationequation (1)
(1)
(1) are
-th moments exponentially stable.
Proof. -th moment exponential stability of
: Combining (5),
,
and
, we can write
Note that
and that
By Holder inequality, it yields that
We note that ; therefore, the last inequality becomes
Recalling Lemma 6.2 and Assumption , as before, we use Holder inequality to get that
Finally, from Lemma 6.3, Assumptions and
and by Holder inequality, we can write
which implies that
If , the last two inequalities hold true with convention
. Substituting (25)–(29) into (24), we obtain
where satisfies (22). Multiplying both sides of (30) by
yields
From properties of the norm (see Appendix), we can write for any
which implies that
Recall that and
, hence
where and
. Using the Gronwall lemma, we get
which implies that
Therefore, from the condition in (23), the result of the -th moment exponential stability of solution
is satisfied.
-th moment exponential stability of
: For any
, we have
multiplying both sides of the last inequality by , we obtain
From (31), we have
Using the Gronwall lemma, we obtain
which results in
The inequality in (33) assures the exponential stability in -th moment of the solution maps
.□
4.1 Almost surely exponential stability
Definition 4.2. The mild solution of (1) is said to be almost surely exponentially stable if the following inequality is guaranteed almost surely
for any measurable initial value
.
We present the main result of this section.
Theorem 4.2. Assume that all the conditions of Theorem 4.1 hold true. Then
for any , which implies that
and
are almost surely exponentially stable.
Proof. Now we show for all
. It follows from Theorem 4.1 that
Let be the interval
, for any
. Define
Since, from assumption, and
, we can use the Markov inequality to write
The rightmost term of (34) is bounded from above by . Therefore, the Borel-Cantelli lemma assures that there exists an integer
such that, for all
,
Thus, if and
, we get
It follows that
which shows that is satisfied. The argument to prove
follows analogous reasoning. Namely, one can use the fact that the solution maps
are
-th moment exponentially stable and can conclude the result by repeating that previous reasoning. The details are omitted.□
Remark 4. The author of (Grimmer, Citation1982) presents sufficient conditions for the exponential stability of the resolvent operator . The paper (Grimmer, Citation1982) shows
and
from the contraction of the
-semigroup
and the properties of the function
by using the infinitesimal generator of the translation semigroup.
5. Example
Set and
. Consider the following neutral stochastic integro-differential equation with infinite delay and Poisson jumps of the form:
Here, denotes the standard
-valued Wiener process;
and
are continuous functions,
and
.
Let with the norm
, we define
by
with domain
. It is well known that
is the infinitesimal generator of a strongly continuous contraction semigroup
on
. Thus,
holds. Let
be the operator defined by
for
and
. The resolvent operator
decays exponentially, i.e.
, see Remark for more details. We now suppose that the next three conditions are valid.
(i) For and
.
(ii) Let , there exist real numbers
such that
for and
, where
denotes the norm of
.
(iii) Let , there exists an integrable function
such that
for
and
.
In addition, we assume that .
For and
define the operators
and
as follows:
which contain a variable delay, and
contains a distributed delay. Set for all
and all
, and
for all
and all
. Then, (36) takes the form of the system (1).
Applying together with the definition of the norm
, we have that
and by using together with the Holder inequality, we obtain
for any and
. It then follows that all the assumptions of Theorem 3.1 are satisfied with
and
. As a result, (36) has a mild solution.
For the case , by the convention
, we have the following constants
Thus, by Theorems 4.1 and 4.2, these solutions and their maps are mean square exponentially stable and almost surely exponentially stable provided that
For a numerical illustration, take and
. It follows that
The mild solutions and their maps are fourth moment exponentially stable and almost surely exponentially stable provided that
6. Conclusion
In this paper, we have studied neutral stochastic integro-differential equations with infinite delay and Poisson jumps under global Lipschitz conditions. In this study, we have used successive approximations to show the existence of mild solutions. We also prove the exponential stability of solutions and their maps. It is uncertain whether our approach copes with weaker conditions, such as local Lipschitz and non-Lipschitz conditions. The results in this paper can be seen as an extension of the ones in (Diop et al., Citation2014) because we consider the infinite-delay case here; in contrast, the authors of (Diop et al., Citation2014) have considered the finite delay case.
PUBLIC INTEREST STATEMENT
Stochastic processes are much used to represent mathematical models for phenomena and systems that vary randomly. These phenomena and systems may include the growth of bacterial population, price changes in the stock market, extinction and persistence of diseases, movement of a gas molecule and the number of phone calls. Representing such random phenomena motivates the study of stochastic differential equations. Characterizing the stability of stochastic systems has become a central topic in systems sciences. In this paper, we focus on the stability of stochastic integro-differential equations with noise.
AppendixThe next result is immediate from the definition of the phase space ![](//:0)
; see (Hino et al., Citation1991) for more details
Proposition 6.1. The phase space ![](//:0)
satisfies the next proprieties.
![](//:0)
If ![](//:0)
is such that ![](//:0)
and ![](//:0)
is continuous on ![](//:0)
, then for every ![](//:0)
, the conditions below hold:
,
.
The corresponding history
of the function
in
is a
-valued continuous function in
.
The space
is complete.
Suppose that
is a Cauchy sequence in
, and if
converges to
for
on any compact subset of the interval
, then
and
as
.
Lemma 6.2. ([20, Thm. 4.36, p. 114]). For a -valued predictable process
and for any
, we have the following inequality
Lemma 6.3. ([18]). Let and suppose that
, for any
. Then, there exists a positive real number
such that
Additional information
Funding
Notes on contributors
![](/cms/asset/2b75d0d6-6863-4a2b-9419-469d99c02849/oama_a_1979733_ilg0001.jpg)
Aziz Zouine
Aziz Zouine is a Ph.D. student at the Ibn Zohr University, Morocco. His main research area is stochastic (hybrid) nonlinear differential equations with delay. He holds a Master’s degree in Mathematics and Applications from the Cadi Ayyad University and a Bachelor's degree in Mathematics from the Ibn Zohr University, Morocco.
Hassane Bouzahir
Dr. Hassane Bouzahir is a Full Professor of Applied Mathematics/Statistics at the National School of Applied Sciences, Ibn Zohr University in Agadir, Morocco. He is the Director of a Research Laboratory on Systems Engineering and Information Technology. His research interests include linear operators theory, functional differential equations, partial functional differential equations, stochastic systems and control, modeling, and engineering simulation.
Alessandro N. Vargas
Alessandro N. Vargas received the BS degree in computer engineering from the Universidade Federal do Esprito Santo, Vitória, in 2002, and Master and Ph.D. degrees in electrical engineering from the School of Electrical and Computer Engineering, University of Campinas, Campinas, Brazil, in 2004 and 2009, respectively. Since 2007, he has been a Control Systems Professor with the Universidade Tecnológica Federal do Paraná, Paraná, Brazil. His research interests include stochastic systems and control with applications in electronics, mechatronics, and electrical engineering.
References
- Bouzahir, H., Benaid, B., & Imzegouan, C. (2017). Some stochastic functional differential equations with infinite delay: A result on existence and uniqueness of solutions in a concrete fading memory space. Chin. Journal Math. (N.Y.), 9(1). https://doi.org/https://doi.org/10.1155/2017/8219175
- Chen, H., Zhu, C., & Zhang, Y. (2014). A note on exponential stability for impulsive neutral stochastic partial functional differential equations. Applied Mathematics and Computation, 227(15), 139–13. https://doi.org/https://doi.org/10.1016/j.amc.2013.10.058
- Cui, J., & Yan, L. (2012). Successive approximation of neutral stochastic evolution equations with infinite delay and Poisson jumps. Appl. Math. Comput, 218(12), 6776–6784. https://doi.org/https://doi.org/10.1016/j.amc.2011.12.045
- Cui, J., Yan, L., & Sun, X. (2011). Exponential stability for neutral stochastic partial differential equations with delays and poisson jumps. Statistics & Probability Letters, 81(12), 1970–1977. https://doi.org/https://doi.org/10.1016/j.spl.2011.08.010
- Dieye, M., Diop, M. A., & Ezzinbi, K. (2017). On exponential stability of mild solutions for some stochastic partial integrodifferential equations. Statistics & Probability Letters, 123(123), 61–76. https://doi.org/https://doi.org/10.1016/j.spl.2016.10.031
- Dieye, M., Diop, M. A., & Ezzinbi, K. (2019). Almost sure asymptotic stability for some stochastic partial functional integrodifferential equations on Hilbert spaces. Cogent Mathematics & Statistics, 6(1), 1–15. https://doi.org/https://doi.org/10.1080/25742558.2019.1602928
- Diop, M. A., Ezzinbi, K., & Lo, M. (2014). Exponential stability for some stochastic neutral partial functional integrodifferential equations with delays and Poisson jumps. Semigroup Forum, 88(3), 595–609. https://doi.org/https://doi.org/10.1007/s00233–013–9555–y
- Govindan, T. E. (2003). Stability of mild solutions of stochastic evolution equations with variable delay. Stochastics Analysis Applications, 5(5), 1059–1077. https://doi.org/https://doi.org/10.1081/SAP–120022863
- Grimmer, R. C. (1982). Resolvent operators for integral equations in a Banach space., Trans. Amer. Math. Soc, 273(1), 333–349. https://doi.org/https://doi.org/10.1090/S0002–9947–1982–0664046–4
- Helge, H., Jan, U., Tusheng, Z., & Bernt, O. (2010). Stochastic partial differential equations: A modeling, white noise functional approach. Springer-Verlag.
- Hino, Y., Naito, T., & Murakami, S. (1991). Functional Differential Equations with Infinite Delay. Springer-Verlag. https://doi.org/https://doi.org/10.1007/BFb0084432
- Intissar., A. (2020). A mathematical study of a generalized seir model of covid-19. SciMedicine Journal, 2(38), 30–67. https://doi.org/https://doi.org/10.28991/SciMedJ –2020–02–SI–4
- Jiang, F., Yang, H., & Shen, Y. (2016). A note on exponential stability for second-order neutral stochastic partial differential equations with infinite delays in the presence of impulses. Applied Mathematics and Computation, 287-288(5), 125–133. https://doi.org/https://doi.org/10.1016/j.amc.2016.04.021
- Kai, L., & Yufeng, S. (2006). Razumikhin-type theorems of infinite dimensional stochastic functional differential equations. IFIP Int. Fed. Inf. Process. Systems, Control, Modeling and Optimization, 90(202), 237–247. https://doi.org/https://doi.org/10.1007/0–387–33882–9–22
- Kostikov, Y. A., & Romanenkov, A. M. (2020). Approximation of the multidimensional optimal control problem for the heat equation (applicable to computational fluid dynamics (cfd)). Civil Engineering Journal, 6(4), 743–768. https://doi.org/https://doi.org/10.28991/cej–2020–03091506
- Mao, W., Hu, L., & Mao, X. (2017). Neutral stochastic functional differential equations with Lévy jumps under the local Lipschitz condition. Adv. Difference Equ, 57(1), 24. https://doi.org/https://doi.org/10.1186/s13662–017–1102–9
- Mao, X. (2007). Stochastic Differential Equations and Applications. second ed. Horvood.
- Ren, Y., & Xia, N. (2009). Existence, uniqueness and stability of the solutions to neutral stochastic functional differential equations with infinite delay. Appl. Math. Comput, 1(1), 72–79. https://doi.org/https://doi.org/10.1016/j.amc.2008.11.009
- Taniguchi, T. (1998). Almost sure exponential stability for stochastic partial functional differential equations. Stochastic Analysis and Applications, 16(5), 965–975. https://doi.org/https://doi.org/10.1080/07362999808809573
- Taniguchi, T., Liu, K., & Truman, A. (2002). Existence, uniqueness, and asymptotic behavior of mild solutions to stochastic functional differential equations in Hilbert spaces. Journal of Differential Equations, 181(1), 72–91. https://doi.org/https://doi.org/10.1006/jdeq.2001.4073
- Trung, T. T. (2020). Smart city and modelling of its unorganized flows using cell machines. Civil Engineering Journal, 6(5), 954–960. https://doi.org/https://doi.org/10.28991/cej–2020–03091520
- Wu, F., Yin, G., & Mei, H. (2017). Stochastic functional differential equations with infinite delay: Existence and uniqueness of solutions, solution maps, Markov properties, and ergodicity. Journal of Differential Equations, 262(3), 1226–1252. https://doi.org/https://doi.org/10.1016/j.jde.2016.10.006
- Yue, C. (2014). Neutral stochastic functional differential equations with infinite delay and poisson jumps in the Cg space. Appl. Math. Comput, 237(15), 595–604. https://doi.org/https://doi.org/10.1016/j.amc.2014.03.079
- Zouine, A., Bouzahir, H., & Imzegouan, C. (2020). Delay-dependent stability of highly nonlinear hybrid stochastic systems with Levy noise. Nonlinear Stud, 27(4), 879–896. http://www.nonlinearstudies.com/index.php/nonlinear/article/view/2403