SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series A Tribute to Professor Ivan Kiguradze.

Open Access Research

Oscillatory behavior of second-order nonlinear neutral differential equations with distributed deviating arguments

Tongxing Li1*, Blanka Baculíková2 and Jozef Džurina2

Author Affiliations

1 Qingdao Technological University, Feixian, Shandong, 273400, P.R. China

2 Department of Mathematics, Faculty of Electrical Engineering and Informatics, Technical University of Košice, Letná 9, Košice, 042 00, Slovakia

For all author emails, please log on.

Boundary Value Problems 2014, 2014:68  doi:10.1186/1687-2770-2014-68


Dedicated to Professor Ivan Kiguradze

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/68


Received:16 January 2014
Accepted:6 March 2014
Published:24 March 2014

© 2014 Li et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

We study oscillatory properties of a class of second-order nonlinear neutral functional differential equations with distributed deviating arguments. On the basis of less restrictive assumptions imposed on the neutral coefficient, some new criteria are presented. Three examples are provided to illustrate these results.

MSC: 34C10, 34K11.

Keywords:
oscillation; neutral differential equation; second-order equation; distributed deviating argument

1 Introduction

This paper is concerned with oscillation of the second-order nonlinear functional differential equation

( r ( t ) | z ( t ) | α 1 z ( t ) ) + a b q ( t , ξ ) | x [ g ( t , ξ ) ] | α 1 x [ g ( t , ξ ) ] d σ ( ξ ) = 0 , (1.1)

where t t 0 > 0 , α 1 is a constant, and z : = x + p x τ . Throughout, we assume that the following hypotheses hold:

(H1) I : = [ t 0 , ) , r , p C 1 ( I , R ) , r ( t ) > 0 , and p ( t ) 0 ;

(H2) q C ( I × [ a , b ] , [ 0 , ) ) and q ( t , ξ ) is not eventually zero on any [ t μ , ) × [ a , b ] , t μ I ;

(H3) g C ( I × [ a , b ] , [ 0 , ) ) , lim inf t g ( t , ξ ) = , and g ( t , a ) g ( t , ξ ) for ξ [ a , b ] ;

(H4) τ C 2 ( I , R ) , τ ( t ) > 0 , lim t τ ( t ) = , and g ( τ ( t ) , ξ ) = τ [ g ( t , ξ ) ] ;

(H5) σ C ( [ a , b ] , R ) is nondecreasing and the integral of (1.1) is taken in the sense of Riemann-Stieltijes.

By a solution of (1.1), we mean a function x C ( [ t x , ) , R ) for some t x t 0 , which has the properties that z C 1 ( [ t x , ) , R ) , r | z | α 1 z C 1 ( [ t x , ) , R ) , and satisfies (1.1) on [ t x , ) . We restrict our attention to those solutions x of (1.1) which exist on [ t x , ) and satisfy sup { | x ( t ) | : t T } > 0 for any T t x . A solution x of (1.1) is termed oscillatory if it is neither eventually positive nor eventually negative; otherwise, it is called nonoscillatory. Equation (1.1) is said to be oscillatory if all its solutions oscillate.

As is well known, neutral differential equations have a great number of applications in electric networks. For instance, they are frequently used in the study of distributed networks containing lossless transmission lines, which rise in high speed computers, where the lossless transmission lines are used to interconnect switching circuits; see [1]. Hence, there has been much research activity concerning oscillatory and nonoscillatory behavior of solutions to different classes of neutral differential equations, we refer the reader to [2-30] and the references cited therein.

In the following, we present some background details that motivate our research. Recently, Baculíková and Lacková [6], Džurina and Hudáková [12], Li et al.[15,18], and Sun et al.[22] established some oscillation criteria for the second-order half-linear neutral differential equation

( r ( t ) | z ( t ) | α 1 z ( t ) ) + q ( t ) | x ( δ ( t ) ) | α 1 x ( δ ( t ) ) = 0 ,

where z : = x + p x τ ,

0 p ( t ) < 1 or p ( t ) > 1 .

Baculíková and Džurina [4,5] and Li et al.[17] investigated oscillatory behavior of a second-order neutral differential equation

( r ( t ) ( x ( t ) + p ( t ) x [ τ ( t ) ] ) ) + q ( t ) x [ σ ( t ) ] = 0 ,

where

0 p ( t ) p 0 < and τ ( t ) τ 0 > 0 . (1.2)

Ye and Xu [26] and Yu and Fu [27] considered oscillation of the second-order differential equation

( x ( t ) + p ( t ) x ( t τ ) ) + a b q ( t , ξ ) x ( g ( t , ξ ) ) d σ ( ξ ) = 0 .

Assuming 0 p ( t ) < 1 , Thandapani and Piramanantham [23], Wang [24], Xu and Weng [25], and Zhao and Meng [30] studied oscillation of an equation

( r ( t ) ( x ( t ) + p ( t ) x ( t τ ) ) ) + a b q ( t , ξ ) f ( x ( g ( t , ξ ) ) ) d σ ( ξ ) = 0 .

As yet, there are few results regarding the study of oscillatory properties of (1.1) under the conditions p ( t ) 1 or lim t p ( t ) = . Thereinto, Li and Thandapani [19] obtained several oscillation results for (1.1) in the case where (1.2) holds, σ ( ξ ) = ξ , and

t 0 d t r 1 / α ( t ) = . (1.3)

In the subsequent sections, we shall utilize the Riccati substitution technique and some inequalities to establish several new oscillation criteria for (1.1) assuming that (1.3) holds or

t 0 d t r 1 / α ( t ) < . (1.4)

All functional inequalities are assumed to hold eventually, that is, they are satisfied for all t large enough.

