JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

Laurent Series, Singularities & Residues

When a function blows up at a point, a Taylor series can't follow it there — but a Laurent series can. We add negative powers, read off how the function misbehaves, and capture the whole story of a pole in a single number: the residue.

Where the Taylor series runs out

By now you know that an analytic function — one that is complex-differentiable on a region — is fabulously rigid: near any interior point it equals its Taylor series, a sum of non-negative powers (z - z0)^n that you first met for real functions in Volume I. But a Taylor series is an optimist. It only knows how to describe a function that is well-behaved at z0; it has no vocabulary for a point where the function blows up. Consider f(z) = 1/z near z0 = 0: there is simply no honest power series a0 + a1 z + a2 z^2 + ... that equals 1/z, because every such series stays finite at 0 while 1/z races off to infinity.

The fix is delightfully simple: allow negative powers too. A Laurent series about z0 is a two-way sum, ... + a_{-2}/(z - z0)^2 + a_{-1}/(z - z0) + a_0 + a_1 (z - z0) + a_2 (z - z0)^2 + ..., running over all integers n from minus infinity to plus infinity. The piece with negative n — called the principal part — is exactly the new vocabulary that can describe a blow-up. For 1/z the whole Laurent series is just the single term 1/z. The Laurent series is the right local picture of a function on a ring (an annulus) that surrounds a bad point but doesn't include it.

Three flavors of singularity

An isolated singularity is a point z0 where f fails to be analytic, but where f IS analytic on a punctured neighborhood around it. The shape of the principal part sorts these into exactly three kinds — this is the classification of singularities. If the principal part is empty (no negative powers at all), the singularity is removable: the function only looked broken, and a single value plugs the hole. The classic is sin(z)/z at z0 = 0, whose Laurent series 1 - z^2/6 + z^4/120 - ... has no negative powers, so defining the value to be 1 there makes it perfectly analytic.

If the principal part has finitely many terms — it stops at some 1/(z - z0)^m with a_{-m} not zero — the singularity is a pole of order m. Here f genuinely runs to infinity, but in a controlled, polynomial way: multiplying by (z - z0)^m tidies it back into a well-behaved analytic function. A pole of order 1 is called simple; 1/z has a simple pole, and 1/z^3 has a pole of order 3. Poles are the friendly singularities — almost everything useful in applied complex analysis happens at a pole.

And if the principal part never stops — infinitely many negative powers — the singularity is essential, and the function behaves wildly. The headline example is e^{1/z} at z0 = 0, whose Laurent series 1 + 1/z + 1/(2! z^2) + 1/(3! z^3) + ... has an endless principal part. Picard's theorem says that near an essential singularity f takes every complex value (with at most one exception) infinitely often, in any tiny ring you choose. These are not friendly, and we mostly steer around them.

The one coefficient that matters: the residue

Among all the Laurent coefficients, one is royalty: a_{-1}, the coefficient of 1/(z - z0). It is called the residue of f at z0, written Res(f, z0). Why this one and not the others? Because of a small miracle you saw building in the previous guides. If you integrate (z - z0)^n once around a small circle enclosing z0, the answer is zero for EVERY integer power n except n = -1, where it equals 2 pi i. So when you integrate a whole Laurent series term by term around the singularity, every term cancels to nothing except the 1/(z - z0) term — and what survives is exactly 2 pi i times a_{-1}.

So the residue is, almost by design, the part of the function that an contour integral can actually feel. Everything else around the singularity integrates to zero; the residue is the irreducible leftover, the 'flux' the pole pushes out through any loop around it. This is the bridge that turns a local algebraic fact — one Laurent coefficient — into a global statement about an integral, and it is the seed of the residue theorem that the next guide will harvest.

Simple pole at z0:      Res(f, z0) = lim_{z->z0} (z - z0) f(z)

Order-m pole at z0:     Res(f, z0) = 1/(m-1)! * lim_{z->z0} d^{m-1}/dz^{m-1} [ (z - z0)^m f(z) ]

Quotient g/h, simple pole where h(z0)=0, h'(z0)!=0:
                        Res(g/h, z0) = g(z0) / h'(z0)
Three workhorse formulas for residues. The first is just the order-m rule with m = 1; the third is the fastest route for a ratio whose denominator has a simple zero.

Reading a residue off a function

You rarely compute a residue by writing out the full Laurent series — that is the slow path. Instead you use shortcuts that extract a_{-1} directly. For a simple pole, multiply away the pole and take a limit: Res(f, z0) = lim as z goes to z0 of (z - z0) f(z). That single multiplication strips off the lone 1/(z - z0) and reveals its coefficient. The trick rests on the same limit machinery from the derivative in Volume I, just carried into the complex plane.

  1. Find the singular point and its order. Factor the denominator; if (z - z0) appears to the first power and the numerator is nonzero there, it is a simple pole.
  2. Worked simple pole: take f(z) = e^z / (z^2 + 1). The denominator factors as (z - i)(z + i), so z0 = i is a simple pole. Multiply: (z - i) f(z) = e^z / (z + i).
  3. Take the limit z -> i: Res(f, i) = e^i / (2i). That is the whole residue — a single complex number summarizing the pole at i.
  4. For a higher-order pole of order m, instead differentiate m-1 times: Res(f, z0) = 1/(m-1)! times the limit of the (m-1)-th derivative of (z - z0)^m f(z). The clean power (z - z0)^m cancels the pole; differentiating digs down to the a_{-1} layer.

When f is a ratio g(z)/h(z) and h has a simple zero at z0 (so h(z0) = 0 but h'(z0) is not zero), there is an even slicker formula: Res(g/h, z0) = g(z0) / h'(z0). In the example above, g = e^z and h = z^2 + 1 with h'(z) = 2z, giving e^i / (2i) instantly — no factoring needed. This shortcut is a direct cousin of the Cauchy integral formula you met earlier, which already taught you that a function's values on a loop control everything inside.

Why this matters, and honest fine print

The reason to care is breathtakingly practical: residues let you compute real integrals that defy every real-variable method. Many a definite integral from 0 to infinity — the kind that has no elementary antiderivative — falls instantly once you close the path into a contour and add up the residues inside. That is the whole point of the evaluation of real integrals technique ahead. Moving into the complex plane to crack a real problem is not a cheat or a trick; it is a fully rigorous method, and the residue is the gear that makes it turn.

A few honest cautions. First, the residue formulas above apply to poles only — an essential singularity has a residue too (it is still a_{-1}), but the limit-and-derivative shortcuts do not work; you must dig the coefficient out of the actual Laurent expansion. Second, get the order right: applying the simple-pole formula to a double pole silently gives the wrong number. Third, residues are about ISOLATED singularities. A branch point — like the one at 0 for the complex logarithm — is not isolated, has no Laurent series, and no residue in this sense; it needs branch cuts and different tools entirely.

Step back and see the shape of the idea. A Laurent series widens the lens of the power series just enough to look at a function right where it breaks; the principal part is a precise field guide to HOW it breaks; and the residue distills that breakage into one number that an integral can read. With this in hand you are ready for the residue theorem — sum the residues inside a loop, multiply by 2 pi i, and you have the integral — the slickest tool in the applied-mathematician's kit.