Mathematicians have traditionally dealt with convergent series and shunned divergent ones. But, long ago, astronomers found that divergent expansions yield valuable results. If these so-called asymptotic expansions are truncated, the error is bounded by the first term omitted. Thus, by stopping just before the smallest term, excellent approximations may be obtained.

**Astronomical Series**

Many of the series that astronomers have used do not converge. Henri Poincaré pointed out the difference in viewpoint between geometers and astronomers — or between pure and applied mathematicians. He considered first the series

This is the Taylor series for , and it is convergent for all values of . However, it is a nightmare to compute, as the terms grow rapidly before decreasing to zero.

In contrast, the series

diverges for all nonzero . However, the terms decrease rapidly at first, and only after exceeds do they begin to increase. The first twenty terms of the two series are shown for .

** Convergence & Divergence **

For the infinite series

we define the sequence of partial sums

The series is convergent if the sequence of partial sums has a finite limit. If the sequence has no such finite limit, the series is divergent. The sequence may oscillate regularly or irregularly or tend to an infinite limit.

Obviously, we cannot compute the sum of a divergent series in the usual way. However, there are several “summation methods” for associating a number with a divergent series. Leonhard Euler was a maestro of divergent series, and he obtained many extraordinary results using ingenious methods.

**An Example**

One of Euler’s examples was the divergent series

or, more explicitly

It is clear that the partial sums oscillate between positive and negative values with ever-increasing amplitude, so the series is divergent. Nevertheless, Euler assigned the value to the sum of the series. He arrived at this value using no fewer than five different summation methods. We will look at one of these, based on the solution of a differential equation. Let us consider the linear first-order ordinary differential equation

**Solution in Series**

There are several ways to solve this differential equation. First, let’s seek a solution in the form of a power series

Substituting this into the equation, and equating the coefficients of different powers of , we find

So the series expansion becomes

Setting , this is just Euler’s series that we started with above:

Applying the ratio test to the series we have

which has the limit for all nonzero , so the series diverges unless . In particular, the value is not defined by the series.

**An Integrating Factor**

Another way to solve the differential equation is by means of an integrating factor. We write the equation as

Now multiply by the exponential of the integral of the coefficient of :

The equation becomes

and the solution follows:

A non-obvious change of variable or converts this to

This can be expressed in terms of the standard exponential integral , but we will not worry about that. Expanding the component by the binomial theorem and integrating term by term, we get

the same series expansion that we obtained above. The series

is an asymptotic expansion. It is divergent but, if we stop just before the smallest term, we get an approximation to the true value.

**Numerical Results**

Suppose ; then the terms decrease until and then increase. Summing the first ten terms and comparing it to the integral solution, we get

accurate to three significant figures. Now, taking , the terms decrease until . Summing the first 99 terms and comparing it to the integral solution, we get

which is accurate to (at least) six significant digits.

Let’s recap: we have a solution of the ode in terms of a divergent series

and a solution in terms of an integral

The integral is well defined and, for , it has the value . This is the value that Euler associates with the divergent series

** Sources **

Étienne Ghys, 2017: *A Singular Mathematical Promenade.* ENS Éditions, 306pp. ISBN: 978-2-8478-8940-6.