Skip to main content

Large time behavior of a linear delay differential equation with asymptotically small coefficient

Abstract

The linear delay differential equation x (t)=p(t)x(tr), t t 0 , is considered, where r>0 and the coefficient p:[ t 0 ,)R is continuous and small in the sense that t t + r |p(s)|ds0, t. It is shown that the large time behavior of the solutions can be described in terms of a special solution of the associated formal adjoint equation and the initial data. In the special case of the Dickman-de Bruijn equation, x (t)= x ( t 1 ) t , t1, our result yields an explicit asymptotic representation of the solutions as t.

MSC:34K06, 34K25, 11A51.

1 Introduction

The linear scalar delay differential equation

x (t)= x ( t 1 ) t ,t1,
(1.1)

plays a significant role in analytic number theory. Its special solution x:[0,)R with initial values x(t)=1 for t[0,1] is called the Dickman-de Bruijn function. The Dickman-de Bruijn function was first studied by actuary Dickman [1] and later by de Bruijn [2, 3] to estimate the proportion of smooth numbers up to a given bound. In [4] (see also [5]) van der Lune proposed some problems regarding the solutions of (1.1). The solutions to these problems, given by Bierkens, appeared in [6]. Suppose that ϕ:[0,1]R is a continuous function and let x= x ϕ denote the unique solution of (1.1) with initial values

x(t)=ϕ(t),0t1.
(1.2)

Among others, it was shown that if the limit

c(ϕ)= lim t [ t x ϕ ( t ) ]
(1.3)

exists and is finite, then its value is given by

c(ϕ)=ϕ(1) 0 1 ϕ(t)dt.
(1.4)

We emphasize that the solution presented in [6] does not imply the existence of the limit (1.3).

Our aim in this paper is twofold. First, we give an alternative proof of the limit relation (1.4) including the existence of the limit in (1.3). Second, we will show that the above result on the Dickman-de Bruijn equation (1.1) can be extended to the general linear equation

x (t)=p(t)x(tr),t t 0 ,
(1.5)

where r>0, t 0 R and p:[ t 0 ,)R is a continuous function such that

t t + r | p ( s ) | ds0as t.
(1.6)

Our main result (see Theorem 3.3 below) provides an asymptotic description of the solutions of (1.5). The asymptotic formula is given in terms of the initial data and a special solution of the associated formal adjoint equation

y (t)=p(t+r)y(t+r),t t 0 .
(1.7)

The special solution of (1.7) is eventually positive, it has bounded growth and it is unique up to a constant multiple.

The large time behavior of the solutions of (1.1) is discussed in Section 2 and our general result on the asymptotic description of the solutions (1.5) is presented in Section 3.

2 Large time behavior of the Dickman-de Bruijn equation

In this section, we prove the existence of the limit (1.3) for the solutions of (1.1).

Theorem 2.1 Suppose that ϕ:[0,1]R is a continuous function and let x ϕ be the unique solution of the initial value problem (1.1) and (1.2). Then

lim t [ t x ϕ ( t ) ] =ϕ(1) 0 1 ϕ(t)dt.
(2.1)

The proof of Theorem 2.1 will be based on the identity

t x ϕ (t) t 1 t x ϕ (s)ds=c(ϕ),t1,
(2.2)

where c(ϕ) is the constant given by (1.4). Using (1.1), it is easily shown that the derivative of the function on the left-hand side of (2.2) is equal to 0 identically on (1,). This, together with (1.4), implies (2.2). Since the proof is straightforward, we omit it.

Now we can give a simple short proof of Theorem 2.1.

Proof of Theorem 2.1 Write x ϕ =x for brevity. First we show that

M= sup t 1 [ t | x ( t ) | ] <.
(2.3)

Define

u(t)=tx(t),t0.
(2.4)

Then (2.2) can be written in the form

u(t)=c+ t 1 t u ( s ) s ds,t1,
(2.5)

where c=c(ϕ). From this, we find for t2,

| u ( t ) | |c|+ t 1 t | u ( s ) | s ds|c|+ 1 t | u ( s ) | s ds=K+ 2 t | u ( s ) | s ds,

where

K=|c|+ 1 2 | u ( s ) | s ds.

From the last inequality, by the application of Gronwall’s lemma (see, e.g., [[7], Chapter 1, Lemma 3.1]), we conclude that

| u ( t ) | Kexp ( 2 t 1 s d s ) = K t 2 ,t2.

Hence

sup t 2 | u ( t ) | t K 2 .

From this and (2.5), we find for t3,

| u ( t ) | |c|+ t 1 t | u ( s ) | s ds|c|+ K 2 .

In view of (2.4), this implies (2.3).

From (2.2) and (2.3), we obtain for t2,

| t x ( t ) c | = | t 1 t x ( s ) d s | t 1 t | x ( s ) | ds t 1 t M s ds M t 1 .