2 Main results

In what follows, we use the following notation for the convenience of the reader:

Q ( t , ξ ) : = min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } , d + ( t ) : = max { 0 , d ( t ) } , ϕ ( t ) : = α p [ h ( t ) ] h ( t ) p [ h ( t ) ] τ ( t ) τ ( t ) , ζ ( t ) : = ρ + ( t ) ρ ( t ) + ϕ ( t ) , φ ( t ) : = ( ρ + ( t ) ρ ( t ) ) α + 1 + p α [ h ( t ) ] ( ζ + ( t ) ) α + 1 τ ( t ) , and δ ( t ) : = η ( t ) d s r 1 / α ( s ) ,

where h, ρ, and η will be specified later.

Theorem 2.1Assume (H1)-(H5), (1.3), and let g ( t , a ) C 1 ( I , R ) , g ( t , a ) > 0 , g ( t , a ) t , and g ( t , a ) τ ( t ) for t I . Suppose further that there exists a real-valued function h C 1 ( I , R ) such that p [ g ( t , ξ ) ] p [ h ( t ) ] for t I and ξ [ a , b ] . If there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that

lim sup t t 0 t ρ ( s ) [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 r [ g ( s , a ) ] φ ( s ) ( α + 1 ) α + 1 ( g ( s , a ) ) α ] d s = , (2.1)

then (1.1) is oscillatory.

Proof Let x be a nonoscillatory solution of (1.1). Without loss of generality, we assume that there exists a t 1 I such that x ( t ) > 0 , x [ τ ( t ) ] > 0 , and x [ g ( t , ξ ) ] > 0 for all t t 1 and ξ [ a , b ] . Then z ( t ) > 0 . Applying (1.1), one has, for all sufficiently large t,

( r ( t ) | z ( t ) | α 1 z ( t ) ) + a b q ( t , ξ ) x α [ g ( t , ξ ) ] d σ ( ξ ) + a b q ( τ ( t ) , ξ ) p α [ h ( t ) ] x α [ g ( τ ( t ) , ξ ) ] d σ ( ξ ) + p α [ h ( t ) ] τ ( t ) ( r [ τ ( t ) ] | z [ τ ( t ) ] | α 1 z [ τ ( t ) ] ) = 0 .

Using the inequality (see [[5], Lemma 1])

( A + B ) α 2 α 1 ( A α + B α ) , for  A 0 , B 0 ,  and  α 1 ,

the definition of z, g ( τ ( t ) , ξ ) = τ [ g ( t , ξ ) ] , and p [ g ( t , ξ ) ] p [ h ( t ) ] , we conclude that

( r ( t ) | z ( t ) | α 1 z ( t ) ) + 1 2 α 1 a b Q ( t , ξ ) z α [ g ( t , ξ ) ] d σ ( ξ ) + p α [ h ( t ) ] τ ( t ) ( r [ τ ( t ) ] | z [ τ ( t ) ] | α 1 z [ τ ( t ) ] ) 0 . (2.2)

By virtue of (1.1), we get

( r ( t ) | z ( t ) | α 1 z ( t ) ) 0 , t t 1 . (2.3)

Thus, r | z | α 1 z is nonincreasing. Now we have two possible cases for the sign of z : (i) z < 0 eventually, or (ii) z > 0 eventually.

(i) Assume that z ( t ) < 0 for t t 2 t 1 . Then we have by (2.3)

r ( t ) | z ( t ) | α 1 z ( t ) r ( t 2 ) | z ( t 2 ) | α 1 z ( t 2 ) < 0 , t t 2 ,

which yields

z ( t ) z ( t 2 ) r 1 / α ( t 2 ) | z ( t 2 ) | t 2 t r 1 / α ( s ) d s .

Then we obtain lim t z ( t ) = due to (1.3), which is a contradiction.

(ii) Assume that z ( t ) > 0 for t t 2 t 1 . It follows from (2.2) and g ( t , ξ ) g ( t , a ) that

( r ( t ) ( z ( t ) ) α ) + p α [ h ( t ) ] τ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) + 1 2 α 1 z α [ g ( t , a ) ] a b Q ( t , ξ ) d σ ( ξ ) 0 . (2.4)

We define a Riccati substitution

ω ( t ) : = ρ ( t ) r ( t ) ( z ( t ) ) α ( z [ g ( t , a ) ] ) α , t t 2 . (2.5)

Then ω ( t ) > 0 . From (2.3) and g ( t , a ) t , we have

z [ g ( t , a ) ] ( r ( t ) / r [ g ( t , a ) ] ) 1 / α z ( t ) . (2.6)

Differentiating (2.5), we get

ω ( t ) = ρ ( t ) r ( t ) ( z ( t ) ) α ( z [ g ( t , a ) ] ) α + ρ ( t ) ( r ( t ) ( z ( t ) ) α ) ( z [ g ( t , a ) ] ) α α ρ ( t ) r ( t ) ( z ( t ) ) α z α 1 [ g ( t , a ) ] z [ g ( t , a ) ] g ( t , a ) ( z [ g ( t , a ) ] ) 2 α . (2.7)

Therefore, by (2.5), (2.6), and (2.7), we see that

ω ( t ) ρ ( t ) ρ ( t ) ω ( t ) + ρ ( t ) ( r ( t ) ( z ( t ) ) α ) ( z [ g ( t , a ) ] ) α α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ω ( α + 1 ) / α ( t ) . (2.8)

Similarly, we introduce another Riccati transformation:

υ ( t ) : = ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ( z [ g ( t , a ) ] ) α , t t 2 . (2.9)

Then υ ( t ) > 0 . From (2.3) and g ( t , a ) τ ( t ) , we obtain

