Fubini theorems

First, I want to note some notation that is a little more precise than the notation in the textbook. The notation in the textbook is very common, and it's usually quite clear what it means, but it's not good enough if you want to be completely unambiguous about what variables you're using and where. So rather than write, for example, ∫∫_D f(x, y) dA, where D is a region in 2 dimensions (formally a relation between 2 variables) and f is a function of 2 variables, I'll write ∫_(x,y)∈D f(x, y) |dx ∧ dy|. So to begin with, since the integrand (all of the stuff after the integral symbols) makes it clear that there are 2 variables of integration, it's not necessary to repeat the integral symbol. But just as we write f(x, y) after that symbol to indicate the value of the function f at particular values of the variables x and y (rather than its value somewhere else), so I write (x, y) ∈ D beneath that symbol to indicate that the point whose coordinates are those values (rather than some other point) belongs to the region D.

At the end, since dA (and dV in 3 dimensions) don't indicate which variables are being used, I use the notation |dx ∧ dy| (or |dx ∧ dy ∧ dz| in 3 dimensions). This notation is more complicated than strictly necessary, but there is a reason for it; just as the notation dy/dx for a derivative is not merely an arbitrary symbol but can be literally understood as the result of dividing expressions (called differentials) obtained by applying an operator d, so the notation |dx ∧ dy| for an area element (or |dx ∧ dy ∧ dz| for a volume element) is not merely an arbitrary symbol but can be literally understood as the absolute value of an expression (called an exterior differential form) involving an operator ∧. However, don't worry about that for now; just treat it as a notation used to indicate precisely which variables are used in the area (or volume) element.

Using that notation, here are the important theorems:

• The integral of a continuous function on a compact (that is closed and bounded) region always exists: ∫_(x,y)∈D f(x, y) |dx ∧ dy| exists if f is continuous and D is compact (and similarly in more variables).
• If two regions D₁ and D₂ are completely disjoint (no overlap at all), or if their overlap is contained within a single point/line/plane of fewer dimensions than the overall number of variables, and if a function f has integrals on both of these regions, then the integral of f on their union (the combined region D₁ ∪ D₂) also exists and is the sum of the separate integrals: ∫_{(x,y)∈D1∪D2} f(x, y) |dx ∧ dy| = ∫_(x,y)∈D1 f(x, y) |dx ∧ dy| + ∫_(x,y)∈D2 f(x, y) |dx ∧ dy| (and similarly in more variables) if the integrals on the right exist and the overlap is small.
• In any double (or higher) integral, if two of the variables are swapped in both the function being integrated and in the region over which it is integrated, then the result is the same (so that if either integral exists, then so does the other, and then they are equal): ∫_(x,y)∈D f(x, y) |dx ∧ dy| = ∫_(y,x)∈D f(y, x) |dx ∧ dy| (note no swap in the area element at the end).
• For a region D in 2 dimensions, if there are constants a and b and continuous functions g and h of 1 variable each such that (x, y) ∈ D if and only if a ≤ x≤ b and g(x) ≤ y ≤ h(x), and if g(x) ≤ h(x) whenever a ≤ x ≤ b, then the integral of any continuous function f on D is the same as the corresponding iterated integral: ∫_(x,y)∈D f(x, y) |dx ∧ dy| = ∫^b_x=a (∫^h(x)_y=g(x) dy) dx.
• For a region D in 3 (or more) variables, if there are a compact region R in 2 variables (or in general a compact region of one fewer dimension) and continuous functions g and h of 2 variables each (or in general with the same number of variables as R has dimensions) such that (x, y, z) ∈ D if and only if (x, y) ∈ R and g(x, y) ≤ z ≤ h(x, y), and if g(x, y) ≤ h(x, y) whenever (x, y) ∈ R, then the integral of any continuous function f on D is the same as the corresponding iterated integral: ∫_(x,y,z)∈D f(x, y, z) |dx ∧ dy ∧ dz| = ∫_(x,y)∈R (∫^h(x,y)_z=g(x,y) f(x, y, z) dz) |dx ∧ dy|.

By itself, the Fubini Theorem only works for regions of particular shapes, but the other theorems combine to make it more useful. First of all, the middle theorem allows us to put the variables in whatever order we like. Even so, the regions still require particular shapes; we can just orient those however we wish. The second theorem, in principle, allows us to divide a region up into smaller regions appropriate for the Fubini Theorem; the only question is whether the integrals exist. The first theorem guarantees this existence for continuous functions.

So using these in order, if you want to integrate over a crazy region, then divide the region into pieces of suitable shape. If the function is continuous and these smaller regions are all compact, then you know that their integrals exist; and if the regions overlap only slightly, then you can recover the answer to the original problem by adding them up. Finally, to get the integrals on these small regions, think of the variables as coming in whichever order works best, and use the Fubini Theorem (possibly more than once) to replace double and triple integrals with iterated integrals. Hopefully, these will be integrals that you can do!

The permanent URI of this web page is http://tobybartels.name/MATH-2080/2016SP/Fubini/.