Triple Integrals

Work in progress. This note is still being written and incomplete.

Let GG be a bounded solid region and ff a real-valued function defined on GG. Partition a rectangular box containing GG into sub-boxes by planes parallel to the coordinate planes, and choose a sample point (xk∗,yk∗,zk∗)(x_k^*, y_k^*, z_k^*) in each resulting subregion of volume ΔVk\Delta V_k. The triple integral of ff over GG is

∭Gf dV=lim⁡n→∞∑k=1nf(xk∗,yk∗,zk∗) ΔVk\iiint_G f\,\text{d}V = \lim_{n \to \infty} \sum_{k=1}^{n} f(x_k^*, y_k^*, z_k^*)\,\Delta V_k

provided the limit exists.

Properties

Suppose ff and gg integrable over GG, and c∈Rc \in \mathbb{R}:

  • ∭G1 dV\displaystyle \iiint_G 1\,\text{d}V: volume of GG
  • ∭Gf dV\displaystyle\iiint_G f\,\text{d}V: mass of GG with ff density
  • ∭Gcf dV=c∭Gf dV\displaystyle\iiint_G cf\,\text{d}V = c\iiint_G f\,\text{d}V
  • ∭G(f+g) dV=∭Gf dV+∭Gg dV\displaystyle\iiint_G (f+g)\,\text{d}V = \iiint_G f\,\text{d}V + \iiint_G g\,\text{d}V
  • If G=G1∪G2G = G_1 \cup G_2 with G1G_1, G2G_2 having disjoint interiors,
    ∭Gf dV=∭G1f dV+∭G2f dV\displaystyle\iiint_G f\,\text{d}V = \iiint_{G_1} f\,\text{d}V + \iiint_{G_2} f\,\text{d}V

Fubini’s Theorem

Let GG be the rectangular box a≤x≤ba \le x \le b, c≤y≤dc \le y \le d, k≤z≤lk \le z \le l.

If ff is continuous on GG:

∭Gf(x,y,z) dV=∫ab∫cd∫klf(x,y,z) dz dy dx\iiint_G f(x,y,z)\,\text{d}V = \int_a^b \int_c^d \int_k^l f(x,y,z)\,\text{d}z\,\text{d}y\,\text{d}x

The box is cut into thin slabs perpendicular to the xx-axis. Each slab is cut into columns parallel to the zz-axis. The inner integral adds up ff along 1 column. The middle integral adds up the columns across 1 slab. The outer integral adds up the slabs from x=ax = a to x=bx = b.

Integrating f along a column gives its total. Integrating across the columns of a slab gives the slab total. Integrating across the slabs gives the integral over the box.

If f(x,y,z)=p(x) q(y) r(z)f(x,y,z) = p(x)\,q(y)\,r(z) is a product of single-variable continuous functions, the iterated integral splits into 3 single integrals. Each inner integral pulls the factors that don’t depend on its variable outside, and the limits are constants because GG is a box.

∭Gp(x) q(y) r(z) dV=∫abp(x) dx⋅∫cdq(y) dy⋅∫klr(z) dz\iiint_G p(x)\,q(y)\,r(z)\,\text{d}V = \int_a^b p(x)\,\text{d}x \cdot \int_c^d q(y)\,\text{d}y \cdot \int_k^l r(z)\,\text{d}z

Regions of Integration

A solid region GG is a simple xyxy-solid region iff:

  • Bounded above by z=v(x,y)z = v(x,y) and below by z=u(x,y)z = u(x,y)
  • uu, vv continuous on the projection RR of GG onto the xyxy-plane
  • u(x,y)≤v(x,y)u(x,y) \le v(x,y) for all (x,y)∈R(x,y) \in R
A solid bounded above by z = v(x, y) and below by z = u(x, y) over its projection R onto the xy-plane. Each vertical line through R enters at u and leaves at v. Drag to rotate.

For ff continuous on GG,

∭Gf(x,y,z) dV=∬R(∫u(x,y)v(x,y)f(x,y,z) dz)dA\iiint_G f(x,y,z)\,\text{d}V = \iint_R \left( \int_{u(x,y)}^{v(x,y)} f(x,y,z)\,\text{d}z \right) \text{d}A