z [ g ( t , a ) ] ( r [ τ ( t ) ] / r [ g ( t , a ) ] ) 1 / α z [ τ ( t ) ] . (2.10)

Differentiating (2.9), we have

υ ( t ) = ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ( z [ g ( t , a ) ] ) α + ρ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ g ( t , a ) ] ) α α ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α z α 1 [ g ( t , a ) ] z [ g ( t , a ) ] g ( t , a ) ( z [ g ( t , a ) ] ) 2 α . (2.11)

Therefore, by (2.9), (2.10), and (2.11), we find

υ ( t ) ρ ( t ) ρ ( t ) υ ( t ) + ρ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ g ( t , a ) ] ) α α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] υ ( α + 1 ) / α ( t ) . (2.12)

Combining (2.8) and (2.12), we get

ω ( t ) + p α [ h ( t ) ] τ ( t ) υ ( t ) ρ ( t ) ( r ( t ) ( z ( t ) ) α ) + p α [ h ( t ) ] τ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ g ( t , a ) ] ) α + ρ ( t ) ρ ( t ) ω ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ω ( α + 1 ) / α ( t ) + p α [ h ( t ) ] τ ( t ) ρ ( t ) ρ ( t ) υ ( t ) p α [ h ( t ) ] τ ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] υ ( α + 1 ) / α ( t ) .

It follows from (2.4) that

ω ( t ) + p α [ h ( t ) ] τ ( t ) υ ( t ) ρ ( t ) 2 α 1 a b Q ( t , ξ ) d σ ( ξ ) + ρ + ( t ) ρ ( t ) ω ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ω ( α + 1 ) / α ( t ) + p α [ h ( t ) ] τ ( t ) ρ + ( t ) ρ ( t ) υ ( t ) p α [ h ( t ) ] τ ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] υ ( α + 1 ) / α ( t ) .

Integrating the latter inequality from t 2 to t, we obtain

ω ( t ) ω ( t 2 ) + p α [ h ( t ) ] τ ( t ) υ ( t ) p α [ h ( t 2 ) ] τ ( t 2 ) υ ( t 2 ) t 2 t ρ ( s ) 2 α 1 a b Q ( s , ξ ) d σ ( ξ ) d s + t 2 t [ ρ + ( s ) ρ ( s ) ω ( s ) α g ( s , a ) ρ 1 / α ( s ) r 1 / α [ g ( s , a ) ] ω ( α + 1 ) / α ( s ) ] d s + t 2 t p α [ h ( s ) ] τ ( s ) { [ ρ + ( s ) ρ ( s ) + ϕ ( s ) ] + υ ( s ) α g ( s , a ) ρ 1 / α ( s ) r 1 / α [ g ( s , a ) ] υ ( α + 1 ) / α ( s ) } d s . (2.13)

Define

A : = [ α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ] α / ( α + 1 ) ω ( t ) and B : = [ α α + 1 ρ + ( t ) ρ ( t ) [ α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ] α / ( α + 1 ) ] α .

Using the inequality

α + 1 α A B 1 / α A ( α + 1 ) / α 1 α B ( α + 1 ) / α , for  A 0  and  B 0 , (2.14)

we get

ρ + ( t ) ρ ( t ) ω ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ω ( α + 1 ) / α ( t ) 1 ( α + 1 ) α + 1 r [ g ( t , a ) ] ( ρ + ( t ) ) α + 1 ( ρ ( t ) g ( t , a ) ) α .

On the other hand, define

A : = [ α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ] α / ( α + 1 ) υ ( t ) and B : = [ α α + 1 ζ + ( t ) [ α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] ] α / ( α + 1 ) ] α .

Then we have by (2.14)

ζ + ( t ) υ ( t ) α g ( t , a ) ρ 1 / α ( t ) r 1 / α [ g ( t , a ) ] υ ( α + 1 ) / α ( t ) 1 ( α + 1 ) α + 1 r [ g ( t , a ) ] ( ζ + ( t ) ) α + 1 ρ ( t ) ( g ( t , a ) ) α .

Thus, from (2.13), we get

ω ( t ) ω ( t 2 ) + p α [ h ( t ) ] τ ( t ) υ ( t ) p α [ h ( t 2 ) ] τ ( t 2 ) υ ( t 2 ) t 2 t ρ ( s ) { a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 r [ g ( s , a ) ] ( α + 1 ) α + 1 ( g ( s , a ) ) α × [ ( ρ + ( s ) ρ ( s ) ) α + 1 + p α [ h ( s ) ] ( ζ + ( s ) ) α + 1 τ ( s ) ] } d s ,

which contradicts (2.1). This completes the proof. □

Assuming (1.2), where p 0 and τ 0 are constants, we obtain the following result.

Theorem 2.2Suppose (H1)-(H5), (1.2), (1.3), and let g ( t , a ) C 1 ( I , R ) , g ( t , a ) > 0 , g ( t , a ) t , and g ( t , a ) τ ( t ) for t I . If there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that

lim sup t t 0 t [ ρ ( s ) a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 1 + p 0 α τ 0 ( α + 1 ) α + 1 r [ g ( s , a ) ] ( ρ + ( s ) ) α + 1 ( ρ ( s ) g ( s , a ) ) α ] d s = , (2.15)

then (1.1) is oscillatory.

Proof As above, let x be an eventually positive solution of (1.1). Proceeding as in the proof of Theorem 2.1, we have z ( t ) > 0 , (2.3), and (2.4) for all sufficiently large t. Using (1.2), (2.3), and (2.4), we obtain

( r ( t ) ( z ( t ) ) α ) + p 0 α τ 0 ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) + 1 2 α 1 z α [ g ( t , a ) ] a b Q ( s , ξ ) d σ ( ξ ) 0 . (2.16)

