Double Series and Fubini

Numbers on a grid get summed by rows first, then by columns, in this scene. With absolute convergence, the two answers agree.

Add along the rows or down the columns: the result is one and the same. The order cannot change it because the absolute values sum to a finite total. With only conditional convergence, swapping the order changes the answer itself. The headline value is the second formula below.

For example, stack πŸ“¦ boxes in a grid and count the items inside each box. Count row by row or column by column β€” the total count comes out identical.

βˆ‘mβˆ‘nam,n=βˆ‘nβˆ‘mam,n=sup⁑finiteΒ Fβˆ‘Fam,n\sum_{m}\sum_{n} a_{m,n} = \sum_{n}\sum_{m} a_{m,n} = \sup_{\text{finite } F} \sum_{F} a_{m,n}
βˆ‘m,nβ‰₯11m2n2=Ο€436β‰ˆ2.7058\sum_{m,n \ge 1} \frac{1}{m^2n^2} = \frac{\pi^4}{36} \approx 2.7058

Examples & Exercises

Prove Fubini for infinite sums (nonnegative terms) completely: if am,nβ‰₯0a_{m,n} \ge 0 then βˆ‘mβˆ‘nam,n=βˆ‘nβˆ‘mam,n\sum_m \sum_n a_{m,n} = \sum_n \sum_m a_{m,n}.

Show both sides equal the supremum SS of sums over finite rectangles.

S=sup⁑M,Nβˆ‘m≀Mβˆ‘n≀Nam,n∈[0,∞]S = \sup_{M,N} \sum_{m \le M} \sum_{n \le N} a_{m,n} \in [0, \infty]
βˆ‘m=1βˆžβˆ‘n=1∞am,n=S=βˆ‘n=1βˆžβˆ‘m=1∞am,n\sum_{m=1}^{\infty}\sum_{n=1}^{\infty} a_{m,n} = S = \sum_{n=1}^{\infty}\sum_{m=1}^{\infty} a_{m,n}
N=5N = 5: sum by rows β‰ˆ2.1422\approx 2.1422, sum by columns β‰ˆ2.1422\approx 2.1422 β€” both head to (Ο€2/6)2=Ο€4/36β‰ˆ2.7058(\pi^2/6)^2 = \pi^4/36 \approx 2.7058.
Show proof
  1. Put T(M,N)=βˆ‘m≀Mβˆ‘n≀Nam,nT(M,N) = \sum_{m \le M}\sum_{n \le N} a_{m,n} (a finite sum, β‰₯0\ge 0) and S=sup⁑M,NT(M,N)∈[0,∞]S = \sup_{M,N} T(M,N) \in [0,\infty].

    TT is monotone in each of M,NM, N.

    For each fixed mm, the inner series βˆ‘n=1∞am,n∈[0,∞]\sum_{n=1}^{\infty} a_{m,n} \in [0,\infty] converges (finitely or to +∞+\infty) since the terms are nonnegative.

  2. Put R=βˆ‘m=1βˆžβˆ‘n=1∞am,nR = \sum_{m=1}^{\infty}\sum_{n=1}^{\infty} a_{m,n}.

    For each M,NM, N: T(M,N)=βˆ‘m≀Mβˆ‘n≀Nam,nβ‰€βˆ‘m≀Mβˆ‘n=1∞am,n≀RT(M,N) = \sum_{m \le M}\sum_{n \le N} a_{m,n} \le \sum_{m \le M}\sum_{n=1}^{\infty} a_{m,n} \le R (inner finite sum ≀\le inner infinite sum, outer finite ≀\le outer infinite β€” all terms nonnegative).

    Taking the supremum over M,NM, N gives S≀RS \le R.

  3. The other direction: for each MM, βˆ‘m≀Mβˆ‘n=1∞am,n=sup⁑NT(M,N)≀S\sum_{m \le M}\sum_{n=1}^{\infty} a_{m,n} = \sup_N T(M,N) \le S (the inner infinite sum is the supremum of the inner finite sums).

    This equality needs justification: for fixed MM, each inner sum βˆ‘n≀Nam,n\sum_{n \le N} a_{m,n} is increasing in NN (since am,nβ‰₯0a_{m,n} \ge 0) and converges to βˆ‘n=1∞am,n\sum_{n=1}^{\infty} a_{m,n}.

    For a finite sum (MM terms), supremum and summation interchange β€” a finite sum of increasing sequences is increasing, and its supremum equals the sum of the suprema (checked termwise).

    Hence βˆ‘m≀Msup⁑Nβˆ‘n≀Nam,n=sup⁑Nβˆ‘m≀Mβˆ‘n≀Nam,n=sup⁑NT(M,N)\sum_{m \le M} \sup_N \sum_{n \le N} a_{m,n} = \sup_N \sum_{m \le M} \sum_{n \le N} a_{m,n} = \sup_N T(M,N).

    Taking the supremum over MM gives R=sup⁑Mβˆ‘m≀Mβˆ‘n=1∞am,n≀SR = \sup_M \sum_{m \le M}\sum_{n=1}^{\infty} a_{m,n} \le S.

    Hence R=SR = S.

  4. Swapping the roles of rows and columns, the identical argument gives C=βˆ‘n=1βˆžβˆ‘m=1∞am,n=SC = \sum_{n=1}^{\infty}\sum_{m=1}^{\infty} a_{m,n} = S (the same supremum of the same T(M,N)T(M,N)).

    Hence R=C=SR = C = S β€” the two iterated sums agree.

    If S<∞S < \infty the common value is finite; if S=∞S = \infty both sides diverge to +∞+\infty.

    QED.