Letting t in the last inequality, we obtain

lim t [ t x ( t ) c ] =0

which is equivalent to the limit relation (2.1). □

We remark that the existence of the limit in (1.3) can also be deduced from the results by Győri and the first author (see [[8], Theorem 3.3] and its proof) and by Diblík (see [[9], Theorem 18 and Example 20]). However, the above results cannot be used to compute the value of the limit explicitly in terms of the initial data.

3 Main result

In this section, we present our main result on the large time behavior of the solutions of (1.5). First we show that under the smallness condition (1.6) the formal adjoint equation has an eventually positive solution with bounded growth.

Theorem 3.1 Suppose condition (1.6) holds. Then (1.7) has a solution y:[ t 0 ,)R which is positive for all large t and such that

lim sup t y ( t + r ) y ( t ) <.
(3.1)

We will prove Theorem 3.1 by applying a technique known from the oscillation theory of delay differential equations (see [[10], Section 2.3]).

Proof Let p + and p denote the positive part and the negative part of p, respectively, defined by

p ± (t)=max { 0 , ± p ( t ) } ,t t 0 .

Since 0 p |p|, by virtue of (1.6), there exists t 1 > t 0 such that

t t + r p (s)ds 1 e ,t t 1 .

Let C([ t 1 ,),R) be the space of continuous functions mapping [ t 1 ,) into with the topology of uniform convergence on compact subsets of [ t 1 ,). Let Ω denote the set of functions from C([ t 1 ,),R) which satisfy the system of inequalities

exp ( e t 1 t p + ( s + r ) d s ) y(t)exp ( e t 1 t p ( s + r ) d s ) ,t t 1 ,

and

y ( t + r ) y ( t ) e,t t 1 .

Clearly, Ω is a nonempty, closed and convex subset of C([ t 1 ,),R). Define the operator F:ΩC([ t 1 ,),R) by

F(y)(t)=exp ( t 1 t p ( s + r ) y ( s + r ) y ( s ) d s ) ,t t 1 ,yΩ.

It is easily verified that F is continuous and F(Ω)Ω. Furthermore, the functions from F(Ω) are uniformly bounded and equicontinuous on each compact subinterval of [ t 1 ,). Therefore, by the Arzela-Ascoli theorem, the closure of F(Ω) is compact in C([ t 1 ,),R). By the application of the Schauder-Tychonoff fixed point theorem, we conclude that there exists yΩ such that F(y)=y. It is easily seen that this fixed point y is a solution of (1.7) on [ t 1 ,) with property (3.1). Clearly, the solution y:[ t 1 ,)R can be extended backward to all t[ t 0 , t 1 ) by the method of steps. □

It should be noted that under the smallness condition (1.6), (1.7) may have a positive solution which does not satisfy condition (3.1). Indeed, the equation

y (t)=2t e 2 t 1 y(t+1),t0,

a special case of (1.7) when r=1, t 0 =0 and p(t)=2(t1) e 2 t + 1 , has the positive solution y(t)= e t 2 for which the ratio y(t+1)/y(t)= e 2 t + 1 is unbounded as t.

In the next theorem, we show that up to a constant multiple the special solution of (1.7) described in Theorem 3.1 is unique.

Theorem 3.2 Suppose condition (1.6) holds. Let y 1 and y 2 be eventually positive solutions of (1.7) satisfying condition (3.1). Then y 2 is a constant multiple of y 1 .

Proof We begin with two simple observations. First, if y is a solution of (1.7), then

y(t)=y( t 1 ) t 1 t p(s+r)y(s+r)dswhenever t t 1 t 0 .
(3.2)

Second, if y is a solution of (1.7) which is positive on [T,) for some T> t 0 and satisfies condition (3.1), then

y(t)y( t 1 )exp ( M t 1 t | p ( s + r ) | d s ) whenever t t 1 T,
(3.3)

where M is an arbitrary constant such that

M sup t T y ( t + r ) y ( t ) .
(3.4)

Indeed, from (1.7) we find for tT,

y (t)=p(t+r) y ( t + r ) y ( t ) y(t).

Hence

y(t)=y( t 1 )exp ( t 1 t p ( s + r ) y ( s + r ) y ( s ) d s ) ,t t 1 T.

This, together with (3.4), implies (3.3).

By assumptions, there exists T> t 0 such that both solutions y 1 and y 2 are positive on [T,) and satisfy condition (3.1). As noted before (see (3.3)), if M>1 is sufficiently large, then

y j (t) y j ( t 1 )exp ( M t 1 t | p ( s + r ) | d s ) whenever t t 1 T,j=1,2.
(3.5)

Since M>1, if q>0 is sufficiently small, then e q M <M. By virtue of (1.6), there exists t 1 >T such that

