Zeros, singularities, and residues

A function that blows up is not automatically a mystery — complex analysis gives a complete, three-way classification of every point where an analytic function misbehaves, and attaches to each such point a single extractable number, the residue, that will later do the work of an entire contour integral. The classification comes for free once you have the Laurent series about the point: look only at its negative-power terms, the principal part, and every isolated singularity sorts itself into removable, pole, or essential. Zeros get a matching, dual classification by order. This article builds both classifications from the Laurent series you already met, derives the shortcut formulas that let you find a residue without ever writing out a full series, and works a Laurent-based example alongside the formula-based one to show the two methods agree. It sets up the Residue Theorem, the topic immediately following.

The core insight

Every isolated singularity of an analytic function is classified — removable, a pole of order mm, or essential — by how many negative-power terms sit in its Laurent series' principal part, and the residue, the coefficient c−1c_{-1}, is the one number from that series that survives contour integration.

The core idea

Zeros and their order. A point z0z_0 is a zero of order mm of an analytic function ff if f(z0)=f′(z0)=⋯=f(m−1)(z0)=0f(z_0) = f'(z_0) = \cdots = f^{(m-1)}(z_0) = 0 but f(m)(z0)≠0f^{(m)}(z_0) \ne 0. Equivalently, near z0z_0 the function factors as

f(z)=(z−z0)mg(z),f(z) = (z-z_0)^m g(z),

where: gg is analytic at z0z_0 and g(z0)≠0g(z_0) \ne 0 — the order mm is exactly the power of (z−z0)(z-z_0) you can pull out before what remains stops vanishing. A zero of order 11 is called a simple zero. This factored form is the mirror image of the pole formula below, and the two are used together constantly: a zero of q(z)q(z) in a denominator is precisely what creates a pole of f=p/qf = p/q.

Isolated singularities and the Laurent test. z0z_0 is an isolated singularity of ff if ff is analytic on a punctured disk 0<∣z−z0∣<r0 < |z-z_0| < r but not analytic (or not even defined) at z0z_0 itself. Expand ff there as a Laurent series,

f(z)=∑n=−∞∞cn(z−z0)n,f(z) = \sum_{n=-\infty}^{\infty} c_n (z-z_0)^n,

where: the terms with n<0n<0 form the principal part. Reading only that principal part sorts every isolated singularity into exactly one of three types. A removable singularity has principal part identically zero — every cn=0c_n=0 for n<0n<0 — so lim⁡z→z0f(z)\lim_{z\to z_0}f(z) exists and defining f(z0)f(z_0) to equal it makes ff analytic there; f(z)=sin⁡z/zf(z)=\sin z/z at z=0z=0 is the standard example. A pole of order mm has a principal part that stops after finitely many terms, with c−m≠0c_{-m}\ne 0 and every cn=0c_n=0 for n<−mn<-m; equivalently f(z)=h(z)/(z−z0)mf(z) = h(z)/(z-z_0)^m with hh analytic and h(z0)≠0h(z_0)\ne 0. A pole of order 11 is a simple pole. An essential singularity has infinitely many nonzero negative-power terms, as in f(z)=e1/zf(z)=e^{1/z} at z=0z=0, whose series ∑n≥01n! zn\sum_{n\ge 0} \frac{1}{n!\,z^n} never terminates going leftward. The figure below lays these three cases out as bars along the Laurent index nn, making the "how far left do the nonzero bars extend" test visual.

Figure
Bar chart of Laurent series coefficients against index n from minus three to three showing a removable singularity with no bars left of zero, a pole of order two with exactly two bars left of zero, and an essential singularity with bars continuing indefinitely to the left
Classifying a singularity from its Laurent principal part removable -3 0 3 no bars for n less than 0 — principal part empty pole of order 2 bars stop at n = -2, the pole order m = 2 essential -3 0 3 bars never stop going left — infinitely many terms Purple and red bars mark the principal part (negative n); blue bars are the analytic part (n greater than or equal to 0)

Residue. The residue of ff at an isolated singularity z0z_0, written Res⁡z=z0f(z)\operatorname{Res}_{z=z_0}f(z), is defined as the Laurent coefficient c−1c_{-1}. It is the only coefficient that matters for contour integration because ∮C(z−z0)n dz=0\oint_C (z-z_0)^n\,dz = 0 for every integer n≠−1n\ne -1 around a small loop about z0z_0, but equals 2πi2\pi i when n=−1n=-1 — so integrating the whole Laurent series term by term kills every term except the c−1c_{-1} one. A removable singularity always has residue 00, since its principal part vanishes entirely.

Worked example