Extend to the absolutely convergent case (short complete proof): if βˆ‘m,n∣am,n∣<∞\sum_{m,n} \lvert a_{m,n}\rvert < \infty, use the decomposition a=a+βˆ’aβˆ’a = a^+ - a^- to apply the first problem to each part and conclude βˆ‘mβˆ‘nam,n=βˆ‘nβˆ‘mam,n\sum_m\sum_n a_{m,n} = \sum_n\sum_m a_{m,n}.

a+=max⁑(a,0),aβˆ’=max⁑(βˆ’a,0),∣a∣=a++aβˆ’,a=a+βˆ’aβˆ’a^+ = \max(a, 0), \quad a^- = \max(-a, 0), \quad \lvert a\rvert = a^+ + a^-, \quad a = a^+ - a^-
βˆ‘mβˆ‘na=βˆ‘mβˆ‘na+βˆ’βˆ‘mβˆ‘naβˆ’=βˆ‘nβˆ‘ma+βˆ’βˆ‘nβˆ‘maβˆ’=βˆ‘nβˆ‘ma\sum_m\sum_n a = \sum_m\sum_n a^+ - \sum_m\sum_n a^- = \sum_n\sum_m a^+ - \sum_n\sum_m a^- = \sum_n\sum_m a
N=5N = 5: signs are present but βˆ‘βˆ£a∣\sum \lvert a\rvert converges, so sum by rows β‰ˆ0.7033\approx 0.7033 and sum by columns β‰ˆ0.7033\approx 0.7033 agree, both heading to Ο€4/144β‰ˆ0.6765\pi^4/144 \approx 0.6765.
Show proof
  1. Put am,n+=max⁑(am,n,0)a^+_{m,n} = \max(a_{m,n}, 0) and am,nβˆ’=max⁑(βˆ’am,n,0)a^-_{m,n} = \max(-a_{m,n}, 0).

    Then 0≀am,nΒ±β‰€βˆ£am,n∣0 \le a^{\pm}_{m,n} \le \lvert a_{m,n}\rvert, with ∣a∣=a++aβˆ’\lvert a\rvert = a^+ + a^- and a=a+βˆ’aβˆ’a = a^+ - a^-.

    By comparison, both βˆ‘am,n+\sum a^+_{m,n} and βˆ‘am,nβˆ’\sum a^-_{m,n} converge, bounded by βˆ‘βˆ£am,n∣<∞\sum \lvert a_{m,n}\rvert < \infty.

  2. Since a+a^+ and aβˆ’a^- are nonnegative, apply the first problem (Fubini, nonnegative terms) to each: βˆ‘mβˆ‘na+=βˆ‘nβˆ‘ma+<∞\sum_m\sum_n a^+ = \sum_n\sum_m a^+ < \infty, and likewise for aβˆ’a^-.

    All iterated sums are finite, so subtraction is legitimate, and βˆ‘mβˆ‘na=βˆ‘mβˆ‘na+βˆ’βˆ‘mβˆ‘naβˆ’=βˆ‘nβˆ‘ma+βˆ’βˆ‘nβˆ‘maβˆ’=βˆ‘nβˆ‘ma\sum_m\sum_n a = \sum_m\sum_n a^+ - \sum_m\sum_n a^- = \sum_n\sum_m a^+ - \sum_n\sum_m a^- = \sum_n\sum_m a.

    QED.

