Skip to content
Open Access Published by De Gruyter Open Access December 31, 2018

# An integral that counts the zeros of a function

Norbert Hungerbühler and Micha Wasem
From the journal Open Mathematics

# Abstract

Given a real function f on an interval [a, b] satisfying mild regularity conditions, we determine the number of zeros of f by evaluating a certain integral. The integrand depends on f, f′ and f″. In particular, by approximating the integral with the trapezoidal rule on a fine enough grid, we can compute the number of zeros of f by evaluating finitely many values of f, f′ and f″. A variant of the integral even allows to determine the number of the zeros broken down by their multiplicity.

MSC 2010: 30C15

## 1 Introduction

Counting the zeros of a given function f in a certain region belongs to the basic tasks in analysis. Iff: ℂ → ℂ is holomorphic, the Argument Principle and Rouché’s Theorem are tools which allow to find the number of zeros of f, counted with multiplicity, in a bounded domain of ℂ with sufficiently regular boundary (see, e.g. [4] for an overview of methods used for analytic functions). Descartes’ Sign Rule is a method of determining the maximum number of positive and negative real roots (counted with multiplicity) of a polynomial. The Fourier-Budan Theorem yields the maximum number of roots (counted with multiplicity) of a polynomial in an interval. Sturm’s Theorem, a refinement of Descartes’ Sign Rule and the Fourier-Budan Theorem, allows to count the exact number of distinct roots of a polynomial on a real interval (see, e.g., [5], [2], [8]). The mentioned methods are restricted to holomorphic functions and polynomials, respectively. On the other end of the regularity spectrum, for a merely continuous function f, the Theorem of Bolzano yields the information that at least one zero exists on an interval [a, b] if f has opposite signs at its endpoints, though, it does not count the zeros. Here, we want to construct a method which gives the number of zeros of a real function under only mild regularity assumptions. More precisely, we want to express the number of zeros of a function f by a certain integral (and boundary terms). The integrand depends on f, f′ and f″. If f is sufficiently regular, the integral (and hence the number of zeros of f) can be expressed by evaluating the integrand on a sufficiently fine partition of [a, b]. Modifications of the integral even allow to determine the number of the zeros broken down by their multiplicity.

To explain the basic idea, we consider the following elementary connection between the number of zeros of a periodic function and the winding number of the related kinematic curve in the state space with respect to the origin:

## Lemma 1.1

Let f : ℝ → ℝ be a 2π-periodic C2function with only simple zeros, i.e. points x withf(x) = 0 ≠ f′(x). Then, the number n of zeros of f in [0,2π)equals twice the winding number of the curve γ: [0, 2π) → ℝ2, x ↦ (f′(x), f(x)) with respect to the origin. Hence

n=1π02πf(x)2f(x)f(x)f(x)2+f(x)2dx.

Figure 1 illustrates a heuristic proof without words: Each colored arc between two zeros of 12 to the winding number of γ. In the sequel, we will rigorously prove much more general versions and variants of this result. We will develop integrals that count the number of zeros with and without multiplicity, and we will even be able to determine the number of zeros of a given multiplicity. As a byproduct, a coherent definition of a fractional multiplicity of zeros will be possible. To start with, it is necessary to analyze the nature of zeros of afunction.

### Figure 1

Number of zeros of f vs. winding number of (f′,f).

## 2 Zeros of Functions

A function f:(a, b) → ℝ may, in general, show a quite pathological behavior in the neighborhood of one of its zeros (see, e.g., Examples 2.2.3 and 2.9 below). To exclude such exotic cases but still be sufficiently general to cover most of the relevant cases, we use the following definition.

## Definition 2.1

A zero x0 ∈ (a, b) of a function fC0(a, b) ∩ C1((a, b)∖{x0}) will be called admissible provided

limxx0f(x)f(x)= and limxx0f(x)f(x)=.(2.1)

If f extends continuously to a (or b) and f(a) = 0 (or f(b) = 0), we will say that f has an admissible zero in a (or b) if

limxaf(x)f(x)=or limxbf(x)f(x)=.

## Remarks

1. An admissible zero is necessarily an isolated zero. In fact, if the zero x0 is an accumulation point of zeros of f then, by Rolle’s Theorem, it is also an accumulation point of zeros of f′ and the limits in Definition 2.1 cannot be plus or minus infinity.