t t + r | p ( s ) | ds<q,t t 1 .
(3.6)

We will show that y 2 (t)=c y 1 (t) for all t t 1 , where c= y 2 ( t 1 )/ y 1 ( t 1 ). In view of the linearity of (1.7), the function y 3 =c y 1 is a solution of (1.7) and, by virtue of (3.5), the quantity

S= sup t t 1 [ | y 2 ( t ) y 3 ( t ) | exp ( M t 1 t | p ( s + r ) | d s ) ]

is finite. Applying (3.2) to both solutions y 2 and y 3 of (1.7) and taking into account that y 3 ( t 1 )= y 2 ( t 1 ), we obtain, for t t 1 ,

| y 2 ( t ) y 3 ( t ) | = | t 1 t p ( s + r ) ( y 2 ( s + r ) y 3 ( s + r ) ) d s | t 1 t | p ( s + r ) | | y 2 ( s + r ) y 3 ( s + r ) | d s S t 1 t | p ( s + r ) | exp ( M t 1 s + r | p ( u + r ) | d u ) d s = S t 1 t | p ( s + r ) | exp ( M t 1 s | p ( u + r ) | d u ) exp ( M s s + r | p ( u + r ) | d u ) d s S e q M t 1 t | p ( s + r ) | exp ( M t 1 s | p ( u + r ) | d u ) d s = S e q M [ 1 M exp ( M t 1 s | p ( u + r ) | d u ) ] t 1 t S e q M M exp ( M t 1 t | p ( u + r ) | d u ) ,

where the last but one inequality is a consequence of (3.6). From the last inequality, we obtain

| y 2 ( t ) y 3 ( t ) | exp ( M t 1 t | p ( u + r ) | d u ) S e q M M ,t t 1 .

Hence

S e q M M S.

Since e q M <M, this implies that S=0 and therefore y 2 (t)= y 3 (t)=c y 1 (t) for all t[ t 1 ,). Finally, by the uniqueness of the backward continuation of the solutions of (1.7), we conclude that y 2 (t)=c y 1 (t) for all t[ t 0 ,). □

Now we can formulate our main result about the large time behavior of the solutions of (1.5).

Theorem 3.3 Suppose condition (1.6) holds. Let x ϕ denote the solution of (1.5) with initial data

x(t)=ϕ(t), t 0 rt t 0 ,
(3.7)

where ϕ:[ t 0 r, t 0 ]R is a continuous function. Then

x ϕ (t)= 1 y ( t ) ( c ( ϕ ) + o ( 1 ) ) ,t,
(3.8)

where y is any eventually positive solution of (1.7) satisfying (3.1) and c(ϕ) is a constant given by

c(ϕ)=ϕ( t 0 )y( t 0 )+ t 0 r t 0 p(s+r)ϕ(s)y(s+r)ds.
(3.9)

As shown in Theorem 3.2, the special solution y of (1.7) in the asymptotic relation (3.8) is unique up to a constant multiple. Thus, (3.8) gives the same asymptotic representation independently of the choice of y.

Theorem 3.3 is a generalization of Theorem 2.1 to (1.5). Indeed, in the special case r=1, t 0 =1 and p(t)= t 1 ; (1.5) reduces to the Dickman-de Bruijn equation (1.1). Its formal adjoint equation

y (t)= y ( t + 1 ) t + 1 ,t1,

has the positive solution y(t)=t satisfying condition (3.1). Therefore, Theorem 3.3 applies and its conclusion reduces to the limit relation (1.3).

For qualitative results similar to Theorem 3.3, see [8, 9, 11, 12] and the references therein.

The proof of Theorem 3.3 will be based on the well-known duality between the solutions of a linear delay differential equation and its formal adjoint equation (see [[7], Section 6.3]). Namely,

x(t)y(t)+ t t + r p(s)x(sr)y(s)ds=constant
(3.10)

for t t 0 whenever x and y are solutions of (1.5) and (1.7), respectively. We will also need the following simple lemma.

Lemma 3.4 Let cR and suppose that a:[ t 0 ,)R is a continuous function such that

t t + r | a ( s ) | ds0,t.
(3.11)

Then every continuous solution of the integral equation

z(t)=c+ t t + r a(s)z(sr)ds,t t 0 ,
(3.12)

converges to c as t.

Proof Let q(0,1). By virtue of (3.11), there exists T> t 0 such that

t t + r | a ( s ) | ds<q,tT.
(3.13)

Define

M= max T r t T | z ( t ) | .

Choose a constant K such that

K>max { M , | c | ( 1 q ) 1 } .
(3.14)

Clearly, |z(t)|M<K for t[Tr,T] and we claim that

| z ( t ) | <Kfor all tTr.
(3.15)

Otherwise, there exists t 1 >T such that