Show the hypothesis is sharp: with an,n=1a_{n,n} = 1, an,n+1=βˆ’1a_{n,n+1} = -1 (nβ‰₯1n \ge 1) and 00 elsewhere, verify completely that every row sum is 00 (so βˆ‘mβˆ‘nam,n=0\sum_m\sum_n a_{m,n} = 0) but column 11 sums to 11 and columns nβ‰₯2n \ge 2 sum to 00 (so βˆ‘nβˆ‘mam,n=1\sum_n\sum_m a_{m,n} = 1), and note βˆ‘βˆ£a∣=∞\sum \lvert a\rvert = \infty.

an,n=1,an,n+1=βˆ’1Β (nβ‰₯1),elseΒ 0a_{n,n} = 1, \quad a_{n,n+1} = -1 \ (n \ge 1), \quad \text{else } 0
βˆ‘mβˆ‘nam,n=βˆ‘m0=0β‰ 1=1+βˆ‘nβ‰₯20=βˆ‘nβˆ‘mam,n\sum_m \sum_n a_{m,n} = \sum_m 0 = 0 \quad\ne\quad 1 = 1 + \sum_{n \ge 2} 0 = \sum_n \sum_m a_{m,n}
The alternating series in the picture above is absolutely convergent, so reordering is safe β€” but as the counterexample below shows, with βˆ‘βˆ£a∣=∞\sum\lvert a\rvert = \infty the order can change the answer.
Show proof
  1. Fix a row (mm fixed).

    The only nonzero entries are am,m=1a_{m,m} = 1 at n=mn = m and am,m+1=βˆ’1a_{m,m+1} = -1 at n=m+1n = m+1.

    Hence the row sum is βˆ‘n=1∞am,n=1+(βˆ’1)=0\sum_{n=1}^{\infty} a_{m,n} = 1 + (-1) = 0 (a finite sum, so no convergence issue).

    This holds for every mm, so βˆ‘m=1βˆžβˆ‘n=1∞am,n=βˆ‘m=1∞0=0\sum_{m=1}^{\infty}\sum_{n=1}^{\infty} a_{m,n} = \sum_{m=1}^{\infty} 0 = 0.

  2. Fix a column (nn fixed).

    The nonzero entries: an,n=1a_{n,n} = 1 (row m=nm = n) and, for nβ‰₯2n \ge 2, anβˆ’1,n=βˆ’1a_{n-1,n} = -1 (row m=nβˆ’1m = n-1, from am,m+1=βˆ’1a_{m,m+1} = -1 with m=nβˆ’1m = n-1).

    Column 11: only a1,1=1a_{1,1} = 1, so the sum is 11.

    Column nβ‰₯2n \ge 2: 1+(βˆ’1)=01 + (-1) = 0.

    Hence βˆ‘n=1βˆžβˆ‘m=1∞am,n=1+βˆ‘n=2∞0=1\sum_{n=1}^{\infty}\sum_{m=1}^{\infty} a_{m,n} = 1 + \sum_{n=2}^{\infty} 0 = 1.

  3. The two iterated sums differ, 0β‰ 10 \ne 1 β€” the order changes the answer.

    The broken hypothesis is absolute convergence: βˆ‘m,n∣am,n∣=βˆ‘n=1∞(1+1)=∞\sum_{m,n} \lvert a_{m,n}\rvert = \sum_{n=1}^{\infty} (1 + 1) = \infty, so the second problem's hypothesis fails.

    Without the nonnegative (or absolutely convergent) hypothesis Fubini fails, and this example shows the hypothesis is sharp.

    QED.

7. Series