Evaluate
I=\int_{x=0}^{x=1} \int_{y=0}^{y=1} \int_{z=0}^{z=1} x+y+z \mathrm{~d} z \mathrm{~d} y \mathrm{~d} xWhat is meant by this expression is
I=\int_{x=0}^{x=1}\left(\int_{y=0}^{y=1}\left(\int_{z=0}^{z=1} x+y+z \mathrm{~d} z\right) \mathrm{d} y\right) \mathrm{d} xwhere, as before, the inner integral is performed first, integrating with respect to z, with x and y being treated as constants. So
\begin{aligned} I & =\int_{x=0}^{x=1}\left(\int_{y=0}^{y=1}\left[x z+y z+\frac{z^2}{2}\right]_0^1 \mathrm{~d} y\right) \mathrm{d} x \\ & =\int_{x=0}^{x=1}\left(\int_{y=0}^{y=1} x+y+\frac{1}{2} \mathrm{~d} y\right) \mathrm{d} x \\ & =\int_{x=0}^{x=1}\left[x y+\frac{y^2}{2}+\frac{1}{2} y\right]_0^1 \mathrm{~d} x \\ & =\int_{x=0}^{x=1} x+\frac{1}{2}+\frac{1}{2} \mathrm{~d} x \\ & =\left[\frac{x^2}{2}+x\right]_0^1 \\ & =\frac{3}{2} \end{aligned}Consideration of the limits of integration shows that the integral is evaluated over a unit cube.