2. The condition on the limits given in (2.1) is in fact equivalent to

limxx0ddxln|f(x)|=.(2.2)

Indeed, if 2.2) holds true, it follows that x0 is an isolated zero of f, hence f does not change its sign on (x0,x0+ε) and on (x0-ε,x0) for ε > 0 small enough. Moreover 0 < |f(x)| < |f′(x)| on a punctured neighborhood of x0. Hence, f′ cannot change sign and the claim follows by distinction of cases. The condition (2.2) is slightly more compact than (2.1), however, (2.1) is easier to handle in the calculations below.

3. A simple zero x0 ∈ (a, b) of fC1(a, b), i.e.f(x0) = 0 and f′(x0) ≠ 0 is admissible. It suffices to consider x0 = 0:

limx0f(x)f(x)=limx0f(0)+o(1)f(0)+xf(0)+o(x)=limx01xf(0)+o(1)f(0)+o(1)=.

The limit x ↗ 0 is analogous.

4. If f(x0) = f′(x0) = 0 and f′ is monotone on (x0,x0+ϵ) and on (x0−ϵ,x0) for some ϵ > 0, then x0 is an admissible zero: Indeed, for x0 < x < x0+ϵ and f′ non-decreasing (if f′ is non-increasing consider -f) on (x0, x0+ϵ), we have f(x) = x0xf′(t) dt ≤ (xx0)f′(x) and thus f(x)f(x)1xx0 for xx0. The argument for the limit xx0 is analogous.

5. If fCk(a, b) and x0 ∈ (a, b) is a zero of multiplicity k > 1, i.e.f()(x0) = 0 for all = 0,…, k−1 and f(k)(x0) ≠ 0, then x0 is admissible. This follows easily by an iterated application of L–Hôpital’s rule. Hence the zeros of real-analytic functions and a fortiori zeros of polynomials are admissible.

6. If f(x) = |xx0|αg(x) for a C1-function g with g(x0) ≠ 0 and 0 < α ∈ ℝ, then x0 is an admissible zero of f.

7. Every fC1([a, b]) can be extended to C1(I), where I ⊃ [a, b] is an open interval and the limits

limxaf(x)f(x) and limxbf(x)f(x)(2.3)

can be defined via , provided f(a),f(b) ≠ 0. If f has an admissible zero in a (or b), f can be extended antisymmetrically with respect to a (or b) to an extension for which a (or b) is an admissible zero. We will henceforth use this particular extension when computing limits like in (2.3).

## Example 2.2

1. The function f1C0(ℝ) ∩ C(ℝ∖{0}), x|x| has an admissible zero in x = 0 (see Remark 6 above).

2. The C-function

f2(x):=exp1x2,x00,x=0,

has an admissible zero of infinite multiplicity at x = 0 (see Remark 4 above).

3. An example of an isolated zero which is not admissible is given by the C-function

f3(x):=f2(x)(sin(1x3)+2),

which vanishes (together with all derivatives) in 0 but the corresponding limits (2.1) do not exist.

## Definition 2.3

A function f : [a, b] → ℝ belongs to 𝓐k([a, b]), k ∈ ℕ, if the following holds:

1. fC0([a, b]).

2. f has only admissible (and therefore finitely many) zeros x1 < … < xn and f|(xi,xi+1) (i = 1,…,n−1), f|(a,x1) and f|(xn,b) are of class Ck+1.

3. There exists a partition a = y1 < y2 < … < ym = b such that f|(yi,yi+1) is of class Ck+2 for all i = 1,…,m−1.

If f ∈ 𝓐0([a, b]), f will be called admissible.

## Remarks

1. Observe that 𝓐k+1([a, b]) ⊂ 𝓐k([a, b]) for all k ∈ ℕ by construction.

2. Every analytic function is in 𝓐([a, b]).

3. f:[−1,1] → ℝ, x|x| is in 𝓐([a, b]).

4. If f is admissible, then x ↦ (f′(x),f(x)) is not necessarily a continuous curve.

As a building block of the intended results we need the following: For σ ∈ [−∞,∞], let