The remainder of the proof is similar to that of Theorem 2.1, and hence it is omitted. □

Theorem 2.3Suppose we have (H1)-(H5), (1.3), and let τ ( t ) t and g ( t , a ) τ ( t ) for t I . Assume also that there exists a real-valued function h C 1 ( I , R ) such that p [ g ( t , ξ ) ] p [ h ( t ) ] for t I and ξ [ a , b ] . If there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that

lim sup t t 0 t ρ ( s ) [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 r [ τ ( s ) ] φ ( s ) ( α + 1 ) α + 1 ( τ ( s ) ) α ] d s = , (2.17)

then (1.1) is oscillatory.

Proof Let x be a nonoscillatory solution of (1.1). Without loss of generality, we assume that there exists a t 1 I such that x ( t ) > 0 , x [ τ ( t ) ] > 0 , and x [ g ( t , ξ ) ] > 0 for all t t 1 and ξ [ a , b ] . As in the proof of Theorem 2.1, we obtain (2.3) and (2.4). In view of (2.3), r | z | α 1 z is nonincreasing. Now we have two possible cases for the sign of z : (i) z < 0 eventually, or (ii) z > 0 eventually.

(i) Suppose that z ( t ) < 0 for t t 2 t 1 . Then, with a proof similar to the proof of case (i) in Theorem 2.1, we obtain a contradiction.

(ii) Suppose that z ( t ) > 0 for t t 2 t 1 . We define a Riccati substitution

ω ( t ) : = ρ ( t ) r ( t ) ( z ( t ) ) α ( z [ τ ( t ) ] ) α , t t 2 . (2.18)

Then ω ( t ) > 0 . From (2.3) and τ ( t ) t , we have

z [ τ ( t ) ] ( r ( t ) / r [ τ ( t ) ] ) 1 / α z ( t ) . (2.19)

Differentiating (2.18), we obtain

ω ( t ) = ρ ( t ) r ( t ) ( z ( t ) ) α ( z [ τ ( t ) ] ) α + ρ ( t ) ( r ( t ) ( z ( t ) ) α ) ( z [ τ ( t ) ] ) α α ρ ( t ) r ( t ) ( z ( t ) ) α z α 1 [ τ ( t ) ] z [ τ ( t ) ] τ ( t ) ( z [ τ ( t ) ] ) 2 α . (2.20)

Therefore, by (2.18), (2.19), and (2.20), we see that

ω ( t ) ρ ( t ) ρ ( t ) ω ( t ) + ρ ( t ) ( r ( t ) ( z ( t ) ) α ) ( z [ τ ( t ) ] ) α α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ω ( α + 1 ) / α ( t ) . (2.21)

Similarly, we introduce another Riccati substitution:

υ ( t ) : = ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ( z [ τ ( t ) ] ) α , t t 2 . (2.22)

Then υ ( t ) > 0 . Differentiating (2.22), we have

υ ( t ) = ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ( z [ τ ( t ) ] ) α + ρ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ τ ( t ) ] ) α α ρ ( t ) r [ τ ( t ) ] ( z [ τ ( t ) ] ) α z α 1 [ τ ( t ) ] z [ τ ( t ) ] τ ( t ) ( z [ τ ( t ) ] ) 2 α . (2.23)

Therefore, by (2.22) and (2.23), we get

υ ( t ) = ρ ( t ) ρ ( t ) υ ( t ) + ρ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ τ ( t ) ] ) α α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] υ ( α + 1 ) / α ( t ) . (2.24)

Combining (2.21) and (2.24), we have

ω ( t ) + p α [ h ( t ) ] τ ( t ) υ ( t ) ρ ( t ) ( r ( t ) ( z ( t ) ) α ) + p α [ h ( t ) ] τ ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) ( z [ τ ( t ) ] ) α + ρ ( t ) ρ ( t ) ω ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ω ( α + 1 ) / α ( t ) + p α [ h ( t ) ] τ ( t ) ρ ( t ) ρ ( t ) υ ( t ) p α [ h ( t ) ] τ ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] υ ( α + 1 ) / α ( t ) .

It follows from (2.4) and g ( t , a ) τ ( t ) that

ω ( t ) + p α [ h ( t ) ] τ ( t ) υ ( t ) ρ ( t ) 2 α 1 a b Q ( t , ξ ) d σ ( ξ ) + ρ + ( t ) ρ ( t ) ω ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ω ( α + 1 ) / α ( t ) + p α [ h ( t ) ] τ ( t ) ρ + ( t ) ρ ( t ) υ ( t ) p α [ h ( t ) ] τ ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] υ ( α + 1 ) / α ( t ) .

Integrating the latter inequality from t 2 to t, we obtain

ω ( t ) ω ( t 2 ) + p α [ h ( t ) ] τ ( t ) υ ( t ) p α [ h ( t 2 ) ] τ ( t 2 ) υ ( t 2 ) t 2 t ρ ( s ) 2 α 1 a b Q ( s , ξ ) d σ ( ξ ) d s + t 2 t [ ρ + ( s ) ρ ( s ) ω ( s ) α τ ( s ) ρ 1 / α ( s ) r 1 / α [ τ ( s ) ] ω ( α + 1 ) / α ( s ) ] d s + t 2 t p α [ h ( s ) ] τ ( s ) { [ ρ + ( s ) ρ ( s ) + ϕ ( s ) ] + υ ( s ) α τ ( s ) ρ 1 / α ( s ) r 1 / α [ τ ( s ) ] υ ( α + 1 ) / α ( s ) } d s . (2.25)

Define