Problem. Find the residue of f(z)=ez−1z3f(z) = \dfrac{e^z-1}{z^3} at z=0z=0, first by expanding a Laurent series directly and then by the order-mm pole formula, and confirm the two answers agree.

Step 1. Expand the numerator with the known Maclaurin series ez=1+z+z22!+z33!+z44!+⋯e^z = 1+z+\dfrac{z^2}{2!}+\dfrac{z^3}{3!}+\dfrac{z^4}{4!}+\cdots, so ez−1=z+z22+z36+z424+⋯e^z-1 = z+\dfrac{z^2}{2}+\dfrac{z^3}{6}+\dfrac{z^4}{24}+\cdots.

Step 2. Divide by z3z^3 term by term: f(z)=1z2+12z+16+z24+⋯f(z) = \dfrac{1}{z^2}+\dfrac{1}{2z}+\dfrac{1}{6}+\dfrac{z}{24}+\cdots. The principal part has exactly two nonzero terms, 1/z21/z^2 and 1/(2z)1/(2z), so z=0z=0 is a pole of order 22, and reading off the coefficient of 1/z1/z gives Res⁡z=0f(z)=12\operatorname{Res}_{z=0}f(z) = \dfrac{1}{2} directly from the series.

Step 3. Cross-check with the order-mm formula (derived below) for m=2m=2: Res⁡z=0f(z)=lim⁡z→0ddz ⁣[z2f(z)]=lim⁡z→0ddz ⁣[ez−1z]\operatorname{Res}_{z=0}f(z) = \lim_{z\to 0}\dfrac{d}{dz}\!\left[z^2 f(z)\right] = \lim_{z\to 0}\dfrac{d}{dz}\!\left[\dfrac{e^z-1}{z}\right]. Writing ez−1z=1+z2+z26+⋯\dfrac{e^z-1}{z} = 1+\dfrac{z}{2}+\dfrac{z^2}{6}+\cdots from Step 1, its derivative is 12+z3+⋯\dfrac{1}{2}+\dfrac{z}{3}+\cdots, which at z=0z=0 gives 12\dfrac{1}{2}. Both routes agree. The residue of f(z)=(ez−1)/z3f(z) = (e^z-1)/z^3 at z=0z=0 is 12\dfrac{1}{2}, and z=0z=0 is a pole of order 22.

Competitive stretch

The same function also has an effective order lower than the denominator's bare factor of 33 suggests, precisely because the numerator itself vanishes at z=0z=0 (a simple zero, since ez−1≈ze^z-1 \approx z there): the pole order is 3−1=23-1=2, not 33. Spotting this cancellation before applying any residue formula is exactly the habit tested in the traps below.

Common traps

  • Assuming every singularity is a pole. Check the Laurent expansion, or at least the vanishing/non-vanishing of a limit, before reaching for a pole-residue formula — essential and removable singularities need entirely different treatment, and the pole formulas give nonsense if misapplied to them.
  • Misjudging pole order from the denominator's factored form alone. If the numerator also vanishes at that point, as in the worked example above, the effective order is lower than the denominator's exponent suggests; always confirm the numerator does not cancel part of the singularity first.
  • Forgetting the 1(m−1)!\frac{1}{(m-1)!} factor or differentiating the wrong number of times. For a pole of order mm you differentiate m−1m-1 times and divide by (m−1)!(m-1)!, not mm times — an easy off-by-one slip.
  • Treating the residue of a removable singularity as nonzero. By definition it is always 00, since there is no principal part to contribute a 1/z1/z term.
Derivation AThe residue formula for a pole of order m

The goal is to derive, purely from the definition of the pole (a Laurent series whose principal part terminates after mm terms) and ordinary differentiation, the shortcut formula Res⁡z=z0f(z)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\operatorname{Res}_{z=z_0}f(z) = \frac{1}{(m-1)!}\lim_{z\to z_0}\frac{d^{m-1}}{dz^{m-1}}\big[(z-z_0)^m f(z)\big] that lets you extract a residue without ever writing out a full Laurent series by hand.

Step 1. Start from the definition of a pole of order mm: near z0z_0, f(z)=c−m(z−z0)m+c−m+1(z−z0)m−1+⋯+c−1z−z0+c0+c1(z−z0)+⋯f(z) = \dfrac{c_{-m}}{(z-z_0)^m} + \dfrac{c_{-m+1}}{(z-z_0)^{m-1}} + \cdots + \dfrac{c_{-1}}{z-z_0} + c_0 + c_1(z-z_0) + \cdots, a Laurent series whose principal part has exactly mm terms. The target coefficient is c−1c_{-1}.