H(x)=σxh(t)dt,(2.4)

where h:ℝ → ℝ is any piecewise continuous function such that the improper integral h(x)dx=1. Then we have the following theorem (recall (2.3) in order to make sense of the limits that appear).

## Theorem 2.4

Let f ∈ 𝓐0([a, b]). The number of zeros n(f) of f in [a, b] is given by

n(f)=abhf(x)f(x)f(x)2f(x)f(x)f(x)2dx+limxbHf(x)f(x)limxaHf(x)f(x)

and the number of zeros n̊(f) of f in (a, b) by

n˚(f)=abhf(x)f(x)f(x)2f(x)f(x)f(x)2dx+limxbHf(x)f(x)limxaHf(x)f(x).

## Proof

Consider first the case, where f(a),f(b) ≠ 0. Then the zeros of f are given by a < x1 < x2 < … < xn(f) < b. The integrand of

abhf(x)f(x)f(x)2f(x)f(x)f(x)2dx=:abI(x)dx

is a priori undefined whenever f vanishes or whenever f″ is undefined. We decompose the integral and compute the resulting improper integrals using unilateral limits. Since f is admissible, we have

xjxj+1I(x)dx=limxxjHxlimxxj+1Hx=1

for all j = 1,…, n(f)−1, where Hx: = H(f′(x)/f(x)). Integrating over a neighborhood of a point y where f″ is undefined does not introduce further boundary terms since limxy Hx − limxyHx = 0. Hence

abI(x)dx=ax1I(x)dx+j=1n(f)1xjxj+1I(x)dx+xn(f)bI(x)dx==Halimxx1Hx+(n(f)1)+limxxn(f)HxHb(2.5)

and therefore

n(f)=abI(x)dx+HbHa.(2.6)

The computation above suggests that n(f) > 1 but one can check that formula (2.6) holds true for n(f) = 1 and n(f) = 0 as well.

If f has zeros in a and b and therefore x1 = a, xn(f) = b, computation (2.5) gives

n(f)=abI(x)dx+1.(2.7)

According to (2.3), limxbHxlimxaHx=1 Hx = 1 and (2.7) becomes

n(f)=abI(x)dx+limxbHxlimxaHx

and hence (2.8) counts the zeros of f in [a, b] since it reduces to (2.6) if f(a),f(b) ≠ 0 and one can check that the remaining cases f(a) = 0≠ f(b) and f(a) ≠ 0 = f(b) are also covered. Let now

n˚(f)=abI(x)dx+limxbHxlimxaHx.(2.8)

Since

n!(f)n˚(f)=limxbHxlimxaHxlimxbHxlimxaHx==0,if f(a),f(b)01,if either f(a)=0 or f(b)=02,if f(a)=f(b)=0

we conclude that (f) counts the zeros of f in (a, b).□

## Remarks

1. Putting g(x) : = f′(x)/f(x), the integrand in Theorem 2.4 reads −(hg)(x) g′(x).With respect to the signed Borel-Lebesgue-Stieltjes-Measure dg(x): = g′(x) dx (see [9]), the integral can be written more compactly as

abh(g)dg.
2. If h(x): = 1/(π(1+x2)), i.e. h equals the Cauchy Density and f is an admissible 2π-periodic function, then the number n of zeros of f in [0,2π) equals

n=1π[02πf(x)2f(x)f(x)f(x)2+f(x)2dx+limx2πarctanf(x)f(x)limx0arctanf(x)f(x)]==1π02πf(x)2f(x)f(x)f(x)2+f(x)2dx,(2.9)

since the integral-free terms cancel out in this case. In this way we obtain Lemma 1.1 as a corollary of Theorem 2.4. Observe that a 2π-periodic C2 function with an odd number of zeros on [0,2π) gives rise to a curve x ↦ (f′(x),f(x)) having a half-integer valued winding number. This idea, further developed, leads to a generalized version of the Residue Theorem (see [3]).

Observe, that for a C2 function f with only zeros of multiplicity one, the integrand in (2.9) is continuous provided h is continuous. This remains true for zeros of higher multiplicity in the following way:

## Proposition 2.5

Let h:ℝ → ℝ be continuous and h(x) Cx2for|x|.Then, the integrand in Theorem 2.4