A : = [ α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ] α / ( α + 1 ) ω ( t ) and B : = [ α α + 1 ρ + ( t ) ρ ( t ) [ α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ] α / ( α + 1 ) ] α .

Using inequality (2.14), we have

ρ + ( t ) ρ ( t ) ω ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ω ( α + 1 ) / α ( t ) 1 ( α + 1 ) α + 1 r [ τ ( t ) ] ( ρ + ( t ) ) α + 1 ( ρ ( t ) τ ( t ) ) α .

On the other hand, define

A : = [ α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ] α / ( α + 1 ) υ ( t ) and B : = [ α α + 1 ζ + ( t ) [ α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] ] α / ( α + 1 ) ] α .

Then, by (2.14), we obtain

ζ + ( t ) υ ( t ) α τ ( t ) ρ 1 / α ( t ) r 1 / α [ τ ( t ) ] υ ( α + 1 ) / α ( t ) 1 ( α + 1 ) α + 1 r [ τ ( t ) ] ( ζ + ( t ) ) α + 1 ρ ( t ) ( τ ( t ) ) α .

Thus, from (2.25), we get

ω ( t ) ω ( t 2 ) + p α [ h ( t ) ] τ ( t ) υ ( t ) p α [ h ( t 2 ) ] τ ( t 2 ) υ ( t 2 ) t 2 t ρ ( s ) { a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 r [ τ ( s ) ] ( α + 1 ) α + 1 ( τ ( s ) ) α × [ ( ρ + ( s ) ρ ( s ) ) α + 1 + p α [ h ( s ) ] ( ζ + ( s ) ) α + 1 τ ( s ) ] } d s ,

which contradicts (2.17). This completes the proof. □

Assuming we have (1.2), where p 0 and τ 0 are constants, we get the following result.

Theorem 2.4Suppose we have (H1)-(H5), (1.2), (1.3), and let τ ( t ) t and g ( t , a ) τ ( t ) for t I . If there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that

lim sup t t 0 t [ ρ ( s ) a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 1 ( α + 1 ) α + 1 ( 1 + p 0 α τ 0 ) r [ τ ( s ) ] ( ρ + ( s ) ) α + 1 ( τ 0 ρ ( s ) ) α ] d s = , (2.26)

then (1.1) is oscillatory.

Proof Assume again that x is an eventually positive solution of (1.1). As in the proof of Theorem 2.1, we have z ( t ) > 0 , (2.3), and (2.4) for all sufficiently large t. By virtue of (1.2), (2.3), and (2.4), we have (2.16) for all sufficiently large t. The rest of the proof is similar to that of Theorem 2.3, and so it is omitted. □

In the following, we present some oscillation criteria for (1.1) in the case where (1.4) holds.

Theorem 2.5Suppose we have (H1)-(H5), (1.2), (1.4), and let g ( t , a ) C 1 ( I , R ) , g ( t , a ) > 0 , g ( t , a ) τ ( t ) t for t I , and g ( t , ξ ) g ( t , b ) for ξ [ a , b ] . Assume further that there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that (2.15) is satisfied. If there exists a real-valued function η C 1 ( I , R ) such that η ( t ) t , η ( t ) g ( t , b ) , η ( t ) > 0 for t I , and

lim sup t t 0 t [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 δ α ( s ) ( 1 + p 0 α τ 0 ) ( α α + 1 ) α + 1 η ( s ) δ ( s ) r 1 / α [ η ( s ) ] ] d s = , (2.27)

then (1.1) is oscillatory.

Proof Let x be a nonoscillatory solution of (1.1). Without loss of generality, we assume that there exists a t 1 I such that x ( t ) > 0 , x [ τ ( t ) ] > 0 , and x [ g ( t , ξ ) ] > 0 for all t t 1 and ξ [ a , b ] . Then z ( t ) > 0 . As in the proof of Theorem 2.1, we get (2.2). By virtue of (1.1), we have (2.3). Thus, r | z | α 1 z is nonincreasing. Now we have two possible cases for the sign of z : (i) z < 0 eventually, or (ii) z > 0 eventually.

(i) Suppose that z ( t ) > 0 for t t 2 t 1 . Then, by the proof of Theorem 2.2, we obtain a contradiction to (2.15).

(ii) Suppose that z ( t ) < 0 for t t 2 t 1 . It follows from (2.2), (2.3), and g ( t , ξ ) g ( t , b ) that

( r ( t ) ( z ( t ) ) α ) + p 0 α τ 0 ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) + 1 2 α 1 z α [ g ( t , b ) ] a b Q ( t , ξ ) d σ ( ξ ) 0 . (2.28)

We define the function u by

u ( t ) : = r ( t ) ( z ( t ) ) α z α [ η ( t ) ] , t t 2 . (2.29)

Then u ( t ) < 0 . Noting that r ( z ) α is nondecreasing, we get

z ( s ) r 1 / α ( t ) r 1 / α ( s ) z ( t ) , s t t 2 .

Integrating this inequality from η ( t ) to l, we obtain

z ( l ) z [ η ( t ) ] + r 1 / α ( t ) z ( t ) η ( t ) l d s r 1 / α ( s ) .

Letting l , we have

0 z [ η ( t ) ] + r 1 / α ( t ) z ( t ) δ ( t ) .

That is,

δ ( t ) r 1 / α ( t ) z ( t ) z [ η ( t ) ] 1 .

Thus, we get by (2.29)

δ α ( t ) u ( t ) 1 . (2.30)

Similarly, we define another function v by

v ( t ) : = r [ τ ( t ) ] ( z [ τ ( t ) ] ) α z α [ η ( t ) ] , t t 2 . (2.31)

Then v ( t ) < 0 . Noting that r ( z ) α is nondecreasing and τ ( t ) t , we get