The projection RR can be:

  • Type I, or vertically simple
    R={(x,y):a≤x≤b, g1(x)≤y≤g2(x)}R = \{(x,y) : a \le x \le b,\ g_1(x) \le y \le g_2(x)\}
  • Type II, or horizontally simple
    R={(x,y):c≤y≤d, h1(y)≤x≤h2(y)}R = \{(x,y) : c \le y \le d,\ h_1(y) \le x \le h_2(y)\}

If RR is Type I,

∭Gf(x,y,z) dV=∫ab∫g1(x)g2(x)∫u(x,y)v(x,y)f(x,y,z) dz dy dx\iiint_G f(x,y,z)\,\text{d}V = \int_a^b \int_{g_1(x)}^{g_2(x)} \int_{u(x,y)}^{v(x,y)} f(x,y,z)\,\text{d}z\,\text{d}y\,\text{d}x
A Type I region R . Each slice at a fixed x runs in y from g₁(x) to g₂(x). The solid above the slice forms a cross-section. Drag to rotate.

If RR is Type II,

∭Gf(x,y,z) dV=∫cd∫h1(y)h2(y)∫u(x,y)v(x,y)f(x,y,z) dz dx dy\iiint_G f(x,y,z)\,\text{d}V = \int_c^d \int_{h_1(y)}^{h_2(y)} \int_{u(x,y)}^{v(x,y)} f(x,y,z)\,\text{d}z\,\text{d}x\,\text{d}y

The region is similar to type 1 but the axes switched.

Change of Variables

Let TT be a transformation from uvwuvw-space to xyzxyz-space mapping a region SS onto a region RR. If:

  • TT is one-to-one on SS
  • TT has continuous first partial derivatives
  • The Jacobian ∂(x,y,z)∂(u,v,w)\dfrac{\partial(x,y,z)}{\partial(u,v,w)} is non-zero and does not change sign on SS

then

∭Rf(x,y,z) dV=∭Sf(x(u,v,w),y(u,v,w),z(u,v,w))∣∂(x,y,z)∂(u,v,w)∣dV\iiint_R f(x,y,z)\,\text{d}V = \iiint_S f(x(u,v,w), y(u,v,w), z(u,v,w)) \left| \frac{\partial(x,y,z)}{\partial(u,v,w)} \right| \text{d}V

Cylindrical Coordinates

In cylindrical coordinates, the volume element is

dV=r dz dr dθ\text{d}V = r\,\text{d}z\,\text{d}r\,\text{d}\theta

Here rr is the Jacobian ∂(x,y,z)∂(r,θ,z)\dfrac{\partial(x,y,z)}{\partial(r,\theta,z)} of the change of variables.

Let GG have upper surface z=v(r,θ)z = v(r,\theta) and lower surface z=u(r,θ)z = u(r,\theta), with projection onto the xyxy-plane a simple polar region R′R' in the rθr\theta-plane. For ff continuous on GG,