Step 2. Multiply through by (z−z0)m(z-z_0)^m to clear every negative power at once: (z−z0)mf(z)=c−m+c−m+1(z−z0)+⋯+c−1(z−z0)m−1+c0(z−z0)m+⋯(z-z_0)^m f(z) = c_{-m} + c_{-m+1}(z-z_0) + \cdots + c_{-1}(z-z_0)^{m-1} + c_0(z-z_0)^m + \cdots. This is now an ordinary Taylor series in (z−z0)(z-z_0), analytic at z0z_0 itself, because every negative power has been absorbed.

Step 3. Recall how Taylor coefficients are extracted: if ϕ(z)=∑k≥0ak(z−z0)k\phi(z) = \sum_{k\ge 0} a_k (z-z_0)^k, then ak=ϕ(k)(z0)k!a_k = \dfrac{\phi^{(k)}(z_0)}{k!}. Apply this to ϕ(z):=(z−z0)mf(z)\phi(z) := (z-z_0)^m f(z) from Step 2, whose coefficient of (z−z0)m−1(z-z_0)^{m-1} is exactly c−1c_{-1}, the residue. So c−1=ϕ(m−1)(z0)(m−1)!c_{-1} = \dfrac{\phi^{(m-1)}(z_0)}{(m-1)!}.

Step 4. Write ϕ(m−1)(z0)\phi^{(m-1)}(z_0) as a limit, since ϕ\phi is analytic there and derivatives are continuous: ϕ(m−1)(z0)=lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\phi^{(m-1)}(z_0) = \lim_{z\to z_0}\dfrac{d^{m-1}}{dz^{m-1}}\big[(z-z_0)^m f(z)\big]. Combining with Step 3 gives the stated formula, Res⁡z=z0f(z)=1(m−1)!lim⁡z→z0dm−1dzm−1[(z−z0)mf(z)]\operatorname{Res}_{z=z_0}f(z) = \dfrac{1}{(m-1)!}\lim_{z\to z_0}\dfrac{d^{m-1}}{dz^{m-1}}\big[(z-z_0)^m f(z)\big]. Setting m=1m=1 collapses this to the simple-pole case, Res⁡z=z0f(z)=lim⁡z→z0(z−z0)f(z)\operatorname{Res}_{z=z_0}f(z) = \lim_{z\to z_0}(z-z_0)f(z), with no derivative needed at all — consistent with 0!=10! = 1 and a zeroth derivative being the function itself.

Ask the tutor: If distinguishing pole order from the numerator or denominator's structure trips you up, or you want to see the order-mm formula applied to a pole of order 33, ask the AI tutor for a few more practice problems.

How do I tell a pole of order 2 from two separate simple poles at the same point?

You cannot have two separate simple poles at the same point — a single point has exactly one Laurent series about it, and that series has one, unambiguous principal part. If the principal part's most negative surviving term is c−2/(z−z0)2c_{-2}/(z-z_0)^2, that is a single pole of order 22, full stop. The confusion usually comes from a denominator like (z−z0)2(z-z_0)^2 arising as a product of two distinct linear factors that happen to be equal, which is exactly what "order 22" already encodes — there is no separate case to worry about.

Why does the residue only care about the c minus 1 term and not the whole principal part?

Because that is exactly what survives contour integration. Integrating (z−z0)n(z-z_0)^n around a small closed loop about z0z_0 gives 00 for every integer nn except n=−1n=-1, where it gives 2πi2\pi i. When you integrate the full Laurent series term by term around such a loop, every term vanishes except the one multiplying c−1c_{-1}. So no matter how complicated the principal part is, only its 1/z1/z coefficient ever shows up in a contour integral — which is precisely why it earns the name "residue," what is left over after everything else cancels.

Do I always need the full order-m formula, or is there a shortcut for simple poles written as a fraction?

For a simple pole of f=p/qf = p/q with q(z0)=0q(z_0)=0 and q′(z0)≠0q'(z_0)\ne 0, you can skip the limit-of-a-product form entirely and use Res⁡z=z0f(z)=p(z0)/q′(z0)\operatorname{Res}_{z=z_0}f(z) = p(z_0)/q'(z_0), which follows from the order-11 case of the formula above by recognizing lim⁡z→z0(z−z0)/q(z)\lim_{z\to z_0}(z-z_0)/q(z) as 1/q′(z0)1/q'(z_0) from the definition of the derivative. It is faster whenever qq factors are awkward to isolate by hand, but it only applies at simple poles — for order 22 or higher you need the full derivative formula.

Before this

Where this leads