I:=hfff2fff2

is continuous if fCn([a, b]), n ≥ 2, has only zeros of multiplicityn.

## Proof

It suffices to show that I is continuous in 0 if x = 0 is a zero of f of multiplicity n. Then, by Taylor expansion, we have

f(x)=f(n)(0)n!+r0(x)xnf(x)=f(n)(0)(n1)!+r1(x)xn1f(x)=f(n)(0)(n2)!+r2(x)xn2

where ri are continuous functions with limx → 0ri(x) = 0. Using these expressions in I, we get

I(x)=hs1(x)xs2(x)x2

for continuous functions si with limx → 0si(x) = n. Thus

I(x)Cx2s12(x)s2(x)x2Cn

for x → 0.

If we only assume that h(x) = O(1/x2) for |x| → ∞ in the previous proposition, the proof shows that then I is at least bounded.□

As a corollary of Proposition 2.5 we obtain that if h is continuous and h(x) h(x)Cx2, then I is in C0 provided f is analytic. Nontheless, the function f may behave in the neighborhood of a zero in such a pathological way, that I becomes unbounded (see Example 2.7.3). This is why, in general, the integrals in Theorem 2.4 have to be interpreted as improper integrals. This means that the concrete computation requires the zeros of f to be known a priori in order to evaluate the improper integrals. It is therefore of practical importance to formulate conditions (see Propositions 2.8 and 2.10) with additional assumptions which guarantee that I is in L1: To this end we will slightly sharpen the admissibility condition for a function and impose some conditions on the behaviour of the zeros of f″ in neighborhoods of the zeros of f. Furthermore we will require h to have at least quadratic decay at infinity.

The proof of Proposition 2.5 for the case C = 1 indicates, how we can generalize the notion of multiplicity of zeros in a natural manner:

## Definition 2.6

The multiplicity μf(x0) of a zero x0 of f ∈ 𝓐0 is defined to be

μf(x0)=limxx0f(x)2f(x)2f(x)f(x).

Since the zeros of functions in 𝓐0 are admissible, it follows that μf(x0) ≥ qslant 0 whenever it exists, however, it can take values in [0,∞] (see Example 2.7.3 and 2.7.4 below). This definition of the multiplicity of a zero will be useful for a variant of Theorem 2.4 that takes the multiplicities of the zeros into account.

## Example 2.7

1. A function fCn, n ≥ 2 with 0 = f(x0) = f′(x0) = … = f(n−1)(x0) ≠ f(n)(x0) has a zero of multiplicity n in x0: the Definition 2.6 is compatible with the usual notion of multiplicity.

2. The function f(x) = |x|r, r > 0 has a zero of multiplicity r in x = 0.

3. The function

f(x)=1ln|x|,x00,x=0

has a zero of multiplicity 0 in x = 0.

4. The function f2 in Example 2.2.2 has a zero in x = 0 with μf2(0) = ∞.

## Proposition 2.8

Let h:ℝ → ℝ be a piecewise continuous function such that h(x) = O(1/x2) for |x| → ∞ and let f ∈ 𝓐0([a, b]) ∩ W2,1(a, b) have only zeros of positive multiplicity in the sense of Definition 2.6. Furthermore we assume that for each zero x0we have a neighborhood U such that either f″(x)≡ 0 on U∖{x0} or

k=1|zkx0|<,

where z1,z2,… denote the countably many zeros of f″ in U. Then

I:=hfff2fff2L1(a,b).

## Proof

Choose neighborhoods U1,…, Un of the n zeros of f, which do not (with the possible exception of the respective zero itself) contain singular points of f″ or zeros of f′ and let

U=i=1nUi.

Since |f| ≥ η for some η > 0 on the complement Uc and W2,1(a, b)↪ C1([a, b]) we can estimate

Uc|I(x)|dx=η2hL(R)(f2C0([a,b])|ba|+fC0([a,b])fL1(a,b))<.

Consider now wlog the neighborhood Ui of the zero xi = 0 and assume Ui = (−ε,ε) for some ε > 0. We need to show that I|(−ε,ε)L1. Since h(x) = O(1/x2) for |x| → ∞, there exists a constant C > 0 such that

