AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 787920 10.1155/2012/787920 787920 Research Article On the Zeroes and the Critical Points of a Solution of a Second Order Half-Linear Differential Equation Almenar Pedro 1 Jódar Lucas 2 Atici Ferhan 1 Division of Network, Vodafone Spain S. A. P. E. Castellana Norte, 28050 Madrid Spain vodafone.es 2 Instituto Universitario de Matemática Multidisciplinar Universitat Politècnica de València Camino de Vera s/n, 46022 Valencia Spain upv.es 2012 28 12 2012 2012 11 09 2012 12 11 2012 28 11 2012 2012 Copyright © 2012 Pedro Almenar and Lucas Jódar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents two methods to obtain upper bounds for the distance between a zero and an adjacent critical point of a solution of the second-order half-linear differential equation (p(x)Φ(y'))'+q(x)Φ(y)=0, with p(x) and q(x) piecewise continuous and p(x)>0, Φ(t)=|t|r-2t and r being real such that r>1. It also compares between them in several examples. Lower bounds (i.e., Lyapunov inequalities) for such a distance are also provided and compared with other methods.

1. Introduction

In a recent paper of the authors (see ), a method to calculate upper bounds for the distance between a zero and an adjacent critical point of a solution of the linear equation (1.1)(p(x)y)+q(x)y=0,x>x0, was introduced. The purpose of this paper is to extend the results described there to the half-linear differential equation (1.2)(p(x)Φ(y))+q(x)Φ(y)=0,x>x0, where p(x),q(x)>0, Φ(t)=|t|r-2t and r is a real number such that r>1.

As happened in the linear case and was commented in , the existing literature on the analysis of the distribution of the zeroes and the critical points of the half-linear equation (1.2) has focused on the determination of the oscillation or nonoscillation nature of the solutions of (1.2) (see, e.g., [2, Sections 5 and 9] for a compilation of such works), on the derivation of formulae for the asymptotic distance between consecutive zeroes (see ), and on the calculation of lower bounds for the distance between consecutive zeroes and between zeroes and adjacent critical points, that is, on extensions of the Lyapunov inequality to the half-linear case ([2, 49] are good examples of that). In contrast, it is really difficult (if not impossible) to find any papers that provide insights on how big the distance between consecutive zeroes or between a zero and a critical point (or, in other words, between a critical point and its adjacent focal points) may be for the half-linear equation. The present paper aims at solving that by providing two alternative bounding methods and comparing weaknesses and advantages of each. Since the results obtained here also allow us to define methods to obtain lower bounds for the aforementioned distance, a comparison of these methods and those of [2, 69] will also be provided.

Following [2, Section 1.2], we will denote by S(x) the solution of the half-linear differential equation (1.3)(Φ(S))+(r-1)Φ(S)=0, given by the initial conditions S(0)=0 and S'(0)=1. If we define (1.4)πr=2πrsin(π/r), it is possible to show (see again [2, Section 1.2]) that the behaviour of this solution is very similar to that of the classical sine function, in the sense that (1.5)|S(x)|r+|S(x)|r=|S(0)|r+|S(0)|r=1,x[0,πr2],S(x) is increasing on ]0,πr/2[ and S(πr/2)=1 (resp., S'(x) decreasing on ]0,πr/2[ and S'(πr/2)=0, according to (1.5)).

This resemblance increases for the generalized sine function Sr(x) defined in the whole real line as the 2πr-periodic odd extension of the function given by (1.6)Sr(x)=S(x),x[0,πr2],Sr(x)=S(πr-x),x[πr2,πr], which obviously satisfies (1.7)|Sr(x)|r+|Sr(x)|r=1,x,|Sr(x)| being increasing on ]nπr,((2n+1)/2)πr[ and decreasing on ]((2n-1)/2)πr,nπr[ for any integer n, with |Sr(((2n+1)/2)πr)|=1 and Sr(nπr)=0 (resp., |Sr(x)| being increasing on ]((2n-1)/2)πr,nπr[ and decreasing on ]nπr,((2n+1)/2)πr[ for any integer n, with Sr(((2n+1)/2)πr)=0 and |Sr(nπr)|=1).

In fact for r=2,Sr(x) is the classical sine function sinx. It is straightforward to show that Sr(x) also satisfies (1.3) for the whole real line.

With this in mind, we introduce the half-linear cotangent function cotr(x) as (1.8)cotrx=Sr(x)Sr(x),x]0,πr[. Throughout the paper s will be the conjugate number of r, that is, (1/r)+(1/s)=1. It is straightforward to show that (1.9)(s-1)(r-1)=1. The organization of the paper is as follows. Section 2 will prove the main results. Section 3 will apply them to devise methods to calculate upper and lower bounds for the distance between a zero and the critical point immediately preceding or following it. Section 4 will provide some examples to illustrate advantages and disadvantages of the methods presented here for the calculation of upper bounds. Section 5 will compare the lower bounds obtained from the main results with other bounds existing in the literature. Finally, Section 6 will state several conclusions.

2. Main Results

The following extended mean value theorem for integrals condenses Theorems 3 and 4 of  and will be the key for the results of this paper, as happened in .

Theorem 2.1 (extended mean value theorem for integrals).

Let a(x), f(x), m(x) be piecewise continuous functions on [a,b] with f(x), and m(x)>0 on ]a,b[. Let K be defined by (2.1)K=aba(x)f(x)dxabf(x)dx. Then, if f(x)/m(x) is increasing on ]a,b[, one has (2.2)abmin{a(t);xtb}m(x)dxabm(x)dxKabmax{a(t);xtb}m(x)dxabm(x)dx. If f(x)/m(x) is decreasing on ]a,b[, one has (2.3)abmin{a(t);atx}m(x)dxabm(x)dxKabmax{a(t);atx}m(x)dxabm(x)dx.

Before replicating the theorems of , we will need to define a similar Prüfer transformation that allows us to obtain a first order differential equation for a sort of “angle” function constructed with the solution of (1.2) and its derivative and with the aid of the cotangent function defined in (1.8). That is the purpose of the next lemma.

Lemma 2.2.

Let y(x) be a nonidentically zero solution of (1.2) with p(x) and q(x) piecewise continuous and positive on an interval I. Let A>0 be a real number. Let R(x) and φ(x) be defined by the Prüfer transformation (2.4)R(x)=|y(x)|r+ps(x)|y(x)|rA2rr, cot rφ(x)=(ps-1(x)y(x)A2y(x)). Then, one has (2.5)φ(x)=p1-s(x)|Sr(φ(x))|rA2+q(x)(r-1)A2(r-1)|Sr(φ(x))|r.

Proof.

Dividing (1.2) by A2(r-1), one gets the equivalent equation (2.6)(p(x)A2(r-1)Φ(y))+q(x)A2(r-1)Φ(y)=0,x>x0. From [2, Section 1.2], (1.9), and (2.6), it is straightforward to show that the angle function φ(x) defined in (2.4) satisfies (2.5).

With the aid of Theorem 2.1 and Lemma 2.2, it is possible to prove the next theorem.

Theorem 2.3.

Let y(x) be a non-identically zero solution of (1.2) with p(x) and q(x) piecewise continuous and positive on an interval I. Let a, b be real numbers such that [a,b]I. Suppose that m(x) is any positive piecewise continuous function on [a,b].

Then, if y'(a)=y(b)=0, one has (2.7)(abmin{1p1/(r-1)(t)m(t),xtb}m(x)dx)1/s·(abmin{q(t)m(t),atx}m(x)dx)1/rπr2r-1r,(2.8)(abmax{1p1/(r-1)(t)m(t),xtb}m(x)dx)1/s·(abmax{q(t)m(t),atx}m(x)dx)1/rπr2r-1r. If y(a)=y'(b)=0, one has (2.9)(abmin{1p1/(r-1)(t)m(t),atx}m(x)dx)1/s·(abmin{q(t)m(t),xtb}m(x)dx)1/rπr2r-1r,(2.10)(abmax{1p1/(r-1)(t)m(t),atx}m(x)dx)1/s·(abmax{q(t)m(t),xtb}m(x)dx)1/rπr2r-1r.

Proof.

Let us focus first on proving (2.7) for the case y(a)=y(b)=0. From Lemma 2.2, it is clear that the Prüfer transformation defined by (2.4) for any A>0 satisfies the first-order differential equation (2.5).

Since p(x),q(x)>0 on I, the right-hand side of the equation is positive. Therefore, φ(x) must be a nondecreasing function with φ(a)=((2n-1)/2)πr and φ(b)=nπr for an integer n. That implies that |Sr(φ(x))| must be decreasing and |Sr(φ(x))| increasing on ]a,b[.

Integrating (2.5) between a and b, one has (2.11)πr2=φ(b)-φ(a)=abφ(x)dx=abp1-s(x)|Sr(φ(x))|rA2dx+abq(x)(r-1)A2(r-1)|Sr(φ(x))|rdx. Since |Sr(φ(x))| is increasing on [a,b], we can apply (2.2) to yield (2.12)abp1-s(x)|Sr(φ(x))|rA2dx=abp1-s(x)|Sr(φ(x))|rA2dxabm(x)|Sr(φ(x))|rdxabm(x)|Sr(φ(x))|rdxA2abmin{(1/ps-1(t)m(t)),xtb}m(x)dxabm(x)dxabm(x)|Sr(φ(x))|rdx. Likewise, we can apply (2.3) to yield (2.13)abq(x)(r-1)A2(r-1)|Sr(φ(x))|rdx=ab(q(x)|Sr(φ(x))|rdx/(r-1)A2(r-1))abm(x)|Sr(φ(x))|rdxabm(x)|Sr(φ(x))|rdx1(r-1)A2(r-1)abmin{q(t)/m(t),atx}m(x)dxabm(x)dxabm(x)|Sr(φ(x))|rdx. Let us take now (2.14)A=(abmin{q(t)/m(t),atx}m(x)dx(r-1)abmin{1/ps-1(t)m(t),xtb}m(x)dx)1/2r. From (2.12) and (2.14), one has (2.15)abp1-s(x)|Sr(φ(x))|rA2dx(abmin{1p1/(r-1)(t)m(t),xtb}m(x)dx)1/s·(abmin{q(t)/m(t),atx}m(x)dxr-1)1/rabm(x)|Sr(φ(x))|rdxabm(x)dx. Likewise, from (2.13) and (2.14), one has (2.16)abq(x)(r-1)A2(r-1)|Sr(φ(x))|rdx(abmin{1p1/(r-1)(t)m(t),xtb}m(x)dx)1/s·(abmin{q(t)/m(t),atx}m(x)dxr-1)1/rabm(x)|Sr(φ(x))|rdxabm(x)dx. From (1.7), (2.11), (2.15), and (2.16), one gets (2.7).

The proof of (2.8)–(2.10) is similar and will not be repeated here.

Remark 2.4.

Theorem 2 of  is a particular case of Theorem 2.3, when r=2.

Remark 2.5.

The conditions p(x),q(x)>0 on [a,b] of Theorem 2.3 are only used to guarantee in (2.5) that φ(x) is increasing on [a,b]. These constraints can be relaxed if it can be proved that their violation will not force φ(x) to cease being increasing on the interval of interest (e.g., because the negativity of either function happens very close to one of the extrema a or b of the interval where the associated term with φ(x) is close to zero).

Setting m(x)=1 and m(x)=q(x)/p1/(r-1)(x), respectively, in Theorem 2.3, it is straightforward to obtain the following corollaries, which will be used in the examples of Section 4.

Corollary 2.6.

Let y(x) be a non-identically zero solution of (1.2) with p(x) and q(x) piecewise continuous and positive on an interval I. Let a, b be real numbers such that [a,b]I.

Then, if y'(a)=y(b)=0, one has (2.17)(abmin{1p1/(r-1)(t),xtb}dx)1/s·(abmin{q(t),atx}dx)1/rπr2r-1r,(abmax{1p1/(r-1)(t),xtb}dx)1/s·(abmax{q(t),atx}dx)1/rπr2r-1r. If y(a)=y'(b)=0, one has (2.18)(abmin{1p1/(r-1)(t),atx}dx)1/s·(abmin{q(t),xtb}dx)1/rπr2r-1r,(abmax{1p1/(r-1)(t),atx}dx)1/s·(abmax{q(t),xtb}dx)1/rπr2r-1r.

Corollary 2.7.

Let y(x) be a non-identically zero solution of (1.2) with p(x) and q(x) piecewise continuous and positive on an interval I. Let a, b be real numbers such that [a,b]I.

Then, if y'(a)=y(b)=0, one has (2.19)(abmin{1p1/(r-1)(t)q(t),xtb}q(x)p1/(r-1)(x)dx)1/s·(abmin{p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/rπr2r-1r,(2.20)(abmax{1p1/(r-1)(t)q(t),xtb}q(x)p1/(r-1)(x)dx)1/s·(abmax{p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/rπr2r-1r. If y(a)=y'(b)=0, one has (2.21)(abmin{1p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/s·(abmin{p1/(r-1)(t)q(t),xtb}q(x)p1/(r-1)(x)dx)1/rπr2r-1r,(2.22)(abmax{1p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/s·(abmax{p1/(r-1)(t)q(t),xtb}q(x)p1/(r-1)(x)dx)1/rπr2r-1r.

The following corollary allows to slightly simplify the expressions on inequalities (2.7) and (2.9).

Corollary 2.8.

Let y(x) be a non-identically zero solution of (1.2) with p(x) and q(x) piecewise continuous and positive on an interval I. Let a, b be real numbers such that [a,b]I. Suppose that m(x) is any positive piecewise continuous function on [a,b].

Then, if y'(a)=y(b)=0, one has (2.23)abmin{1m(t)p1/(r-1)(t),xtb}smin{q(t)m(t),atx}rm(x)dxπr2r-1r. If y(a)=y'(b)=0, one has (2.24)abmin{1m(t)p1/(r-1)(t),atx}smin{q(t)m(t),xtb}rm(x)dxπr2r-1r.

Proof.

From Hölder inequality, it is straightforward to prove that (2.25)abmin{1m(t)p1/(r-1)(t),xtb}smin{q(t)m(t),atx}rm(x)dx(abmin{1p1/(r-1)(t)m(t),xtb}m(x)dx)1/s·(abmin{q(t)m(t),atx}m(x)dx)1/r,(2.26)abmin{1m(t)p1/(r-1)(t),atx}smin{q(t)m(t),xtb}rm(x)dx(abmin{1p1/(r-1)(t)m(t),atx}m(x)dx)1/s·(abmin{q(t)m(t),xtb}m(x)dx)1/r. Now, from (2.7) and (2.25), one gets (2.23). Likewise, from (2.9) and (2.26), one gets (2.24).

To complete this section we will use Corollary 2.7 to establish a result on the distance between zeroes and consecutive critical points of the solutions of (1.2) when p1/(r-1)(x)q(x) is monotonic.

Theorem 2.9.

Let p(x), q(x) be positive piecewise continuous functions on an interval I. Let y1(x), y2(x) be non-identically zero solutions of (1.2) such that y1(a)=y2(a)=0; with aI. Let b1 and b2 be defined such that y1(b1)=y2(b2)=0, that is, b1 is the critical point of y1(x) consecutive to a and b2 the zero of y2(x) consecutive to a, and let one suppose that [a,b1] and [a,b2]I.

If p1/(r-1)(x)q(x) is decreasing on [a,max(b1,b2)], then b2b1.

On the contrary, if p1/(r-1)(x)q(x) is increasing on [a,max(b1,b2)], then b2b1.

Proof.

Let us first suppose that p1/(r-1)(x)q(x) is decreasing on [a,max(b1,b2)]. From (2.19), one has (2.27)(ab2dxp1/(r-1)(x))1/s(ab2q(x)dx)1/r=(ab21p1/(r-1)(x)q(x)q(x)p1/(r-1)(x)dx)1/s×(ab2p1/(r-1)(x)q(x)q(x)p1/(r-1)(x)dx)1/r=(ab2min{1p1/(r-1)(t)q(t),xtb2}q(x)p1/(r-1)(x)dx)1/s·(ab2min{p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/rπr2r-1r. From (2.22), one has (2.28)(ab1dxp1/(r-1)(x))1/s(ab1q(x)dx)1/r=(ab11p1/(r-1)(x)q(x)q(x)p1/(r-1)(x)dx)1/s×(ab1p1/(r-1)(x)q(x)q(x)p1/(r-1)(x)dx)1/r=(ab1max{1p1/(r-1)(t)q(t),atx}q(x)p1/(r-1)(x)dx)1/s·(ab1max{p1/(r-1)(t)q(t),xtb1}q(x)p1/(r-1)(x)dx)1/rπr2r-1r. From (2.27) and (2.28) one gets (2.29)(ab2dxp1/(r-1)(x))1/s(ab2q(x)dx)1/r(ab1dxp1/(r-1)(x))1/s(ab1q(x)dx)1/r, which obviously gives b2b1.

The case p1/(r-1)(x)q(x) being increasing on [a,max(b1,b2)] can be easily proved in the same manner.

3. Methods to Obtain Upper and Lower Bounds

The results presented in the previous section allow us to define two methods to calculate upper and lower bounds for the distance between a zero and a critical point of a solution of (1.2).

Thus, on one hand, Lemma 2.2 can be leveraged to derive a method similar to that of Moore (see [11, Theorem 8]) for the linear equation (1.1), as the following theorem shows.

Theorem 3.1.

Let y(x) be a non-identically zero solution of (1.2) with p(x) and q(x) piecewise continuous and p(x)>0 on an interval I. Let a, b be real numbers such that [a,b]I. Let A be any positive real number.

If y(a)=y'(b)=0 or y'(a)=y(b)=0, one has (3.1)abmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}dxπr2,(3.2)abmax{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}dxπr2.

Proof.

Integrating (2.5) from a to b and applying (1.5), one has (3.3)abmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}dx=abmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}|Sr(φ(x))|rdx+abmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}|Sr(φ(x))|rdxabA2p1/(r-1)(x)|Sr(φ(x))|rdx+abq(x)(r-1)A2(r-1)|Sr(φ(x))|rdx=πr2, regardless of a being a zero or a critical point. Equation (3.2) can be obtained in the same manner.

Remark 3.2.

Note that unlike what happens in Theorems 2.3 and 2.9 and Corollaries 2.6, 2.7, and 2.8, Theorem 3.1 allows q(x) to be negative or zero on the interval I, which makes it applicable to a wider set of problems of the form (1.2).

Since the previous theorem is applicable regardless of the concrete value of A, it is evident that if we search for an A that minimizes the value of the extreme b that yields equality in (3.1) for fixed a and A, we will found the best upper bound for b obtainable with this method. Likewise, if we search for the A that maximizes the value of the extreme b that yields equality in (3.2) for fixed a and A, we will found the best lower bound for b obtainable with this method.

On the other hand, Theorem 2.3 and Corollary 2.6 can also be used to establish upper and lower bounds for the distance between a zero and a critical point of a solution of (1.2). This is quite evident in the case of the lower bound, since the left-hand sides of (2.8) and (2.10) are both increasing functions of the extreme b. Therefore, the smallest b that gives an equality on (2.8)—for the case y'(a)=y(b)=0—or on (2.10)—for the case y(a)=y'(b)=0—will be the lower bound for the searched zero or critical point of the solution of (1.2), as any smaller b will not satisfy (2.8) or (2.10), respectively, and therefore cannot be a zero or a critical point.

In the case of upper bounds, the left hand sides of (2.7) and (2.9) are not necessarily increasing with b, regardless of them being products of positive integrals. The reason for that lies on the fact that the integrands of all integrals appearing in such inequalities contain minima of functions on an interval: if the length of the interval grows (as b grows), the minima can become smaller.

However, if one increases b from a starting point a (e.g., a critical point of a solution of (1.2)), until one finds a b* such that (3.4)(ab*min{1p1/(r-1)(t)m(t),xtb*}m(x)dx)1/s·(ab*min{q(t)m(t),atx}m(x)dx)1/r=πr2r-1r, then such a value b* will be an upper bound for the next zero of the solution of (1.2). To prove that let us pick any b>b* and let us select (3.5)A(b*)=(ab*min{q(t)/m(t),atx}m(x)dx(r-1)ab*min{1/ps-1(t)m(t),xtb*}m(x)dx)1/2r, according to (2.14) in Theorem 2.3. From (2.12), (2.13), (3.4), and (3.5), one has (3.6)abA2(b*)p1/(r-1)(x)|Sr(φ(x))|rdx+abq(x)(r-1)A2(r-1)(b*)|Sr(φ(x))|rdx=ab*A2(b*)p1/(r-1)(x)|Sr(φ(x))|rdx+ab*q(x)(r-1)A2(r-1)(b*)|Sr(φ(x))|rdx+b*bA2(b*)p1/(r-1)(x)|Sr(φ(x))|rdx+b*bq(x)(r-1)A2(r-1)(b*)|Sr(φ(x))|rdx>πr2, given that p,q>0 on [a,b]. It is clear that (3.6) violates (2.11), and therefore b cannot be the next zero of y(x). A similar result can be proved for the case a being a zero and b being a critical point.

The previous reasoning gives also a possible approach for determining upper bounds when the left-hand sides of (2.7) or (2.9), depending on the case, find a maximum M lower than (πr/2)r-1r for a value b*. Fixed b* basically consists in determining the lowest b such that (3.7)b*bmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}dx=πr2-Mr-1r, with (3.8)A(b*)=(ab*min{q(t)/m(t),atx}m(x)dx(r-1)ab*min{1/ps-1(t)m(t),xtb*}m(x)dx)1/2r, for the case y'(a)=y(b)=0, and (3.9)A(b*)=(ab*min{q(t)/m(t),xtb*}m(x)dx(r-1)ab*min{1/ps-1(t)m(t),atx}m(x)dx)1/2r, for the case y(a)=y'(b)=0. From (2.12), (2.13), and (3.7)–(3.9), one has (3.10)abA2p1/(r-1)(x)|Sr(φ(x))|rdx+abq(x)(r-1)A2(r-1)|Sr(φ(x))|rdxMr-1r+b*bmin{A2p1/(r-1)(x),q(x)(r-1)A2(r-1)}dx=πr2. It is important to remark, in any case, that this mechanism does not guarantee success in the search of the upper bound b in all cases, as the integral between b* and may not get to be large enough to fulfil the difference between πr/2 and M/r-1r.

4. Some Examples

Throughout this section, we will introduce examples where upper bounds for the distance between a zero and a critical point of a solution of (1.2) will be provided by means of the two methods presented in the previous section: the one based on Theorem 2.3 (or Corollary 2.6, depending on the case) and the extension of Moore's results (see ) to the half-linear equation introduced in Theorem 3.1. All the examples will address the cases r=1.5 and r=5 to cope with half-linear equations whose exponent r is lower and greater (resp.) than the value r of the linear case (r=2). To simplify the comparison, the analysis will fix the value of the starting point a (either a zero or a critical point) and will search for an upper bound of the corresponding point b (next critical point or zero, resp., adjacent to a).

Example 4.1.

Let us consider the following half-linear differential equation: (4.1)Φ(y)+1xΦ(y)=0,x1,a=1. Corollary 2.6 allows us to find an upper bound for b by searching for the lowest value B that satisfies (4.2)(B-1)1/s(1Bdxx)1/rπr2r-1r, in the case y'(1)=y(b)=0, and (4.3)(B-1)1/s(1B(B-1))1/r=B-1B2rπr2r-1r, in the case y(1)=y'(b)=0.

Likewise, Theorem 3.1 allows us to find an upper bound for b by searching for the lowest value B that satisfies (4.4)1Bmin{A2,1(r-1)A2(r-1)x}dxπr2 across all possible values of A.

The results of the mentioned calculations are summarized in Table 1.

As it can be seen in Table 1, for both values r=1.5 and r=5 the bounds obtained using Corollary 2.6 improve those of Theorem 3.1 in the case y'(1)=y(b)=0, and are worse than those provided by Theorem 3.1 in the case y(1)=y'(b)=0.

Comparison of upper bounds for b in Example 4.1.

Corollary 2.6 Theorem 3.1
r = 1.5 y ( 1 ) = y ' ( b ) = 0 b 3.26 b 2.99
r = 1.5 y ' ( 1 ) = y ( b ) = 0 b 2.86 b 2.99
r = 5.0 y ( 1 ) = y ' ( b ) = 0 b 2.55 b 2.54
r = 5.0 y ' ( 1 ) = y ( b ) = 0 b 2.49 b 2.54

Example 4.2.

Let us consider the following half-linear differential equation (4.5)Φ(y)+xΦ(y)=0,x1,a=1. The application of Corollary 2.6 and Theorem 3.1 to this example yields Table 2.

As it can be seen in Table 2, for both values r=1.5 and r=5, the bounds obtained using Corollary 2.6 improve those of Theorem 3.1 in the case y(1)=y'(b)=0 and are worse than those provided by Theorem 3.1 in the case y'(1)=y(b)=0.

Comparison of upper bounds for b in Example 4.2.

Corollary 2.6 Theorem 3.1
r = 1.5 y ( 1 ) = y ' ( b ) = 0 b 2.14 b 2.20
r = 1.5 y ' ( 1 ) = y ( b ) = 0 b 2.53 b 2.20
r = 5.0 y ( 1 ) = y ' ( b ) = 0 b 2.28 b 2.37
r = 5.0 y ' ( 1 ) = y ( b ) = 0 b 2.42 b 2.37

Example 4.3.

Let us consider the following half-linear differential equation: (4.6)(1xΦ(y))+xΦ(y)=0,x1,a=1. The application of Theorem 2.3 (with m(x)=x1/2(r-1)) and Theorem 3.1 to this example yields Table 3.

As it can be seen in Table 3, for the value r=5, the bounds obtained using Theorem 2.3 improve that of Theorem 3.1 in both cases y(1)=y'(b)=0 and y'(1)=y(b)=0. In the case of the exponent r=1.5, one gets exactly the opposite, the bounds obtained with Theorem 2.3 being slightly worse than those provided by Theorem 3.1 in both cases y(1)=y'(b)=0 and y'(1)=y(b)=0. It is worth remarking that in this latter case the formulae obtained by means of Theorem 2.3 are the same for both cases y(1)=y'(b)=0 and y'(1)=y(b)=0.

Comparison of upper bounds for b in Example 4.3.

Corollary 2.6 Theorem 3.1
r = 1.5 y ( 1 ) = y ' ( b ) = 0 b 2.012 b 2.01
r = 1.5 y ' ( 1 ) = y ( b ) = 0 b 2.012 b 2.01
r = 5.0 y ( 1 ) = y ' ( b ) = 0 b 2.32 b 3.09
r = 5.0 y ' ( 1 ) = y ( b ) = 0 b 2.43 b 3.09

Example 4.4.

Let us consider the following half-linear differential equation: (4.7)(exΦ(y))+exΦ(y)=0,x0,a=0. The application of Corollary 2.6 and Theorem 3.1 to this example yields Table 4.

As it can be seen in Table 4, for both values r=1.5 and r=5, the bounds obtained using Corollary 2.6 improve those of Theorem 3.1 in the case y(1)=y'(b)=0 and are worse than those provided by Theorem 3.1 in the case y'(1)=y(b)=0 (in fact in that case for the exponent r=1.5, Corollary 2.6 does not give any valid upper bound at all).

Comparison of upper bounds for b in Example 4.4.

Corollary 2.6 Theorem 3.1
r = 1.5 y ( 1 ) = y ' ( b ) = 0 b 1.325 b 1.33
r = 1.5 y ' ( 1 ) = y ( b ) = 0 No valid bound b 1.33
r = 5.0 y ( 1 ) = y ' ( b ) = 0 b 1.39 b 1.83
r = 5.0 y ' ( 1 ) = y ( b ) = 0 b 2.19 b 1.83

5. Comparison of Lower Bounds

Although the focus of this paper so far has been the calculation of upper bounds for the distance between adjacent zeroes and critical points, it is interesting to compare the lower bounds that can be obtained by means of Corollary 2.6 (for the sake of simplicity we will use such a corollary instead of Theorem 2.3) with similar results obtained by other authors (see [2, 49]). To facilitate this, we will concentrate on the case p(x)=1 in (1.2).

In all the cases under comparison, the lower bound can be determined from a Lyapunov inequality of the form AB, where A is an integral whose integrand contains the function q(x) of (1.2) and whose limits are the zero and the critical point, and B is a real value which depends only on r. Bearing this in mind, we can summarize the cases under analysis in Table 5.

Comparison of different lower bounds.

Reference A B Constraints
Corollary 2.6 ( b - a ) 1 / s ( a b max { q ( s ) ,    a s x } d x ) 1 / r , if y'(a)=y(b)=0. ( π r / 2 ) r - 1 r q ( x ) > 0
( b - a ) 1 / s ( a b max { q ( s ) ,    x s b } d x ) 1 / r , if y(a)=y'(b)=0.

[2, 7, 8], [9, Lemma  4.1] ( b - a ) 1 / s ( a b q ( x ) d x ) 1 / r 1 N/A

[9, Proposition  3.3] ( b - a ) 1 / s ( a b q ( x ) d x ) 1 / r π r / 2 q ( x ) > 0 , increasing

 a b | Q ( x ) | s ( b - x ) r - 1 d x with Q(x)=axq(s)ds, if y'(a)=y(b)=0. 1 N/A
a b | Q ( x ) | s ( x - a ) r - 1 d x with Q(x)=xbq(s)ds, if y(a)=y'(b)=0.

A first glance at Table 5 allows to notice the difficulty to compare Corollary 2.6 with the rest of methods due to the presence of max functions in the integrand of A. Such difficulty disappears (in general) in the case y(a)=y'(b)=0 when q(x) is decreasing, and in the case y'(a)=y(b)=0 when q(x) is increasing, since in both cases the max function becomes q(x) and the integral A which results is the same as the integral A which appears in other methods. If we focus on those cases, one can easily check that the bound obtained from Corollary 2.6 always improves those of [2, 7, 8] and [9, Lemma 4.1] (in all these latter cases, the value B is 1, which is lower than (πr/2)r-1r; this allows Corollary 2.6 to yield greater—and therefore better—lower bounds). It also improves that of [9, Proposition 3.3] in the case r>2 (for r2, one has r-1r1). That cannot be taken as a general rule if the character of q(x) is not as stated before, and in that case, it is relatively easy to find examples where the inequalities given in [2, 79] improve those obtained with Corollary 2.6.

In turn, the comparison with Yang’s extension (see, ) of Brown and Hinton’s version of the Lyapunov inequality to the half-linear case is more difficult. One can argue that in most cases  will give better (i.e., greater) lower bounds than those of Corollary 2.6. To see that, let us pick the case y'(a)=y(b)=0 and q(x) positive and increasing. Integrating by parts and applying Theorem 2.1 one can prove (5.1)ab|Q(x)|s(b-x)r-1dx=ab|Q(x)|s-1q(x)(b-x)rrdxab|Q(x)|s-1q(x)dxb-a(b-a)r+1r(r+1)dx=|Q(b)|s(b-a)rsr(r+1). Yang’s formula for such a case is (5.2)ab|Q(x)|s(b-x)r-1dx1, which combined with the previous equation gives (5.3)|Q(b)|1/r(b-a)1/s[r(r+1)r-1](r-1)/r2. The inequality (5.3) shares the same form of the integral A of Corollary 2.6 and a value B that is greater for 1r1,92 (in consequence for these cases, Yang’s formula will give much better lower bounds than Corollary 2.6) and slightly lower (2,4 percent lower in the worst case, which is a very low value) than the value B associated with Corollary 2.6 (in these cases, one can only guarantee that Corollary 2.6 will give better lower bounds than (5.3), but not that those lower bounds will be better than the ones calculated with Yang’s formula). That does not mean, of course, that Corollary 2.6 cannot ever give better lower bounds than , and in fact the case q(x) constant is a good counter example. However, it is true that the task of finding functions q(x) where Corollary 2.6 improves  is not easy at all.

Remark 5.1.

We have decided not to include the Lyapunov inequality (3.2) of Theorem 3.1 in the Table 5 to simplify the comparison among lower bounds. The reason behind is that, despite it sharing the form AB with A being an integral depending on q(x), there is no case where the integrand of A can be easily compared with the integrand of A of the rest of the cases, unlike what happens with Corollary 2.6. In consequence, the comparison has to be made numerically with concrete examples, which, although certainly interesting, would have increased the length of the paper excessively.

6. Conclusions

The methods described in this paper provide upper and lower bounds for the distance between a zero and a critical point of a solution of the half-linear differential equation (1.2).

The advantage of them over other methods ([2, 79]) which provide lower bounds for such a distance is quite clear if the monotonic conditions stated in Section 5 for q(x) are satisfied (it can be far from being really an advantage otherwise). In contrast, their advantage over the extension of Brown and Hinton’s method described by  in the same case is lesser, although it does exist for some (in fact for infinitely many) functions q(x).

In turn, their value in the determination of upper bounds is high, since no other alternatives to estimate such bounds seem to exist in the literature associated with (1.2), as far as the authors are aware.

As for the comparison between both methods of calculating upper bounds, the main disadvantage of the method associated with Theorem 2.3 is that it imposes the need for both functions p(x) and q(x) of (1.2) to be strictly positive on the interval where the calculation is performed. That does not happen in the case of the method associated with Theorem 3.1, which can be used even with p(x) or q(x) being negative. However, the application of this latter method is usually much more difficult due to the need to minimize b as a function of the parameter A, minimization which can be tricky and difficult to calculate, specially when compared with the method of Theorem 2.3.

The provided examples do not favour either method over the other: basically, it depends on the concrete case, and in fact in almost all of them (Example 4.3 is an exception), the method that was better in the case y(a)=y'(b)=0 proved to be worse in the case y'(a)=y(b)=0. We cannot guarantee that this is a general rule, anyway, although one can conjecture, just by simple examination of the resulting formulae, that cases where p(x) and q(x) are monotonic in the same sense may be likely to follow that pattern. More work is required in this area to find underlying rules, therefore.

Acknowledgment

This work has been supported by the Spanish Ministry of Science and Innovation Project DPI2010-C02-01.

Almenar P. Jódar L. An upper bound for the distance between a zero and a critical point of a solution of a second order linear differential equation Computers & Mathematics with Applications 2012 63 1 310 317 10.1016/j.camwa.2011.11.023 2863502 Došlý O. Řehák P. Half-Linear Differential Equations 2005 202 Amsterdam, The Netherland North-Holland xiv+517 Mathematics Studies 2158903 Elbert Á. Kusano T. Tanigawa T. An oscillatory half-linear differential equation Archivum Mathematicum 1997 33 4 355 361 1601353 Elbert Á. A half-linear second order differential equation Colloquia Mathematica Societatis János Bolyai 1979 30 158 180 Li H. J. Yeh C. C. Sturmian comparison theorem for half-linear second-order differential equations Proceedings of the Royal Society of Edinburgh A 1995 125 6 1193 1204 10.1017/S0308210500030468 1362999 Yang X. On inequalities of Lyapunov type Applied Mathematics and Computation 2003 134 2-3 293 300 10.1016/S0096-3003(01)00283-1 1931539 Lee C.-F. Yeh C.-C. Hong C.-H. Agarwal R. P. Lyapunov and Wirtinger inequalities Applied Mathematics Letters of Rapid Publication 2004 17 7 847 853 10.1016/j.aml.2004.06.016 2072845 Pinasco J. P. Lower bounds for eigenvalues of the one-dimensional p-Laplacian Abstract and Applied Analysis 2004 2 147 153 10.1155/S108533750431002X 2058270 Pinasco J. P. Comparison of eigenvalues for the p-Laplacian with integral inequalities Applied Mathematics and Computation 2006 182 2 1399 1404 10.1016/j.amc.2006.05.027 2282582 Almenar P. Jódar L. Improving explicit bounds for the solutions of second order linear differential equations Computers & Mathematics with Applications 2009 57 10 1708 1721 10.1016/j.camwa.2009.03.076 2528415 Moore R. A. The behavior of solutions of a linear differential equation of second order Pacific Journal of Mathematics 1955 5 125 145 0068690