∭Gf(x,y,z) dV=∬R′(∫u(r,θ)v(r,θ)f(rcos⁡θ,rsin⁡θ,z) dz)r dA\iiint_G f(x,y,z)\,\text{d}V = \iint_{R'} \left( \int_{u(r,\theta)}^{v(r,\theta)} f(r\cos\theta, r\sin\theta, z)\,\text{d}z \right) r\,\text{d}A

If for each fixed θ∈[α,β]\theta \in [\alpha, \beta], rr varies from r1(θ)r_1(\theta) to r2(θ)r_2(\theta),

∭Gf(x,y,z) dV=∫αβ∫r1(θ)r2(θ)∫u(r,θ)v(r,θ)f(rcos⁡θ,rsin⁡θ,z) r dz dr dθ\iiint_G f(x,y,z)\,\text{d}V = \int_\alpha^\beta \int_{r_1(\theta)}^{r_2(\theta)} \int_{u(r,\theta)}^{v(r,\theta)} f(r\cos\theta, r\sin\theta, z)\,r\,\text{d}z\,\text{d}r\,\text{d}\theta

Spherical Coordinates

In spherical coordinates, the volume element is

dV=ρ2sin⁡ϕ dρ dϕ dθ\text{d}V = \rho^2 \sin\phi\,\text{d}\rho\,\text{d}\phi\,\text{d}\theta

Here ρ2sin⁡ϕ\rho^2\sin\phi is the Jacobian ∂(x,y,z)∂(ρ,ϕ,θ)\dfrac{\partial(x,y,z)}{\partial(\rho,\phi,\theta)} of the change of variables.

For ff continuous on GG, with appropriate limits (see common spherical regions),

∭Gf(x,y,z) dV=∭f(ρsin⁡ϕcos⁡θ,ρsin⁡ϕsin⁡θ,ρcos⁡ϕ) ρ2sin⁡ϕ dρ dϕ dθ\iiint_G f(x,y,z)\,\text{d}V = \iiint f(\rho\sin\phi\cos\theta, \rho\sin\phi\sin\theta, \rho\cos\phi)\,\rho^2\sin\phi\,\text{d}\rho\,\text{d}\phi\,\text{d}\theta

Worked Examples

Box

Evaluate ∭Gxyz dV\iiint_G xyz\,\text{d}V over G:0≤x≤1, 0≤y≤2, 0≤z≤3G: 0 \le x \le 1,\ 0 \le y \le 2,\ 0 \le z \le 3.

The integrand xyzxyz is a product of single-variable functions and GG is a box, so the integral separates into 3 single integrals.

∫01x dx⋅∫02y dy⋅∫03z dz=12⋅2⋅92=92\int_0^1 x\,\text{d}x \cdot \int_0^2 y\,\text{d}y \cdot \int_0^3 z\,\text{d}z = \frac{1}{2} \cdot 2 \cdot \frac{9}{2} = \frac{9}{2}

Tetrahedron

Find the volume of the solid bounded by the planes x=0x = 0, y=0y = 0, z=0z = 0 and x+y+z=1x + y + z = 1.

The solid lies between z=0z = 0 and z=1−x−yz = 1 - x - y. The projection onto the xyxy-plane is the triangle 0≤x≤10 \le x \le 1, 0≤y≤1−x0 \le y \le 1 - x.

V=∫01∫01−x∫01−x−ydz dy dxV = \int_0^1 \int_0^{1-x} \int_0^{1-x-y} \text{d}z\,\text{d}y\,\text{d}x

Integrate with respect to zz.

V=∫01∫01−x(1−x−y) dy dxV = \int_0^1 \int_0^{1-x} (1 - x - y)\,\text{d}y\,\text{d}x

Integrate with respect to yy.

V=∫01[(1−x)y−y22]01−xdx=∫01(1−x)22 dxV = \int_0^1 \left[ (1-x)y - \frac{y^2}{2} \right]_0^{1-x} \text{d}x = \int_0^1 \frac{(1-x)^2}{2}\,\text{d}x

Integrate with respect to xx.

V=[−(1−x)36]01=16V = \left[ -\frac{(1-x)^3}{6} \right]_0^1 = \frac{1}{6}

Paraboloid

Find the volume of the solid between the plane z=0z = 0 and the paraboloid z=4−x2−y2z = 4 - x^2 - y^2.

The paraboloid meets z=0z = 0 on the circle x2+y2=4x^2 + y^2 = 4, so the projection of the solid onto the xyxy-plane is a disk. In Cartesian coordinates, the disk gives limits with square roots.

V=∫−22∫−4−x24−x2∫04−x2−y2dz dy dxV = \int_{-2}^{2} \int_{-\sqrt{4-x^2}}^{\sqrt{4-x^2}} \int_0^{4 - x^2 - y^2} \text{d}z\,\text{d}y\,\text{d}x

The surface and the disk both depend on xx and yy only through x2+y2x^2 + y^2, which is r2r^2 in cylindrical coordinates. Cylindrical coordinates turn the disk into the rectangle 0≤r≤20 \le r \le 2, 0≤θ≤2π0 \le \theta \le 2\pi and the height into z=4−r2z = 4 - r^2, so every limit is simple. The volume element is dV=r dz dr dθ\text{d}V = r\,\text{d}z\,\text{d}r\,\text{d}\theta.

V=∫02π∫02∫04−r2r dz dr dθV = \int_0^{2\pi} \int_0^2 \int_0^{4 - r^2} r\,\text{d}z\,\text{d}r\,\text{d}\theta

Integrate with respect to zz.

V=∫02π∫02(4r−r3) dr dθV = \int_0^{2\pi} \int_0^2 (4r - r^3)\,\text{d}r\,\text{d}\theta

Integrate with respect to rr.

V=∫02π[2r2−r44]02dθ=∫02π4 dθ=8πV = \int_0^{2\pi} \left[ 2r^2 - \frac{r^4}{4} \right]_0^2 \text{d}\theta = \int_0^{2\pi} 4\,\text{d}\theta = 8\pi

Ball

Evaluate ∭G(x2+y2+z2) dV\iiint_G (x^2 + y^2 + z^2)\,\text{d}V over the unit ball x2+y2+z2≤1x^2 + y^2 + z^2 \le 1.

In Cartesian coordinates, the ball gives limits with nested square roots.

∫−11∫−1−x21−x2∫−1−x2−y21−x2−y2(x2+y2+z2) dz dy dx\int_{-1}^{1} \int_{-\sqrt{1-x^2}}^{\sqrt{1-x^2}} \int_{-\sqrt{1-x^2-y^2}}^{\sqrt{1-x^2-y^2}} (x^2 + y^2 + z^2)\,\text{d}z\,\text{d}y\,\text{d}x

The region and the integrand both depend on xx, yy and zz only through x2+y2+z2=ρ2x^2 + y^2 + z^2 = \rho^2, the squared distance from the origin. Spherical coordinates turn the ball into the box 0≤ρ≤10 \le \rho \le 1, 0≤ϕ≤π0 \le \phi \le \pi, 0≤θ≤2π0 \le \theta \le 2\pi and the integrand into ρ2\rho^2, so every limit is constant. The volume element is dV=ρ2sin⁡ϕ dρ dϕ dθ\text{d}V = \rho^2 \sin\phi\,\text{d}\rho\,\text{d}\phi\,\text{d}\theta.

∭G(x2+y2+z2) dV=∫02π∫0π∫01ρ2⋅ρ2sin⁡ϕ dρ dϕ dθ\iiint_G (x^2 + y^2 + z^2)\,\text{d}V = \int_0^{2\pi} \int_0^\pi \int_0^1 \rho^2 \cdot \rho^2 \sin\phi\,\text{d}\rho\,\text{d}\phi\,\text{d}\theta

The integral separates.

∫02πdθ⋅∫0πsin⁡ϕ dϕ⋅∫01ρ4 dρ=2π⋅2⋅15=4π5\int_0^{2\pi} \text{d}\theta \cdot \int_0^\pi \sin\phi\,\text{d}\phi \cdot \int_0^1 \rho^4\,\text{d}\rho = 2\pi \cdot 2 \cdot \frac{1}{5} = \frac{4\pi}{5}

Cone

Find the volume of the solid enclosed laterally by the cone ϕ=π4\phi = \frac{\pi}{4} and above by the plane z=1z = 1.

From the last region in common spherical regions with a=1a = 1, the limits are 0≤ρ≤sec⁡ϕ0 \le \rho \le \sec\phi, 0≤ϕ≤π40 \le \phi \le \frac{\pi}{4}, 0≤θ≤2π0 \le \theta \le 2\pi.

V=∫02π∫0π/4∫0sec⁡ϕρ2sin⁡ϕ dρ dϕ dθV = \int_0^{2\pi} \int_0^{\pi/4} \int_0^{\sec\phi} \rho^2 \sin\phi\,\text{d}\rho\,\text{d}\phi\,\text{d}\theta

Integrate with respect to ρ\rho.

V=∫02π∫0π/4sec⁡3ϕsin⁡ϕ3 dϕ dθ=2π3∫0π/4tan⁡ϕsec⁡2ϕ dϕV = \int_0^{2\pi} \int_0^{\pi/4} \frac{\sec^3\phi \sin\phi}{3}\,\text{d}\phi\,\text{d}\theta = \frac{2\pi}{3} \int_0^{\pi/4} \tan\phi \sec^2\phi\,\text{d}\phi

Integrate with respect to ϕ\phi.

V=2π3[tan⁡2ϕ2]0π/4=2π3⋅12=π3V = \frac{2\pi}{3} \left[ \frac{\tan^2\phi}{2} \right]_0^{\pi/4} = \frac{2\pi}{3} \cdot \frac{1}{2} = \frac{\pi}{3}
Written by October 9, 2026 4 min read
Was this helpful?