|I(x)|C1+f(x)f(x)f(x)2.(2.10)

Note that ff″/f2L1(−ε,ε) if and only if N ∈ BV(−ε,ε), where N(x) = xf(x)/f′(x) denotes the Newton-Operator of f and BV(−ε,ε) denotes the space of functions g:(−ε,ε) → ℝ of bounded variation. It follows from the admissibility of the zero that N:(−ε,ε)∖{0} → ℝ can be continuously extended to N(0) = 0 and it holds that

N(x)=f(x)f(x)f(x)2,

for x ≠ 0. Let μ > 0 denote the multiplicity of the zero according to Definition 2.6. It holds that

limx0N(x)=μ1μ,μ<1,μ=.

According to the mean value theorem we have N(x)/x = N′(ξ) for some ξ between 0 and x and deduce that N ∈ C1(−ε,ε). The Taylor expansion of N around x = 0 is given by

N(x)=μ1μx+o(x),μ<x+o(x),μ=.

In any case there exists a constant K > 0 such that

|N(x)|K|x|,|x|<ε.(2.11)

We will now show that N ∈ BV([0,ε)), the argument on (−ε,0] being similar. We start by noticing that N is absolutely continuous on [δ,ε) for every 0 < δ < ε since x, f(x) and f′(x) are absolutely continuous and f′(x) ≠ 0 on [δ,ε). In particular, N ∈ BV([δ,ε)) for every 0 < δ < ε.

We will now distinguish two cases: If f″≡ 0 on (0,ε), then N ≡ 0 and we are done. In the remaining case we first consider the case when the set of zeros of f″ in (0,ε) is empty: Then N is monotone on [0,ε) and hence N ∈ BV([0,ε)). Otherwise the zeros of f″ in [0,ε) are given by z1 > z2 > … and we may set δ: = z1. According to (2.11) and since the zeros of f″ are precisely the zeros of N′ we can estimate the total variation of N on (zk+1,zk) by

zk+1zk|N(x)|dx2Kzk.

The total variation of N on [0,ε) is bounded by

k=1zk+1zk|N(x)|dx+δε|N(x)|dx2Kk=1zk+δε|N(x)|dx,

where the series converges by assumption and the integral is finite since N ∈ BV([δ,ε)). We conclude that N ∈ BV([0,ε)), which finishes the proof.□

## Remark

The key estimate (2.11) in the proof above follows from the admissibility and the positive multiplicity of the zeros. We will however formulate a variant of Proposition 2.8 below (Proposition 2.10), which covers admissible functions that have zeros of ill-defined multiplicity for which (2.11) still holds true: Take e.g. the C1 function f:xx3 ≤ ft(sin(1/x) + 2) + x which has an admissible zero in x = 0, but for which μf(0) does not exist, however, (2.11) holds true since f(x)/(xf′(x)) is bounded near 0 – in fact

limx0f(x)xf(x)=1.

Example 2.7.3. shows an admissible function for which (2.11) does not hold true. In the mentioned example, the first derivative is unbounded. But even functions with higher regularity may behave in such a pathological way near an admissible zero, that (2.11) does not hold true, as the following example shows:

## Example 2.9

Let

k(x)=x3+|x|7x3cosπlog2|x|,ifx00,ifx0.

Then f(x) = 0xk(t) dt is of class C3 and has an admissible zero in x = 0 but f(x)/(xf′(x)) is unbounded near 0.

## Proposition 2.10

Let h:ℝ → ℝ be a piecewise continuous function such that h(x) = O(1/x2) for |x| → ∞ and let f ∈ 𝓐0([a, b])∩ W2,1(a, b) be such that that for every zero x0of f there exists a relatively open neighborhood U ⊂ [a, b] such that

0<f(x)(xx0)f(x)<K~(2.12)

on U ∖{x0} and such that either f″≡ 0 on U∖{x0}, or

k=1|zkx0|<,

where z1,z2,… denote the countably many zeros of fin U∖{x0}. Then

I:=hfff2fff2L1(a,b).

## Proof

