It is easy to prove that the sum of the infinite geometric series $$\frac{D(1+g)^n}{(1+r)^n}$$ is $$S = \frac{D(1+g)}{r-g}$$ where n = 1, 2, 3, ... and 0 < g < r < 1.
Now, if the infinite series becomes:
$$\frac{D(1+g_0)^n}{(1+r)^n}, n = 1, 2, ... , A$$
$$\frac{D(1+g_0)^A(1+g_1)^{n-A}}{(1+r)^n}, n = A+1, A+2, ...$$
Show that the sum approximates the following (assume $0< g_1 < g_0 < r < 1$):
$$S = \frac{D}{r-g_1}[(1+g_1) + A(g_0-g_1)]$$


I have figured out graphically and suspect that the authors did so similarly.
We are basically trying to approximate the area under the solid blue and black curves (leaving out the discounting by $(1+r)^n$ first, hence the exponential growth).
This area is approximately the same as (a bit smaller than) the area under the blue curve (solid and dashed). To calculate that area, we need to figure out the value at which the blue curve intersects the vertical axis.
Since we already know the value at which the black curve intersects the vertical axis: $D(1+g_0)$, we just need to calculate the gap as marked on the graph. The gap is approximately the difference between $a_A$ and $a_A^*$ if $a_1$ had grown at $g_1$. Then
$$Gap = \frac{D(1+g_0)^A}{(1+g_1)^A}$$
Linear approximation gives:
$$Gap = D[(1 + Ag_0) - (1 + Ag_1)] = DA(g_0 - g_1)$$
Going back to the Gordon Growth Model, all the calculation we just did were for the purpose of adjusting the numerator, which corresponds to the vertical axis intersect on our graph.
Voila! The numerator is now approximated as $D(1+g_0) + DA(g_0 - g_1)$
Finally we have the approximation of $S$:
$$\frac{D[(1 + g_0) + A(g_0 - g_1)]}{r - g_1}$$