r ( t ) ( z ( t ) ) α r [ τ ( t ) ] ( z [ τ ( t ) ] ) α .

Thus, 0 < v ( t ) u ( t ) . Hence, by (2.30), we see that

δ α ( t ) v ( t ) 1 . (2.32)

Differentiating (2.29), we obtain

u ( t ) = ( r ( t ) ( z ( t ) ) α ) z α [ η ( t ) ] + α r ( t ) ( z ( t ) ) α z α 1 [ η ( t ) ] z [ η ( t ) ] η ( t ) z 2 α [ η ( t ) ] .

By (2.3) and η ( t ) t , we have z [ η ( t ) ] ( r ( t ) / r [ η ( t ) ] ) 1 / α z ( t ) , and so

u ( t ) ( r ( t ) ( z ( t ) ) α ) z α [ η ( t ) ] α η ( t ) r 1 / α [ η ( t ) ] ( u ( t ) ) ( α + 1 ) / α . (2.33)

Similarly, we see that

v ( t ) ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) z α [ η ( t ) ] α η ( t ) r 1 / α [ η ( t ) ] ( v ( t ) ) ( α + 1 ) / α . (2.34)

Combining (2.33) and (2.34), we get

u ( t ) + p 0 α τ 0 v ( t ) ( r ( t ) ( z ( t ) ) α ) z α [ η ( t ) ] + p 0 α τ 0 ( r [ τ ( t ) ] ( z [ τ ( t ) ] ) α ) z α [ η ( t ) ] α η ( t ) r 1 / α [ η ( t ) ] ( u ( t ) ) ( α + 1 ) / α α p 0 α τ 0 η ( t ) r 1 / α [ η ( t ) ] ( v ( t ) ) ( α + 1 ) / α . (2.35)

Using (2.28), (2.35), and g ( t , b ) η ( t ) , we obtain

u ( t ) + p 0 α τ 0 v ( t ) a b Q ( t , ξ ) d σ ( ξ ) 2 α 1 α η ( t ) r 1 / α [ η ( t ) ] ( u ( t ) ) ( α + 1 ) / α α p 0 α τ 0 η ( t ) r 1 / α [ η ( t ) ] ( v ( t ) ) ( α + 1 ) / α . (2.36)

Multiplying (2.36) by δ α ( t ) and integrating the resulting inequality from t 2 to t, we have

u ( t ) δ α ( t ) u ( t 2 ) δ α ( t 2 ) + α t 2 t δ α 1 ( s ) η ( s ) u ( s ) r 1 / α [ η ( s ) ] d s + α t 2 t η ( s ) δ α ( s ) r 1 / α [ η ( s ) ] ( u ( s ) ) ( α + 1 ) / α d s + p 0 α τ 0 v ( t ) δ α ( t ) p 0 α τ 0 v ( t 2 ) δ α ( t 2 ) + α p 0 α τ 0 t 2 t δ α 1 ( s ) η ( s ) v ( s ) r 1 / α [ η ( s ) ] d s + α p 0 α τ 0 t 2 t η ( s ) δ α ( s ) r 1 / α [ η ( s ) ] ( v ( s ) ) ( α + 1 ) / α d s + t 2 t a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 δ α ( s ) d s 0 .

Set

A : = [ η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ] ( α + 1 ) / α u ( t ) and B : = [ α α + 1 δ α 1 ( t ) η ( t ) r 1 / α [ η ( t ) ] [ η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ] α / ( α + 1 ) ] α .

Using inequality (2.14), we get

δ α 1 ( t ) η ( t ) u ( t ) r 1 / α [ η ( t ) ] + η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ( u ( t ) ) ( α + 1 ) / α 1 α ( α α + 1 ) α + 1 η ( t ) δ ( t ) r 1 / α [ η ( t ) ] .

Similarly, we set

A : = [ η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ] ( α + 1 ) / α v ( t ) and B : = [ α α + 1 δ α 1 ( t ) η ( t ) r 1 / α [ η ( t ) ] [ η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ] α / ( α + 1 ) ] α .

Then we have by (2.14)

δ α 1 ( t ) η ( t ) v ( t ) r 1 / α [ η ( t ) ] + η ( t ) δ α ( t ) r 1 / α [ η ( t ) ] ( v ( t ) ) ( α + 1 ) / α 1 α ( α α + 1 ) α + 1 η ( t ) δ ( t ) r 1 / α [ η ( t ) ] .

Thus, from (2.30) and (2.32), we find

t 2 t [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 δ α ( s ) ( 1 + p 0 α τ 0 ) ( α α + 1 ) α + 1 η ( s ) δ ( s ) r 1 / α [ η ( s ) ] ] d s u ( t 2 ) δ α ( t 2 ) + p 0 α τ 0 v ( t 2 ) δ α ( t 2 ) + 1 + p 0 α τ 0 ,

which contradicts (2.27). This completes the proof. □

With a proof similar to the proof of Theorems 2.4 and 2.5, we obtain the following result.

Theorem 2.6Suppose we have (H1)-(H5), (1.2), (1.4), and let τ ( t ) t , g ( t , a ) τ ( t ) for t I , and g ( t , ξ ) g ( t , b ) for ξ [ a , b ] . Assume also that there exists a real-valued function ρ C 1 ( I , ( 0 , ) ) such that (2.26) is satisfied. If there exists a real-valued function η C 1 ( I , R ) such that η ( t ) t , η ( t ) g ( t , b ) , η ( t ) > 0 for t I , and (2.27) holds, then (1.1) is oscillatory.

3 Applications and discussion

In this section, we provide three examples to illustrate the main results.

Example 3.1 Consider the second-order neutral functional differential equation

