This paper focuses on the stability and oscillations of Euler-Maclaurin method for linear differential equations with piecewise constant arguments u′(t)=au(t)+bu([t]). The necessary and sufficient conditions under which the numerical stability region contains the analytical stability region are given. Furthermore, the conditions of oscillation for the Euler-Maclaurin method are obtained. We prove that the Euler-Maclaurin method preserves the oscillations of the analytic solution. Moreover, the relationships between stability and oscillations are discussed for analytic solution and numerical solution, respectively. Finally, some numerical experiments for verifying the theoretical analysis are also provided.
1. Introduction
In the present paper we will consider the following differential equations with piecewise constant arguments (EPCA):
(1)u′(t)=au(t)+bu([t]),t≥0,u(0)=u0,
where a,b,u0∈R and [·] denotes the greatest integer function. Moreover, we assume that b<0. The general form of this type of equations is
(2)u′(t)=f(t,u(t),u(β(t))),t≥0,u(0)=u0,
where the argument β(t) has intervals of constancy.
It is well known that studies of EPCA were motivated by the fact that they represent a hybrid of continuous and discrete dynamical systems and combine the properties of both the differential and difference equations. The local asymptotic stabilization of a class of instable autonomous nonlinear difference equations was studied in [1]. Convergence and stability of adaptation algorithm for a class of nominally first-order hybrid time-invariant linear systems were considered in [2], and stability for hybrid difference-differential systems was analyzed in [3]. The first work devoted to EPCA is the paper of Shah and Wiener [4] in 1983. In 1984 Cooke and Wiener [5] studied these equations with delay. The investigation of EPCA is an important subject because they play an important role in numerous applications [6–8]. In the literature, there are many papers dealing with the properties of EPCA, including the existence and uniqueness of solutions [9], the stability [10], the oscillations [11], and the existence of integral manifolds [12]. Moreover, the concept of EPCA has been generalized in [13–15] by introducing arbitrary piecewise constant functions as arguments. Significant parts of pioneer results for EPCA can be found in [16]. For more details of EPCA, the reader can see [17–22] and the references cited therein.
In recent ten years, much work has been done in the numerical solution of EPCA. For instance, in [23, 24], the stability of numerical solution in Runge-Kutta methods for (1) and the mixed type EPCA was studied, respectively. In [25], the numerical stability of θ-methods was considered for EPCA of advanced type. For EPCA with variable coefficients, Song and Liu [26] studied the asymptotic stability of the analytic and numerical solutions. The oscillations of θ-methods and Runge-Kutta methods were investigated in [27, 28] for EPCA of retarded type, respectively. However, as far as we know, very few results were obtained on combining the stability with the oscillations of the numerical solution in the same article except for [29]. Different from [29], the novel idea of our paper is that we will study both stability and oscillations of the numerical solution in the Euler-Maclaurin method for (1), and their relationships are analyzed quantitatively.
Euler-Maclaurin method is an important tool of numerical analysis which was discovered independently and almost simultaneously by Euler and Maclaurin in the first half of 18th century. Rota [30] called Euler-Maclaurin “one of the most remarkable formulas of mathematics.” After that, it shows us how to trade a finite sum for an integral. It works much like Taylor’s formula. The equation involves an infinite series that may be truncated at any point, leaving an error term that can be bounded.
It is well known that the trapezoidal rule can be derived from the Euler-Maclaurin formula. See, for example, Munro’s paper [31]. In [32], the author indicated how the Newton-Cotes quadrature formulas and various other quadrature formulas can be developed from special cases of periodic Euler-Maclaurin formula. With respect to the numerical analysis of delay differential equations (DDEs), few results on the Euler-Maclaurin method were obtained except for [33]. In our present paper, we will apply this method to (1) and investigate numerical stability and oscillations.
The rest of this paper is organized as follows. In Section 2, some necessary notations and results are given. In Section 3, we propose Euler-Maclaurin method to solve (1), and the asymptotical stability of numerical solution is discussed. In Section 4, we investigate the oscillations of the numerical solution and prove that the oscillations of the analytic solution are inherited by the Euler-Maclaurin method. In Section 5, we establish some relationships between stability and oscillations, and several numerical experiments are shown in Section 6. Finally, concluding remarks are given in Section 7.
2. Preliminaries
In this section, we will focus our attention on some preliminary results which will be used in the further analysis.
Proposition 1 (see [16]).
A solution of (1) on [0,∞) is a function u(t) which satisfies the following conditions:
u(t)is continuous on[0,∞);
the derivativeu′(t)exists at each pointt∈[0,∞), with the possible exception of the points[t]∈[0,∞), where one-sided derivatives exist;
equation (1) is satisfied on each interval[n,n+1)⊂[0,∞)with integral end points.
Theorem 2 (see [16]).
Equation (1) has a unique solution on [0,∞):(3)u(t)=m0({t})b0[t]u0,
where {t} is the fractional part of t({t}=t-[t]) and
(4)b0=m0(1),m0(t)=eat+(eat-1)a-1b.
Theorem 3 (see [16]).
The solution of (1) is asymptotically stable (u(t)→0 as t→∞) for all u0, if and only if
(5)-a(ea+1)ea-1<b<-a,a≠0,-2<b<0,a=0.
Definition 4.
A non-trivial solution of (1) is said to be oscillatory if there exists a sequence {tk}k=1∞ such that tk→∞ as k→∞ and u(tk)u(tk-1)<0; otherwise it is called non-oscillatory. We say (1) is oscillatory if all the non-trivial solutions of (1) are oscillatory; we say (1) is non-oscillatory if all the non-trivial solutions of (1) are non-oscillatory.
Theorem 5 (see [16]).
Every solution of (1) is non-oscillatory if and only if
(6)b>aea1-ea.
Remark 6.
By Definition 4, we can see that if a solution u(t) of (1) is non-oscillatory and continuous, then it must be eventually positive or eventually negative. That is, there exists a T0∈R such that u(t)>0 or u(t)<0 for t≥T0.
3. Numerical Stability3.1. Bernoulli’s Number and Bernoulli’s Polynomial
It is well known that
(7)zez-1=∑j=0∞Bjj!zj,|z|<2π,zexzez-1=∑j=0∞Bj(x)j!zj,|z|<2π,
where Bj and Bj(x), j=0,1,…, are called Bernoulli’s number and Bernoulli’s polynomial, respectively.
Proposition 7 (see [33]).
The Bernoulli’s number Bj satisfies the following properties:
The Bernoulli’s polynomial Bj(x) satisfies the following properties:
B0(x)=1,B1(x)=x-1/2, Bk(x)=∑j=0kCkjBjxk-j,
B2k+1(1)=B2k+1(0)=B2k+1=0,
B2k(1)=B2k(0)=B2k,
Bk(x)=(1/(k+1))Bk+1′(x),k=1,2,….
3.2. The Euler-Maclaurin Method and Discretization
Let q>0 be an integer, and assume that the function f(t) is at least 2q+2 times continuously differentiable on [c,d]. Further assume that h evenly divides c and d; then Atkinson’s version of the Euler-Maclaurin formula [34] is as follows:
(8)∫cdf(t)dt=h∑n=cdf(n)-h2(f(c)+f(d))∫cdf(t)dt=-∑i=1qh2iB2i(2i)!(D2i-1f(d)-D2i-1f(c)),
where Bj denotes the jth Bernoulli number and D denotes the differentiation operator.
Let h be a given stepsize; m>1 is a given integer and satisfies h=1/m. The grid points ti are defined by ti=ih(i=0,1,2,…); let i=km+l(l=0,1,…,m-1). Applying formula (8) to (1), we have
(9)ui+1=ui+ha2(ui+1+ui)-∑j=1nB2j(ha)2j(2j)!(ui+1-ui)+hbukm.
Lemma 9 (see [35]).
Assume that f(x) has 2n+3rd continuous derivative on the interval [ti,ti+1]; then we have
(10)|∫titi+1f(t)dt-h2(f(ti+1)+f(ti))+∑j=1nB2jh2j(2j)!(f(2j-1)(ti+1)-f(2j-1)(ti))+∑j=1nB2jh2j(2j)!(f(2j-1)(ti+1)-f(2j-1)(ti))|=O(h2n+3).
According to Theorem 2.2 in [33], we obtain the following theorem for convergence.
Theorem 10.
For any given n∈N, the Euler-Maclaurin method is of order 2n+2.
Equation (9) can be written as
(11)ukm+l+1=α(x)ukm+l+ba(α(x)-1)ukm,whena≠0,ukm+l+1=ukm+l+hbukm,whena=0,
where
(12)x=ha,α(x)=1+xΦ(x),(13)Φ(x)=1-x2+∑j=1nB2jx2j(2j)!.
It is not difficult to see that (11) is equivalent to
(14)u(k+1)m=(α(x)m+ba(α(x)m-1))ukm,whena≠0,u(k+1)m=(1+b)ukm,whena=0,(15)ukm+l=(α(x)l+ba(α(x)l-1))ukm,whena≠0,ukm+l=(1+hbl)ukm,whena=0,
where 0≤l≤m-1.
3.3. Stability AnalysisDefinition 11.
Process (9) for (1) is called asymptotically stable at (a,b) for all sufficiently small h if and only if there exists a constant h0, such that for any given u0, relation (15) defines ui that satisfies ui→0 as i→∞ whenever 0<h<h0.
Definition 12.
The set, denoted by S, of all points (a,b) at which process (9) for (1) is asymptotically stable, for all sufficiently small h, is called the asymptotical stability region.
Lemma 13 (see [33]).
If |x|≤1, then φ(x)≥1/2 for x>0 and φ(x)≥1 for x≤0.
By relation (15), we can easily prove the following lemma.
Lemma 14.
For any k and 0≤l≤m-1, there exists a constant C>0 independent of k and l such that
(16)|ukm+l|≤C|ukm|.
So we have the following corollary.
Corollary 15.
ui→0 as i→∞ if and only if ukm→0 as k→∞.
According to Definition 11, Corollary 15, and relation (14), the following theorem is obtained.
Theorem 16.
The numerical solution of (1) is asymptotically stable if and only if
(17)-a(α(x)m+1)α(x)m-1<b<-a,whena≠0,-2<b<0,whena=0.
For convenience, we divide S into two parts S=S*∪S0, where
(18)S*={(a,b)∈S:a≠0},S0={(a,b)∈S∖S*:a=0},
therefore, by Definition 12 and Theorem 16, we have
(19)S*={(a,b)∈R2:-a(α(x)m+1)α(x)m-1<b<-a,a≠0},S0={(0,b)∈R2:-2<b<0}.
Let us introduce the following sets of all points (a,b)∈R2 which satisfy condition (5) by H; in a similar way, we denote H=H*∪H0, where
(20)H*={(a,b)∈R2:-a(ea+1)ea-1<b<-a,a≠0},H0={(0,b)∈R2:-2<b<0}.
We will study which conditions lead to H*⊆S* and H0⊆S0.
In the following, we assume M>|a|, which implies that |x|<1 for the stepsize h=1/m with m≥M. The following lemma will be useful to determine the stability and oscillatory conditions.
Lemma 17 (see [33]).
If |x|≤1, then
(21)Φ(x)≤xex-1
when n is even and
(22)Φ(x)≥xex-1
when n is odd.
Then the first main result for stability is as follows.
Theorem 18.
For the Euler-Maclaurin method, H*⊆S* if and only if n is odd.
Proof.
In view of (19) and (20), we know that H*⊆S* if and only if
(23)-a(α(x)m+1)α(x)m-1≤-a(ea+1)ea-1,
that is,
(24)aα(x)m-1≥aea-1,
which is equivalent to
(25)α(x)≤ex,a>0,α(x)≥ex,a<0.
In view of α(x)=1+(x/Φ(x)), thus (25) leads to
(26)Φ(x)≥xex-1.
As a consequence of Lemma 17, the proof is complete.
Theorem 19.
For all the Euler-Maclaurin method, we have H0=S0.
4. Oscillations Analysis
Similar to Definition 4, we can easily get the definitions of oscillations and nonoscillations of numerical solution; for the sake of simplicity, we omit them. Then the relationships of the nonoscillations between the integer nodes and any nodes are as follows.
Theorem 20.
{ui} and {ukm} are given by (15) and (14), respectively; then {ui} is non-oscillatory if and only if {ukm} is non-oscillatory.
Proof.
The necessity is easy to prove, so we only consider the sufficiency. If {ukm} is non-oscillatory, without loss of generality, we can assume that {ukm} is an eventually negative solution of (14); that is, there exists a K0∈R such that ukm<0 for k>K0. We will prove ukm+l<0 for all k>K0+1 and l=0,1,…,m-1.
If a<0, then 0<α(x)<1 and α(x)m≤α(x)l; hence
(27)ukm+l=(α(x)l+ba(α(x)l-1))ukmukm+l≤(α(x)m+ba(α(x)m-1))ukmukm+l=ukm+m<0.
If a>0, then α(x)>1 and α(x)-m≤α(x)-l; thus
(28)α(x)-lukm+l=(1+ba(1-α(x)-l))ukmα(x)-lukm+l≤(1+ba(1-α(x)-m))ukmα(x)-lukm+l=α(x)-mukm+m<0.
Therefore ukm+l<0. This completes the proof.
From Theorem 20, we get the following corollary.
Corollary 21.
{ui} and {ukm} are given by (15) and (14), respectively; then {ui} is oscillatory if and only if {ukm} is oscillatory.
Theorem 22.
Equation (14) is oscillatory if and only if
(29)b<aα(x)m1-α(x)m,whena≠0,b<-1,whena=0.
Proof.
Equation (14) is oscillatory if and only if the corresponding characteristic equation has no positive roots, that is,
(30)α(x)m+ba(α(x)m-1)<0,whena≠0,1+b<0,whena=0,
which is equivalent to
(31)b<aα(x)m1-α(x)m,whena≠0,(32)b<-1,whena=0,
thus (29) is obtained. The proof is finished.
Let
(33)Π=aea1-ea,Π(m)=aα(x)m1-α(x)m,
we have the following lemma.
Lemma 23.
(I)Π(m)→Π as h→0;
(II)Π≥Π(m) if either of the following conditions is satisfied:
(i)a>0 and ex≥α(x),
(ii)a<0 and ex≤α(x);
(III)Π<Π(m) if either of the following conditions is satisfied:
(i)a>0 and ex<α(x),
(ii)a<0 and ex>α(x).
Proof.
(II) If a>0 and ex≥α(x), then ea≥α(x)m, which is equivalent to
(34)11-ea≥11-α(x)m,that is,
(35)ea1-ea≥α(x)m1-α(x)m,
so we have Π≥Π(m). Similarly, we can prove the other cases.
Definition 24.
We say the Euler-Maclaurin method preserves the oscillations of (1) if (1) oscillates, then there is an h0>0, such that (15) oscillates for h<h0. Similarly, we say the Euler-Maclaurin method preserves the nonoscillations of (1) if (1) nonoscillates then there is an h0>0, such that (15) non-oscillates for h<h0.
By Theorems 5, 20, and 22 and Corollary 21, we have the following theorem.
Theorem 25.
(i) The Euler-Maclaurin method preserves the oscillations of (1) if and only if Π<Π(m).
(ii) The Euler-Maclaurin method preserves the non-oscillations of (1) if and only if Π≥Π(m).
Applying Theorem 25 and Lemmas 17 and 23, we can obtain the following main theorems for oscillations.
Theorem 26.
The Euler-Maclaurin method preserves the oscillations of (1) if n is even.
Theorem 27.
The Euler-Maclaurin method preserves the non-oscillations of (1) if n is odd.
5. Relationships between Stability and Oscillations
Let
(36)Π~=-a(ea+1)ea-1,Π~(m)=-a(α(x)m+1)α(x)m-1.
By Theorems 3, 5, 16, and 22, we can get the following other main theorems.
Theorem 28.
When a≠0, the analytic solution of (1) is
oscillatory and unstable if b∈(-∞,Π~),
oscillatory and asymptotically stable if b∈(Π~,Π),
non-oscillatory and asymptotically stable if b∈(Π,-a),
non-oscillatory and unstable if b∈(-a,0).
Theorem 29.
When a≠0, the numerical solution of (1) is
oscillatory and unstable if b∈(-∞,Π~(m)),
oscillatory and asymptotically stable if b∈(Π~(m),Π(m)),
non-oscillatory and asymptotically stable if b∈(Π(m),-a),
non-oscillatory and unstable if b∈(-a,0).
Theorem 30.
When a=0, the analytic solution and the numerical solution of (1) both are
oscillatory and unstable if b∈(-∞,-2),
oscillatory and asymptotically stable if b∈(-2,-1),
non-oscillatory and asymptotically stable if b∈(-1,0).
Remark 31.
In Theorem 30, we can see that the solution of (1) is not non-oscillatory and unstable at the same time. Otherwise, if the solution of (1) is both non-oscillatory and unstable, we can get b∈(0,+∞); it is inconsistent with b<0.
6. Numerical Experiments
In this section, we will use the following four equations to demonstrate the main theorems
(37)u′(t)=u(t)-1.5u([t]),u(0)=1,(38)u′(t)=-2u(t)-u([t]),u(0)=1,(39)u′(t)=3u(t)-3.2u([t]),u(0)=1,(40)u′(t)=0.8u(t)-0.1u([t]),u(0)=1.
For (37) and (38), it is not difficult to see that (1,-1.5),(-2,-1)∈H by straightforward computing. Let the stepsize h=1/m; we shall use the Euler-Maclaurin method with n=2, the θ-methods with θ=0.5, and the 2-Radau IA method to get the numerical solution at t=10. On the other hand, the theoretical solutions are u(10)≈3.0750×10-9 and u(10)≈5.3398×10-6 for (37) and (38), respectively. In Tables 1 and 2, we have listed the absolute errors (AE), the relative errors (RE) at t=10, and the ratio of the errors of the case m=20 over that of m=40. We can see from these tables that the Euler-Maclaurin method with n=2 is of order 6; that is, the method preserves its order of convergence. Furthermore, the errors of the θ-methods and the 2-Radau IA method are both larger than that of the Euler-Maclaurin method. Therefore, compared with the two mentioned numerical methods, the Euler-Maclaurin method has higher accuracy.
Errors and ratio for (37).
Euler-Maclaurin method
θ-methods
2-Radau IA method
AE
RE
AE
RE
AE
RE
m=2
1.5555E-13
5.0586E-5
2.7882E-9
9.0670E-1
6.5320E-10
2.1240E-1
m=3
1.3547E-14
4.4054E-6
1.8944E-9
6.1610E-1
1.7228E-10
5.6000E-2
m=5
6.2942E-16
2.0469E-7
8.6318E-10
2.8070E-1
3.5051E-11
1.1400E-2
m=10
9.8172E-18
3.1926E-9
2.3894E-10
7.7700E-2
4.2380E-12
1.4000E-3
m=20
1.5337E-19
4.9876E-11
6.1287E-11
1.9900E-2
5.2217E-13
1.6981E-4
m=40
2.7140E-21
8.8259E-13
1.5420E-11
5.0000E-3
6.4824E-14
2.1081E-5
Ratio
56.5107
56.5110
3.9745
3.9800
8.0552
8.0551
Errors and ratio for (38).
Euler-Maclaurin method
θ-methods
2-Radau IA method
AE
RE
AE
RE
AE
RE
m=2
2.5585E-9
4.7914E-4
1.1595E-5
2.1715E+0
8.9866E-7
1.6830E-1
m=3
2.1747E-10
4.0727E-5
3.5420E-6
6.6330E-1
2.6576E-7
4.9800E-2
m=5
9.9800E-12
1.8690E-6
1.0695E-6
2.0030E-1
5.9321E-8
1.1100E-2
m=10
1.5485E-13
2.9000E-8
2.4905E-7
4.6600E-2
7.7177E-9
1.4000E-3
m=20
2.4153E-15
4.5232E-10
6.1184E-8
1.1500E-2
9.8799E-10
1.8502E-4
m=40
3.7816E-17
7.0819E-12
1.5230E-8
2.9000E-3
1.2508E-10
2.3424E-5
Ratio
63.8698
63.8699
4.0173
3.9655
7.8989
7.8987
In Figures 1 and 3, we draw the numerical solution with m=20 for (37) and (38), respectively. It is easy to see that the numerical solution is asymptotically stable. In Figure 2, we draw the stability region for (37); it can be seen that the analytical stability region is contained in the numerical stability region, which is in accordance with Theorem 18.
The numerical solution of (37) with n=2 and m=20.
The stability region of (37): (a) the analytical stability region, (b) the numerical stability region with n=3 and m=20.
The numerical solution of (38) with n=3 and m=20.
For (39) and (40), the analytic solution of (39) is oscillatory, and the analytic solution of (40) is non-oscillatory according to Theorem 5. From Figures 4 and 5, we can see that the numerical solution of (39) is oscillatory and the numerical solution of (40) is non-oscillatory, which are in accordance with Theorems 26 and 27.
The analytic solution and the numerical solution of (39) with n=2 and m=40.
The analytic solution and the numerical solution of (40) with n=3 and m=40.
In Figure 4, let m=40; it is easy to calculate that Π≈-3.1572, Π~≈-3.3144, Π(m)≈-3.1572, and Π~(m)≈-3.3144. Obviously, b=-3.2∈(Π~,Π) and b=-3.2∈(Π~(m),Π(m)). Therefore, the analytic solution and the numerical solution of (39) are both oscillatory and asymptotically stable according to Theorems 28 and 29, which are in agreement with Figure 4. For (37), (38), and (40), we can verify them in the same way (see Figures 1, 3, and 5).
From Figures 1–5, we can see that the numerical results coincide with our theoretical analysis.
7. Conclusions
In this paper, the Euler-Maclaurin method is presented for solving an important EPCA. Asymptotical stability and oscillations of numerical solution are considered, respectively. Furthermore, the relationships between stability and oscillations are also revealed. Results of numerical examples have shown that the Euler-Maclaurin method has good convergence and stability, so it is applicable for the study of EPCA. That is, the Euler-Maclaurin method is an alternative numerical method for DDEs except for the classical θ-methods and Runge-Kutta methods. We will consider the nonlinear and multidimensional problems in the future work.
Acknowledgments
The first author’s work was supported by the National Natural Science Foundation of China (no. 11201084). The authors would like to thank the editor Chengming Huang and the reviewers for their useful comments which helped to improve the original version of the paper. The first author is also grateful to Professor Mingzhu Liu, Minghui Song, and Doctor Zhanwen Yang for their helpful comments and constructive suggestions.
de la SenM.About the stabilization of a nonlinear perturbed difference equation201220121932030210.1155/2012/320302MR2904563ZBL1244.39015de La SenM.On the robust adaptive stabilization of a class of nominally first-order hybrid systems199944359760210.1109/9.751359MR1680171ZBL1056.93616MarchenkoV. M.LuazoJ. J.On the stability of hybrid differential-difference systems200945572874010.1134/S0012266109050139MR2596705ShahS. M.WienerJ.Advanced differential equations with piecewise constant argument deviations19836467170310.1155/S0161171283000599MR729389ZBL0534.34067CookeK. L.WienerJ.Retarded differential equations with piecewise constant delays198499126529710.1016/0022-247X(84)90248-8MR732717ZBL0557.34059DaiL.SinghM. C.On oscillatory motion of spring-mass systems subjected to piecewise constant forces1994173221723210.1006/jsvi.1994.1227MuroyaY.Persistence, contractivity and global stability in logistic equations with piecewise constant delays2002270260263510.1016/S0022-247X(02)00095-1MR1916600ZBL1012.34076LiangH. H.WangG. Q.Oscillation criteria of certain third-order differential equation with piecewise constant argument201220121849807310.1155/2012/498073AkhmetM. U.BüyükadalıC.Differential equations with state-dependent piecewise constant argument201072114200421010.1016/j.na.2010.01.050MR2606777ZBL1195.34103AkhmetM. U.AruğaslanD.YılmazE.Stability in cellular neural networks with a piecewise constant argument201023392365237310.1016/j.cam.2009.10.021MR2577773ZBL1191.68484FuX.LiX.Oscillation of higher order impulsive differential equations of mixed type with constant argument at fixed time2008485-677678610.1016/j.mcm.2007.11.006MR2451113ZBL1156.34318AkhmetM. U.On the reduction principle for differential equations with piecewise constant argument of generalized type2007336164666310.1016/j.jmaa.2007.03.010MR2348532ZBL1134.34048AkhmetM. U.Integral manifolds of differential equations with piecewise constant argument of generalized type200766236738310.1016/j.na.2005.11.032MR2279531ZBL1122.34054AkhmetM. U.Stability of differential equations with piecewise constant arguments of generalized type200868479480310.1016/j.na.2006.11.037MR2382298ZBL1173.34042AkhmetM. U.ArugaslancD.YılmazE.Stability analysis of recurrent neural networks with piecewise constant argument of generalized type201023780581110.1016/j.neunet.2010.05.006WienerJ.1993River Edge, NJ, USAWorld Scientific Publishingxiv+410MR1233560Ait DadsE.LhachimiL.Pseudo almost periodic solutions for equation with piecewise constant argument2010371284285410.1016/j.jmaa.2010.06.032MR2670161ZBL1206.34094MuroyaY.New contractivity condition in a population model with piecewise constant arguments20083461658110.1016/j.jmaa.2008.05.025MR2428273ZBL1161.34048AkhmetM. U.YılmazE.Impulsive Hopfield-type neural network system with piecewise constant argument20101142584259310.1016/j.nonrwa.2009.09.003MR2661925ZBL1202.92001DimbourW.Almost automorphic solutions for differential equations with piecewise constant argument in a Banach space20117462351235710.1016/j.na.2010.11.038MR2781763ZBL1216.34075ZhangC.JiangL.Remotely almost periodic solutions to systems of differential equations with piecewise constant argument200821876176810.1016/j.aml.2007.08.007MR2436161ZBL1188.34095AkhmetM. U.Asymptotic behavior of solutions of differential equations with piecewise constant arguments200821995195610.1016/j.aml.2007.10.008MR2436530ZBL1152.34348LiuM. Z.SongM. H.YangZ. W.Stability of Runge-Kutta methods in the numerical solution of equation u'(t)=au(t)+a0u([t])2004166236137010.1016/j.cam.2003.04.002MR2041187ZBL1052.65070WangQ.ZhuQ.Stability analysis of Runge-Kutta methods for differential equations with piecewise continuous arguments of mixed type20118851052106610.1080/00207160.2010.492425MR2775098ZBL1220.65105SongM. H.YangZ. W.LiuM. Z.Stability of θ-methods for advanced differential equations with piecewise continuous arguments2005499-101295130110.1016/j.camwa.2005.02.002MR2149481SongM.LiuM. Z.Stability of analytic and numerical solutions for differential equations with piecewise continuous arguments201220121425832910.1155/2012/258329MR2914886ZBL1244.34097LiuM. Z.GaoJ.YangZ. W.Oscillation analysis of numerical solution in the θ-methods for equation x'(t)+ax(t)+a1x([t-1])=02007186156657810.1016/j.amc.2006.07.119MR2314517LiuM. Z.GaoJ. F.YangZ. W.Preservation of oscillations of the Runge-Kutta method for equation x'(t)+ax(t)+a1x([t-1])=020095861113112510.1016/j.camwa.2009.07.030MR2554343ZBL1189.65143WangQ.ZhuQ.LiuM.Stability and oscillations of numerical solutions for differential equations with piecewise continuous arguments of alternately advanced and retarded type201123551542155210.1016/j.cam.2010.08.041MR2728110ZBL1207.65103RotaG.-C.Combinatorial snapshots199921281410.1007/BF03024840MR1684850ZBL1052.05500MunroW. D.Classroom notes: note on the Euler-MacLaurin formula195865320120310.2307/2310068MR1529879BerndtB. C.SchoenfeldL.Periodic analogues of the Euler-Maclaurin and Poisson summation formulas with applications to number theory19752812368MR0384713LvW. J.YangZ. W.LiuM. Z.Stability of the Euler-Maclaurin methods for neutral differential equations with piecewise continuous arguments200718621480148710.1016/j.amc.2006.07.158MR2316941ZBL1175.65088AtkinsonK. E.19892ndNew York, NY, USAJohn Wiley & Sonsxvi+693MR1007135StoerJ.BulirschR.1993122ndNew York, NY, USASpringerxiv+660Texts in Applied MathematicsMR1295246