| z ( t ) | <Kfor t[Tr, t 1 ) and  | z ( t 1 ) | =K.

From this and (3.12), we find that

K = | z ( t 1 ) | | c | + t 1 t 1 + r | a ( s ) | | z ( s r ) | d s | c | + K t 1 t 1 + r | a ( s ) | d s | c | + K q ,

the last inequality being a consequence of (3.13). Hence K|c| ( 1 q ) 1 , contradicting (3.14). Thus, (3.15) holds.

From (3.12) and (3.15), we find for tT,

| z ( t ) c | t t + r | a ( s ) | | z ( s r ) | dsK t t + r | a ( s ) | ds.

Letting t in the last inequality and using (3.11), we conclude that z(t)c as t. □

Now we are in a position to give a proof of Theorem 3.3.

Proof of Theorem 3.3 Write x ϕ =x for brevity and let y be a solution of (1.7) which is positive on [ t 1 ,) for some t 1 > t 0 and satisfies condition (3.1). By virtue of (3.7) and (3.10), we have

x(t)y(t)+ t t + r p(s)x(sr)y(s)ds=c(ϕ)
(3.16)

for t t 0 with c(ϕ) as in (3.9). If we let

z(t)=x(t)y(t),t t 0 ,

then (3.16) can be written in the form (3.12) with

c=c(ϕ),a(t)=p(t) y ( t ) y ( t r ) ,t t 1 +r,

and t 0 replaced with t 1 +r. Clearly, conditions (1.6) and (3.1) imply that assumption (3.11) of Lemma 3.4 is satisfied. By the application of Lemma 3.4, we conclude that

lim t z(t)= lim t [ x ( t ) y ( t ) ] =c

which is only a reformulation of the limit relation (3.8). □

Finally, we remark that applying a transformation technique described in [13] and [14], Theorem 3.3 can possibly be extended to a class of equations with time-varying delays.

References

  1. Dickman K: On the frequency of numbers containing prime factors of a certain relative magnitude. Ark. Mat. Astron. Fys. 1930, 22A(10):1-14.

    Google Scholar 

  2. de Bruijn NG: On the number of positive integers ≤ x and free of prime factors > y . I. Indag. Math. 1951, 13: 50-60.

    Article  Google Scholar 

  3. de Bruijn NG: On the number of positive integers ≤ x and free of prime factors > y . II. Indag. Math. 1966, 2: 239-247.

    Article  MathSciNet  Google Scholar 

  4. Nieuw Archief voor Wiskunde, Problem Section, 4/14 no. 3 Nov. 1996, p.429

  5. Nieuw Archief voor Wiskunde, Problem Section, 5/9 no. 2 June 2008, p.232

  6. Nieuw Archief voor Wiskunde, Problem Section, 5/11 no. 1 March 2010, p.76

  7. Hale JK: Theory of Functional Differential Equations. Springer, New York; 1977.

    Book  Google Scholar 

  8. Győri I, Pituk M: Stability criteria for linear delay differential equations. Differ. Integral Equ. 1997, 10: 841-852.

    Google Scholar 

  9. Diblík J: Behaviour of solutions of linear differential equations with delay. Arch. Math. 1998, 34: 31-47.

    Google Scholar 

  10. Győri I, Ladas G: Oscillation Theory of Delay Differential Equations with Applications. Oxford University Press, New York; 1991.

    Google Scholar 

  11. Arino O, Pituk M: More on linear differential systems with small delays. J. Differ. Equ. 2001, 170: 381-407. 10.1006/jdeq.2000.3824

    Article  MathSciNet  Google Scholar 

  12. Diblík J, Růžičková M: Asymptotic behavior of solutions and positive solutions of differential delayed equations. Funct. Differ. Equ. 2007, 14: 83-105.

    MathSciNet  Google Scholar 

  13. Čermák J: Note on canonical forms for functional differential equations. Math. Pannon. 2000, 11: 29-39.

    MathSciNet  Google Scholar 

  14. Neuman F: On transformations of differential equations and systems with deviating argument. Czechoslov. Math. J. 1981, 31: 87-90.

    MathSciNet  Google Scholar 

Download references

Acknowledgements

M Pituk was supported in part by the Hungarian National Foundation for Scientific Research (OTKA) Grant No. K101217. G Röst was supported in part by European Union and co-funded by the European Social Fund under the project ‘Telemedicine focused research activities on the field of Mathematics, Informatics and Medical sciences’ of project number TÁMOP-4.2.2.A-11/1/KONV-2012-0073, ERC Starting Grant No. 259559 and OTKA K109782.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mihály Pituk.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pituk, M., Röst, G. Large time behavior of a linear delay differential equation with asymptotically small coefficient. Bound Value Probl 2014, 114 (2014). https://doi.org/10.1186/1687-2770-2014-114

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-114

Keywords