[ x ( t ) + x ( t 2 π ) ] + 5 π 2 π 2 x [ t + ξ ] d ξ = 0 , t 10 . (3.1)

Let α = 1 , a = 5 π / 2 , b = π / 2 , r ( t ) = 1 , p ( t ) = 1 , τ ( t ) = t 2 π , q ( t , ξ ) = 1 , g ( t , ξ ) = t + ξ , σ ( ξ ) = ξ , and ρ ( t ) = 1 . Then Q ( t , ξ ) = min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } = 1 , g ( t , a ) = 1 , g ( t , a ) = t 5 π / 2 t + ξ for ξ [ 5 π / 2 , π / 2 ] , and g ( t , a ) τ ( t ) t . Moreover, letting τ 0 = 1 , then

lim sup t t 0 t [ ρ ( s ) a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 1 ( α + 1 ) α + 1 ( 1 + p 0 α τ 0 ) r [ g ( s , a ) ] ( ρ + ( s ) ) α + 1 ( ρ ( s ) g ( s , a ) ) α ] d s = 3 π lim sup t 10 t d s = .

Hence, by Theorem 2.2, (3.1) is oscillatory. As a matter of fact, one such solution is x ( t ) = sin t .

Example 3.2 Consider the second-order neutral functional differential equation

[ x ( t ) + t x ( t β ) ] + 0 1 ξ + 1 t x [ t + ξ ] d ξ = 0 , t 1 , (3.2)

where β 0 is a constant. Let α = 1 , a = 0 , b = 1 , r ( t ) = 1 , p ( t ) = t , τ ( t ) = t β , q ( t , ξ ) = ( ξ + 1 ) / t , g ( t , ξ ) = t + ξ , σ ( ξ ) = ξ , and ρ ( t ) = 1 . Then Q ( t , ξ ) = min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } = ( ξ + 1 ) / t , g ( t , a ) = g ( t , 0 ) = t t + ξ for ξ [ 0 , 1 ] , τ ( t ) = t β t , and g ( t , a ) τ ( t ) for t 1 . Further, setting h ( t ) = t + 1 ,

ϕ ( t ) = α p [ h ( t ) ] h ( t ) p [ h ( t ) ] τ ( t ) τ ( t ) = 1 t + 1 , ζ ( t ) = ρ + ( t ) ρ ( t ) + ϕ ( t ) = 1 t + 1 ,

and

φ ( t ) = ( ρ + ( t ) ρ ( t ) ) α + 1 + p α [ h ( t ) ] ( ζ + ( t ) ) α + 1 τ ( t ) = 1 t + 1 .

Therefore, we have

lim sup t t 0 t ρ ( s ) [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 r [ τ ( s ) ] φ ( s ) ( α + 1 ) α + 1 ( τ ( s ) ) α ] d s = lim sup t 1 t [ 0 1 ξ + 1 s d ξ 1 4 ( s + 1 ) ] d s = lim sup t 1 t [ 3 2 s 1 4 ( s + 1 ) ] d s = .

Hence, (3.2) is oscillatory due to Theorem 2.3.

Example 3.3 Consider the second-order neutral functional differential equation

[ t 2 ( x ( t ) + p ( t ) x ( t β ) ) ] + 0 1 ( ξ + 1 ) x [ t + ξ ] d ξ = 0 , t 1 , (3.3)

where 0 p ( t ) p 0 , p 0 and β are positive constants. Let α = 1 , a = 0 , b = 1 , r ( t ) = t 2 , τ ( t ) = t β , q ( t , ξ ) = ξ + 1 , g ( t , ξ ) = t + ξ , σ ( ξ ) = ξ , ρ ( t ) = 1 , and η ( t ) = t + 1 . Then Q ( t , ξ ) = min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } = ξ + 1 , τ 0 = 1 , g ( t , a ) = g ( t , 0 ) = t t + ξ for ξ [ 0 , 1 ] , τ ( t ) = t β t , g ( t , a ) τ ( t ) for t 1 , and δ ( t ) = 1 / t . Further,

lim sup t t 0 t [ ρ ( s ) a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 1 ( α + 1 ) α + 1 ( 1 + p 0 α τ 0 ) r [ τ ( s ) ] ( ρ + ( s ) ) α + 1 ( τ 0 ρ ( s ) ) α ] d s = 3 2 lim sup t 1 t d s =

and

lim sup t t 0 t [ a b Q ( s , ξ ) d σ ( ξ ) 2 α 1 δ α ( s ) ( 1 + p 0 α τ 0 ) ( α α + 1 ) α + 1 η ( s ) δ ( s ) r 1 / α [ η ( s ) ] ] d s = ( 3 2 1 + p 0 4 ) lim sup t 1 t d s s + 1 = , if  p 0 < 5 .

Hence, by Theorem 2.6, (3.3) is oscillatory when 0 p ( t ) p 0 < 5 .

Remark 3.1 In this paper, we establish some new oscillation theorems for (1.1) in the case where p is finite or infinite on . The criteria obtained extend the results in [22] and improve those reported in [19]. Similar results can be presented under the assumption that 0 < α 1 . In this case, using [[5], Lemma 2], one has to replace Q ( t , ξ ) : = min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } with Q ( t , ξ ) : = 2 α 1 min { q ( t , ξ ) , q ( τ ( t ) , ξ ) } and proceed as above. It would be interesting to find another method to investigate (1.1) in the case where g ( τ ( t ) , ξ ) τ [ g ( t , ξ ) ] .

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to this work. They all read and approved the final version of the manuscript.

Acknowledgements

The authors express their sincere gratitude to the anonymous referees for the careful reading of the original manuscript and useful comments that helped to improve the presentation of the results and accentuate important details.

