Triple integrals

Remember how double integrals can be written as iterated integrals. Triple integrals are essentially the same thing as double integrals. (We just add a third dimension.) We will turn triple integrals into (triple) iterated integrals.

Just as with double integrals, the only trick is determining the limits on the iterated integrals. (Unfortunately, it's harder to draw in three dimensions.)

Defined by Riemann sums

Let f (x, y, z) be the density of some three-dimensional solid W. We want to define the triple integral of f over W to be the total mass of W.

As with double integrals, we define the integral with Riemann sums. We chop up the solid W into small boxes, say with dimensions $\Delta$x, $\Delta$y, $\Delta$z. If W happened to be a cube, this chopping might look something like this.

The volume of each small box is

$$\Delta \Delta \Delta \Delta$$

Think of the boxes as being arranged in layers, with each layer arranged into rows and columns. We can then index the boxes so that box ijk refers to the box in the ith row, the jth column, and the kth layer.

For each box, we pick a point in the box to represent that box. For box ijk, we call that point (x_ijk, y_ijk, z_ijk). Pretend that the density of box ijk is constant, i.e., that the density is f (x_ijk, y_ijk, z_ijk) everywhere in that box. The mass of box ijk is its density times it volume:

$$\Delta$$

We sum up these approximate masses to estimate the total mass of the solid W. We obtain the Riemann sum

$$\sum_{{ijk}}^{} \Delta$$

where the sum is over all small boxes.

Let $\Delta$x$\to$ 0, $\Delta$y$\to$ 0, and $\Delta$z$\to$ 0, (and the number of small boxes go to infinity). The Riemann sum approaches the triple integral over the solid W,

$$\iiint_{W}^{} \lim_{{\Delta x, \Delta y, \Delta z \to 0}}^{} \sum_{{ijk}}^{} \Delta$$

assuming f is continuous. The triple integral is the actual mass of W.

Triple iterated integrals

If the solid W is a cube defined by a$\le$x$\le$b, c$\le$y$\le$d, and p$\le$z$\le$q, then we can easily write the triple integral as an iterated integral. We could first integrate x from a to b, then integrate y from c to d, and finally integrate z from p to q,

$$\iiint_{W}^{} \int_{p}^{q} \left(\vphantom{\int_c^d \left(\int_a^b f(x,y,z) dx \right ) dy }\right. \int_{c}^{d} \left(\vphantom{\int_a^b f(x,y,z) dx }\right. \int_{a}^{b} \left.\vphantom{\int_a^b f(x,y,z) dx }\right) \left.\vphantom{\int_c^d \left(\int_a^b f(x,y,z) dx \right ) dy }\right)$$

This order of integration corresponds to taking the Riemann sum and first summing over rows i, then summing over columns j, and finally summing over layers k.

Just as with double integrals, other orders of integration are possible. We could, for example, first integrate z, then integrate x, and lastly integrate y,

$$\iiint_{W}^{} \int_{c}^{d} \left(\vphantom{\int_a^b \left(\int_p^q f(x,y,z) dz \right ) dx }\right. \int_{a}^{b} \left(\vphantom{\int_p^q f(x,y,z) dz }\right. \int_{p}^{q} \left.\vphantom{\int_p^q f(x,y,z) dz }\right) \left.\vphantom{\int_a^b \left(\int_p^q f(x,y,z) dz \right ) dx }\right)$$

The iterated integral is simple when the solid W is a rectangular solid (like a cube, but where all edges aren't necessarily the same length). For more complicated shapes, finding the limits of integration can be tricky. We'll do some examples in class.

For now, just remember these rules, which are analogous to the rules we had for limits on double iterated integrals.

1. The outer limits have to be constant. They cannot depend on any of the variables.
2. The middle limits can depend on the variable from the outer integral only. They cannot depend on the variable from the inner integral.
3. The inner limits can depend on the variable from the outer integral and the variable from the middle integral.

So, for example, the following integral makes sense

$$\iiint_{W}^{} \int_{2}^{3} \left(\vphantom{\int_{1-z}^0 \left(\int_{-y^2-z^2}^{y^2+z^2} f(x,y,z) dx \right ) dy }\right. \int_{{1-z}}^{0} \left(\vphantom{\int_{-y^2-z^2}^{y^2+z^2} f(x,y,z) dx }\right. \int_{{-y^2-z^2}}^{{y^2+z^2}} \left.\vphantom{\int_{-y^2-z^2}^{y^2+z^2} f(x,y,z) dx }\right) \left.\vphantom{\int_{1-z}^0 \left(\int_{-y^2-z^2}^{y^2+z^2} f(x,y,z) dx \right ) dy }\right)$$

It describes the integral of f over the region W defined by

$$\le \le$$

$$\le \le$$

$$\le \le$$

(not that you should be able to visualize what W looks like).

The following integral does not make sense

$$\iiint_{W}^{} \int_{y}^{x} \left(\vphantom{\int_{1}^{2x} \left(\int_0^1 f(x,y,z) dx \right ) dy }\right. \int_{{1}}^{{2x}} \left(\vphantom{\int_0^1 f(x,y,z) dx }\right. \int_{0}^{1} \left.\vphantom{\int_0^1 f(x,y,z) dx }\right) \left.\vphantom{\int_{1}^{2x} \left(\int_0^1 f(x,y,z) dx \right ) dy }\right)$$

Can you see why? The outer integral depends on both x and y (but y isn't defined until you go inside the middle integral, and x isn't defined until you go inside the inner integral). Also, the middle integral depends on x.

Click here for triple integral examples.