Choose neighborhoods U1,…, Un of the n zeros of f, which do not (with the possible exception of the respective zero itself) contain singular points of f″ or zeros of f′ such that (2.12) holds on each punctured neighborhood. As in the proof of Proposition 2.8 we obtain ∥I∥L1(Uc) < ∞, where U = U1∪ … ∪ Un and the estimate (2.10). Let wlog 0 be a zero of f and let (−ε,ε) be its respective neighborhood for some ε > 0. As in the proof of Proposition 2.8, we are done if we show that N ∈ BV([0,ε)). The condition 0 < |f(x)/(xf′(x))| < on (−ε,ε)∖{0} implies that

0<f(x)f(x)<K~|x|,(2.13)

from which we conclude that N extends continuously to [0,ε) (where N(0) = 0) and

|N(x)|(K~+1)x,x[0,ε).(2.14)

This is just estimate (2.11) with K = + 1. The rest of the proof is exactly the same as the one of Proposition 2.8.□

## 3 Counting Zeros with Multiplicities

Let again h:ℝ → ℝ be a piecewise continuous function such that h(x), dx = 1 and define H as before in (2.4). Moreover, let

Ig(x)=hf(x)f(x)g(x)f(x)2f(x)f(x)f(x)2Hf(x)f(x)g(x),g1(x)=f(x)2f(x)2f(x)f(x)+cf(x)2,g2(x)=expf(x)2f(x)f(x)f(x)2+f(x)2,

where c ∈ ℝ. Note that if x0 is a zero of multiplicity μf(x0), then g1(x) → μf(x0) as xx0 for every value c in the definition of g1 and if μf(x0) > 0, then g2(x) → exp (1μf(x0)) as xx0.

## Lemma 3.1

Let all the zeros of f ∈ 𝓐0([a, b])∩ C2([a, b]) have well-defined multiplicities. Then there exists c ∈ ℝ such that g1has no poles.

## Proof

If x0 is a zero of f, we have that g1(x) → μf(x0) as xx0. In other words g1 extends continuously to the zeros of f. Hence there are open neighborhoods of the zeros of f, where g1 has no poles. On the complement of these neighborhoods, there exists a number δ > 0 such that |f(x)| ≥ δ. Hence f′(x)2+cf(x)2f′(x)2+2. If we choose c > δ−2ff″∥C0([a, b]), then g1 has no poles. In particular, if f is analytic, this choice of c ensures that g1 is analytic as well.□

We have the following theorem for analytic functions f:[a, b] → ℝ:

## Theorem 3.2

Let f:[a, b] → ℝ be an analytic function and choose c in the definition of g1such that g1is analytic. If h(x) = O(1/x2) for |x| → ∞, then Ig1, Ig2L(a, b) and if f has n zeros of multiplicity ℓ in [a, b] and n̊ zeros of mutliplicity ℓ in (a, b), then

abIg1(x)dx+limxbHf(x)f(x)g1(x)limxaHf(x)f(x)g1(x)==1n,abIg1(x)dx+limxbHf(x)f(x)g1(x)limxaHf(x)f(x)g1(x)==1n˚,abIg2(x)dx+limxbHf(x)f(x)g2(x)limxaHf(x)f(x)g2(x)==1nexp1,abIg2(x)dx+limxbHf(x)f(x)g2(x)limxaHf(x)f(x)g2(x)==1n˚exp1.

## Proof

We first prove the L-bounds: It suffices to show that Ig1 and Ig2 are bounded near the zeros of f. Let x0 be a zero of multiplicity k and write (locally) f(x) = (xx0)kj(x), where j is analytic and j(x0)≠0. Since

limxx0g1(x)=2j(x0)j(x0)

we find the limits

limxx0Hf(x)f(x)g1(x)=2j(x0)j(x0)limxx0Hf(x)f(x)g1(x)=0.

If

hf(x)f(x)g(x)f(x)2f(x)f(x)f(x)2