References

  1. Hale, JK: Theory of Functional Differential Equations, Springer, New York (1977)

  2. Agarwal, RP, Bohner, M, Li, W-T: Nonoscillation and Oscillation Theory for Functional Differential Equations, Dekker, New York (2004)

  3. Agarwal, RP, Grace, SR, O’Regan, D: Oscillation Theory for Difference and Functional Differential Equations, Kluwer Academic, Dordrecht (2000)

  4. Baculíková, B, Džurina, J: Oscillation theorems for second order neutral differential equations. Comput. Math. Appl.. 61, 94–99 (2011). Publisher Full Text OpenURL

  5. Baculíková, B, Džurina, J: Oscillation theorems for second order nonlinear neutral differential equations. Comput. Math. Appl.. 62, 4472–4478 (2011). Publisher Full Text OpenURL

  6. Baculíková, B, Lacková, D: Oscillation criteria for second order retarded differential equations. Stud. Univ. Zilina Math. Ser.. 20, 11–18 (2006)

  7. Baculíková, B, Li, T, Džurina, J: Oscillation theorems for second-order superlinear neutral differential equations. Math. Slovaca. 63, 123–134 (2013). Publisher Full Text OpenURL

  8. Candan, T: The existence of nonoscillatory solutions of higher order nonlinear neutral equations. Appl. Math. Lett.. 25, 412–416 (2012). Publisher Full Text OpenURL

  9. Candan, T, Dahiya, RS: Existence of nonoscillatory solutions of higher order neutral differential equations with distributed deviating arguments. Math. Slovaca. 63, 183–190 (2013). Publisher Full Text OpenURL

  10. Candan, T, Karpuz, B, Öcalan, Ö: Oscillation of neutral differential equations with distributed deviating arguments. Nonlinear Oscil.. 15, 65–76 (2012)

  11. Dix, JG, Karpuz, B, Rath, R: Necessary and sufficient conditions for the oscillation of differential equations involving distributed arguments. Electron. J. Qual. Theory Differ. Equ.. 2011, 1–15 (2011)

  12. Džurina, J, Hudáková, D: Oscillation of second order neutral delay differential equations. Math. Bohem.. 134, 31–38 (2009)

  13. Hasanbulli, M, Rogovchenko, YuV: Oscillation criteria for second order nonlinear neutral differential equations. Appl. Math. Comput.. 215, 4392–4399 (2010). Publisher Full Text OpenURL

  14. Karpuz, B, Öcalan, Ö, Öztürk, S: Comparison theorems on the oscillation and asymptotic behaviour of higher-order neutral differential equations. Glasg. Math. J.. 52, 107–114 (2010). Publisher Full Text OpenURL

  15. Li, T, Agarwal, RP, Bohner, M: Some oscillation results for second-order neutral differential equations. J. Indian Math. Soc.. 79, 97–106 (2012)

  16. Li, T, Han, Z, Zhang, C, Li, H: Oscillation criteria for second-order superlinear neutral differential equations. Abstr. Appl. Anal.. 2011, Article ID 367541 (2011)

  17. Li, T, Rogovchenko, YuV, Zhang, C: Oscillation of second-order neutral differential equations. Funkc. Ekvacioj. 56, 111–120 (2013)

  18. Li, T, Sun, S, Han, Z, Han, B, Sun, Y: Oscillation results for second-order quasi-linear neutral delay differential equations. Hacet. J. Math. Stat.. 42, 131–138 (2013)

  19. Li, T, Thandapani, E: Oscillation of second-order quasi-linear neutral functional dynamic equations with distributed deviating arguments. J. Nonlinear Sci. Appl.. 4, 180–192 (2011)

  20. Li, WN: Oscillation of higher order delay differential equations of neutral type. Georgian Math. J.. 7, 347–353 (2000)

  21. Şenel, MT, Candan, T: Oscillation of second order nonlinear neutral differential equation. J. Comput. Anal. Appl.. 14, 1112–1117 (2012)

  22. Sun, S, Li, T, Han, Z, Li, H: Oscillation theorems for second-order quasilinear neutral functional differential equations. Abstr. Appl. Anal.. 2012, Article ID 819342 (2012)

  23. Thandapani, E, Piramanantham, V: Oscillation criteria of second order neutral delay dynamic equations with distributed deviating arguments. Electron. J. Qual. Theory Differ. Equ.. 2010, 1–15 (2010)

  24. Wang, P: Oscillation criteria for second order neutral equations with distributed deviating arguments. Comput. Math. Appl.. 47, 1935–1946 (2004). Publisher Full Text OpenURL

  25. Xu, Z, Weng, P: Oscillation of second-order neutral equations with distributed deviating arguments. J. Comput. Appl. Math.. 202, 460–477 (2007). Publisher Full Text OpenURL

  26. Ye, L, Xu, Z: Interval oscillation of second order neutral equations with distributed deviating argument. Adv. Dyn. Syst. Appl.. 1, 219–233 (2006)

  27. Yu, Y, Fu, X: Oscillation of second order neutral equation with continuous distributed deviating argument. Rad. Mat.. 7, 167–176 (1991)

  28. Zhang, C, Agarwal, RP, Bohner, M, Li, T: Oscillation of second-order nonlinear neutral dynamic equations with noncanonical operators. Bull. Malays. Math. Sci. Soc. (2014, in press)

  29. Zhang, C, Baculíková, B, Džurina, J, Li, T: Oscillation results for second-order mixed neutral differential equations with distributed deviating arguments. Math. Slovaca (2014, in press)

  30. Zhao, J, Meng, F: Oscillation criteria for second-order neutral equations with distributed deviating argument. Appl. Math. Comput.. 206, 485–493 (2008). Publisher Full Text OpenURL