is bounded near x0, the claim follows. Since ≤ |h≤ (f′(x)}/{f(x))|≤ Cf(x)2}/{f′(x)2 and

limxx0C|g1(x)|f(x)2+|f(x)f(x)|f(x)2=C(2k1),

we obtain Ig1L(a, b). For Ig2, observe that

limxx0g2(x)=2exp1kj(x0)k2j(x0)

and therefore

limxx0Hf(x)f(x)g2(x)=2exp1kj(x0)k2j(x0)limxx0Hf(x)f(x)g2(x)=0.

Proceeding as for g1 we find

limxx0C|g2(x)|f(x)2+|f(x)f(x)|f(x)2=Cexp1k2k1k

and hence Ig2L(a, b). The computation of the integrals is done as in the proof of Theorem 2.4.□

## Remark

If f ∈ 𝓐1([a, b])∩ C2([a, b]) only has zeros of well-defined multiplicities and if the set of zeros of f in (a, b) is given by and the set of zeros of f in [a, b] by N, then

abIg1(x)dx+limxbHf(x)f(x)g1(x)limxaHf(x)f(x)g1(x)=xNμf(x),abIg1(x)dx+limxbHf(x)f(x)g1(x)limxaHf(x)f(x)g1(x)=xN˚μf(x).

## Lemma 3.3

Let 𝓝 be the set of sequences with natural entries of which only finitely many are non-zero. Then the map 𝓕:𝓝 → ℝ defined by 𝓕(k1,…) = =1kexp1is injective.

## Proof

The difference 𝓕(k1,… )−𝓕( k1,… ) is equal to the finite sum

=1(kk)exp1.

If this sum vanishes, k = k for all by the von Lindemann-Weierstrass theorem (see [7, §3]).□

## Corollary 3.4

Let f:[a, b] → ℝ be analytic. If f has nℓ zeros of multiplicity ℓ in [a, b] and n̊ zeros of mutliplicity ℓ in (a, b), then

(n1,)=F1abIg2(x)dx+limxbHf(x)f(x)g2(x)limxaHf(x)f(x)g2(x)(n˚1,)=F1abIg2(x)dx+limxbHf(x)f(x)g2(x)limxaHf(x)f(x)g2(x).

## Example 3.5

Let f(x) = cos(2 x) + x2 Sin (2x) 12ex+x24. Using Theorem 2.4 and 3.2 on [0,2π] we obtain

=1n=3,=1n˚=2,=1n=4,=1n˚=2.

and we conclude that f has two zeros in (0,2π) and a double zero on the boundary of [0,2π].

## Example 3.6

Let f(x) = x7−2x6+x5x3+2x2x have n zeros of multiplicity on ℝ. By Theorem 2.4 and 3.2 on ℝ (observe that the boundary terms of the integrals cancel out in this case) we find that

=1n=3 and =1n=5.

Hence (n1,…) either equals (1,2,0,…) or (2,0,1,…). In particular n = 0, for ≥ 4. Using again Theorem 3.2 we get

i=13niexp1i6.8322.

Since 1⋅ e + 2 ⋅ e ≈ 6.0157 and 2⋅ e + 1⋅ e3 ≈ 6.8322 we conclude that f has two simple zeros and one of multiplicity 3.

## 4 Numerical Aspects

The number of zeros of a function f in a given interval [a, b] is of course an integer. Therefore is suffices to compute the integral in Theorem 2.4 with an error ε < 12. In particular, for the trapezoidal rule

TN(I):=baN(I(a)+I(b)2+k=1N1I(a+kbaN))

with N + 1 equidistant grid points, the error ε(N) is estimated by

ε(N)=|abI(x)dxTN(I)|(ba)312N2IL

(see, e.g., [6] or [1]). Thus we have

## Theorem 4.1

Let f satisfy the assumptions in Theorem 2.4. If

N>(ba)36IL,

then one can replace the integral in Theorem 2.4 by the finite sum TN(I) and round the result to the closest integer to get the values n(f) and n̊(f),respectively.

This theorem is quite remarkable: It allows to compute the number of zeros of a function f on [a, b] by evaluating finitely many values of f, f′ and f″.

## Example 4.2

Let f:ℝ → ℝ, xJ0(x), be the zeroth Bessel function of the first kind. If h is the Cauchy density, one can verify that IL<1π. We want to compute the number of zeros of J0 on [0,2π] by Theorem 4.1. It suffices to employ the trapezoidal rule with only

N=2π3=4

equidistant intervals. We find

T4(I)=π2I(0)+I(2π)2+k=13I(kπ2)1.76479

and thus

T4(I)1πarctanJ1(2π)J0(2π)1.76479+0.24419=2.00898

and hence, J0 has two zeros on [0,2π].

If we compute the number of zeros of J0 on [0,100π], we have to choose

N=50063π=1283.

(Actually, a finer analysis shows that a much smaller N suffices). In this case, we get

T1283(I)=100π1283I(0)+I(100π)2+k=11282I(k100π1283)99.75013

and

T1283(I)1πarctanJ1(100π)J0(100π)99.75013+0.24987=100,

hence we conclude that J0 has n = 100 zeros on [0,100π], in accordance with the well known distribution of zeros of J0. Surprisingly, the routine CountRoots of Mathematica is giving up on this simple problem after giving it some thought.

From a practical point of view, it is desirable to keep ∥I″∥L (and hence N) as small as possible. This can be achieved in several ways: First of all, we have the freedom to choose the function h. Below there is a small table of possible choices of h and the resulting function H in Theorem 2.4 (in each case, the integrand I turns out rather nicely).

Moreover, with smooth functions γ and κ that satisfy signγ(x) = signκ(x) = sign x for all x ≠ 0 and γ(x)∼ C1 |x|αsgn x and κ(x)∼ C2 |x|β sgn x as x → 0, where 0 < αβ, one can modify the integrand I as follows and the proof of Theorem 2.4 still goes through:

I(x)=hγ(f(x))κ(f(x))γ(f(x))f(x)κ(f(x))γ(f(x))f(x)κ(f(x))κ(f(x))2.

In this case the boundary terms in a and b have to be taken with the function

Hγ(f(x))κ(f(x)).

h(x)H(x)
1π(1+x2)arctanxπ
12x2+13/2x2x2+1
expx2π12erf(x)
14x214x24x2+14x2+114x
sech(2x)212tanh(2x)
ex(1+ex)211+ex
UnitBox(x)0if2x<1x+12if12<x121if2x>1
UnitTriangle(x)0ifx112(1+x)2if1<x012(1x)2+1if0<x11if1<x

# Acknowledgement

We would like to thank the referees for their valuable remarks which greatly helped to improve this article.

### References

[1] Gautschi W., Numerical analysis. An introduction, 1997, Birkhäuser Boston, Inc., Boston.Search in Google Scholar

[2] Henrici P., Applied and computational complex analysis. Volume 1: Power series—integration—conformal mapping—location of zeros, Pure and Applied Mathematics, 1974, Wiley-Interscience [John Wiley & Sons], New York-London-Sydney.Search in Google Scholar

[3] Hungerbühler N., Wasem M., A generalized version of the residue theorem, ArXiv e-prints 1808.00997, August 2018.Search in Google Scholar

[4] Kravanja P., Van Barel M., Computing the zeros of analytic functions, volume 1727 of Lecture Notes in Mathematics, 2000, Springer-Verlag, Berlin.10.1007/BFb0103927Search in Google Scholar

[5] Obreschkoff N., Verteilung und Berechnung der Nullstellen reeller Polynome, 1963, VEB Deutscher Verlag der Wissenschaften, Berlin.Search in Google Scholar

[6] Schwarz H.R., Numerical analysis. John Wiley & Sons, Ltd., Chichester, 1989, A comprehensive introduction, with a contribution by Jörg Waldvogel.Search in Google Scholar

[7] Weierstraß, K., Zu Lindemann’s Abhandlung: „Über die Ludolph’sche Zahl“. Sitzungsberichte der Königlich Preußischen Akademie der Wissenschaften zu Berlin, 1885, 5, 1067–1085.Search in Google Scholar

[8] Sturmfels B., Solving systems of polynomial equations, volume 97 of CBMS Regional Conference Series in Mathematics. Published for the Conference Board of the Mathematical Sciences, 2002, Washington, DC; by the American Mathematical Society, Providence, RI.10.1090/cbms/097Search in Google Scholar

[9] Dshalalow J.H., Real Analysis: An Introduction to the Theory of Real Functions and Integration, Studies in Advanced Mathematics, 2000. ISBN=9781584880738, Taylor & Francis.10.1201/9781420036893Search in Google Scholar

Received: 2018-09-19
Accepted: 2018-11-29
Published Online: 2018-12-31

© 2018 Hungerbühler and Wasem, published by De Gruyter

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.