Multiple Integrals - Theory & Concepts

Multiple Integrals

Just as a definite integral of a single variable function calculates the area under a curve, multiple integrals extend this concept to calculate the volume under a surface (double integral) or higher-dimensional hypervolumes (triple integral). They are essential for computing physical properties like mass, center of mass, and moments of inertia for multi-dimensional objects with variable density.

Double Integrals

A function of two variables,

f(x, y)

, generates a surface in 3D space. The double integral of this function over a two-dimensional region

R

in the xy-plane calculates the signed volume between the surface

z=f(x,y)

and the region

R

Double Integral Notation

$\iint$ indicates two consecutive integrations.
$R$ is the region of integration in the xy-plane.
$f(x, y)$ is the integrand (the height of the surface).
$dA$ is the differential area element ( $dx\,dy$ or $dy\,dx$ in rectangular coordinates).

$$
\\iint_R f(x, y) \\, dA
$$

Double Integral & Volume Visualizer

Function $z = f(x, y)$

Grid Resolution (n x n): 10 x 10

Increase resolution to see the Riemann sum converge to the exact volume. Note: High resolution may affect performance.

Calculated Volume (Approx.)

25.190

V \approx \sum_{i=1}^n \sum_{j=1}^n f(x_i^*, y_j^*) \Delta A

Drag to rotate • Scroll to zoom

Iterated Integrals and Fubini's Theorem

To evaluate a double integral, we usually express it as an iterated integral, performing the integrations one variable at a time while holding the other constant.

Fubini's Theorem

f(x,y)

is continuous on a rectangular region

R = [a,b] \times [c,d]

, then the double integral can be evaluated as either iterated integral.

$$
\\iint_R f(x, y) \\, dA = \\int_a^b \\int_c^d f(x, y) \\, dy \\, dx = \\int_c^d \\int_a^b f(x, y) \\, dx \\, dy
$$

Double Integrals over General Regions

When the region

R

is not a simple rectangle, the limits of the inner integral become functions of the outer variable.

Type I vs. Type II Regions

Type I Region (Vertical Slices): The region is bounded above and below by continuous functions $y = g_2(x)$ and $y = g_1(x)$ , and lies between $x=a$ and $x=b$ . Integrate with respect to $y$ first.
Type II Region (Horizontal Slices): The region is bounded left and right by continuous functions $x = h_1(y)$ and $x = h_2(y)$ , and lies between $y=c$ and $y=d$ . Integrate with respect to $x$ first.

$$
\\iint_R f(x, y) \\, dA = \\int_a^b \\left[ \\int_{g_1(x)}^{g_2(x)} f(x, y) \\, dy \\right] \\, dx
$$

$$
\\iint_R f(x, y) \\, dA = \\int_c^d \\left[ \\int_{h_1(y)}^{h_2(y)} f(x, y) \\, dx \\right] \\, dy
$$

Double Integrals in Polar Coordinates

If a region has circular symmetry (like a circle, annulus, or a sector), it is often much simpler to evaluate the double integral using polar coordinates

(r, \theta)

instead of rectangular coordinates

(x, y)

Changing to Polar Coordinates

STEP-BY-STEP

Substitute $x = r \cos \theta$ and $y = r \sin \theta$ into the integrand.
Determine the new limits of integration for $r$ and (\theta) based on the region.
Replace the area differential $dA = dx\,dy$ with $dA = r \, dr \, d\theta$ . Do not forget the extra $r$ !

Polar Double Integral

The transformation introduces the scaling factor

r

to the differential area element.

$$
\\iint_R f(x,y) \\, dA = \\int_{\\alpha}^{\\beta} \\int_{h_1(\\theta)}^{h_2(\\theta)} f(r\\cos\\theta, r\\sin\\theta) \\cdot r \\, dr \\, d\\theta
$$

Applications of Double Integrals

Double integrals have numerous physical and geometric applications beyond calculating volume.

Checklist

Area of a Planar Region: If the integrand $f(x,y) = 1$ , the volume calculated is numerically equal to the area of the base region $R$ .
Mass of a Lamina: If a flat object (lamina) has a variable density $\rho(x,y)$ , its total mass $m$ is the double integral of the density function over the object's region.
Center of Mass: The coordinates of the center of mass $(\bar{x}, \bar{y})$ are found using moments.

$$
A = \\iint_R 1 \\, dA
$$

$$
m = \\iint_R \\rho(x,y) \\, dA
$$

$$
\\bar{x} = \\frac{1}{m} \\iint_R x \\rho(x,y) \\, dA
$$

$$
\\bar{y} = \\frac{1}{m} \\iint_R y \\rho(x,y) \\, dA
$$

Triple Integrals

A triple integral extends integration to functions of three variables,

f(x, y, z)

, over a three-dimensional solid region

E

. If

f(x,y,z) = 1

, the triple integral calculates the volume of the solid

E

. If

f

represents density, the integral calculates the total mass of the solid.

Triple Integral in Rectangular Coordinates

The volume differential in rectangular coordinates is

dV = dx \, dy \, dz

(or any of the 6 possible orderings).

$$
\\iiint_E f(x,y,z) \\, dV = \\int_{x=a}^{x=b} \\int_{y=g_1(x)}^{y=g_2(x)} \\int_{z=u_1(x,y)}^{z=u_2(x,y)} f(x,y,z) \\, dz \\, dy \\, dx
$$

3D Mass and Moments of Inertia

Triple integrals allow us to analyze physical solids with varying densities

\rho(x,y,z)

Checklist

Total Mass ( $m$ ): The integral of the density function over the solid.
Center of Mass Coordinates: Calculated using the moments about the coordinate planes.
Moments of Inertia: To find the rotational resistance of a solid around an axis, integrate the density times the square of the distance from that axis.

$$
m = \\iiint_E \\rho(x,y,z) \\, dV
$$

$$
\\bar{x} = \\frac{1}{m} \\iiint_E x \\rho(x,y,z) \\, dV
$$

$$
I_z = \\iiint_E (x^2 + y^2) \\rho(x,y,z) \\, dV
$$

Cylindrical Coordinates

Cylindrical coordinates

(r, \theta, z)

are an extension of polar coordinates into 3D space. They are particularly useful for regions

E

that possess an axis of symmetry, such as cylinders or cones.

Cylindrical Coordinate Conversions

$x = r \cos \theta$
$y = r \sin \theta$
$z = z$
Volume differential: $dV = r \, dz \, dr \, d\theta$

Spherical Coordinates

Spherical coordinates

(\rho, \theta, \phi)

define a point in 3D space by its distance from the origin

(\rho)

, its angle in the xy-plane (

\theta

, same as polar), and its angle from the positive z-axis (

\phi

). They are ideal for integrating over spheres or portions of spheres.

Spherical Coordinate Conversions

$x = \rho \sin \phi \cos \theta$
$y = \rho \sin \phi \sin \theta$
$z = \rho \cos \phi$
Volume differential: $dV = \rho^2 \sin \phi \, d\rho \, d\phi \, d\theta$

Surface Area and Jacobians

Surface Area of a Parametric Surface

Just as we use a single integral to find the arc length of a curve, we can use a double integral to find the surface area of a 3D surface. While we previously discussed surfaces defined by

z = f(x, y)

, more complex surfaces in engineering (like airplane fuselages or terrain meshes) are defined parametrically using vectors:

\mathbf{r}(u, v) = x(u,v)\mathbf{i} + y(u,v)\mathbf{j} + z(u,v)\mathbf{k}

Parametric Surface Area Formula

The area of a surface defined parametrically over a region

D

in the

uv

-plane is found by integrating the magnitude of the cross product of the partial derivative vectors

\mathbf{r}_u

and

\mathbf{r}_v

. This cross product geometrically represents the normal vector, and its magnitude represents the area of a differential parallelogram on the surface.

$$
A(S) = \\iint_D |\\mathbf{r}_u \\times \\mathbf{r}_v| \\, dA
$$

Surface Area for z = f(x, y)

If the surface is defined simply as a function of

x

and

y

(i.e.,

z = f(x, y)

over a region

R

), the parametric formula simplifies to a specific form.

$$
S = \\iint_R \\sqrt{1 + \\left(\\frac{\\partial f}{\\partial x}\\right)^2 + \\left(\\frac{\\partial f}{\\partial y}\\right)^2} \\, dA
$$

Change of Variables: The Jacobian

When transforming a double integral from one coordinate system to another (e.g., from rectangular to polar, or to a custom generic

u,v

system), the differential area

dA

scales by a factor called the Jacobian.

The Jacobian Determinant

For a transformation given by

x = g(u, v)

and

y = h(u, v)

, the Jacobian is defined as the determinant of the partial derivatives:

$$
J(u, v) = \\frac{\\partial(x, y)}{\\partial(u, v)} = \\begin{vmatrix} \\frac{\\partial x}{\\partial u} & \\frac{\\partial x}{\\partial v} \\\\ \\frac{\\partial y}{\\partial u} & \\frac{\\partial y}{\\partial v} \\end{vmatrix} = \\frac{\\partial x}{\\partial u}\\frac{\\partial y}{\\partial v} - \\frac{\\partial x}{\\partial v}\\frac{\\partial y}{\\partial u}
$$

Why Polar Uses 'r'

The extra

r

in the polar area differential

dA = r \, dr \, d\theta

is actually the absolute value of the Jacobian for the transformation

x = r \cos \theta, y = r \sin \theta

Bridge to Vector Calculus

Multiple integrals form the foundation for evaluating complex fields and flows in engineering. The capstone theorems of Vector Calculus relate integration over boundaries (curves/surfaces) to integration over their interiors (areas/volumes).

Fundamental Theorems of Vector Calculus

Green's Theorem: Relates a line integral over a closed plane curve to a double integral over the bounded planar region. Useful for calculating macroscopic circulation and fluid flux.
Stokes' Theorem: The 3D extension of Green's Theorem. It relates a line integral around a boundary curve to a surface integral (which is evaluated using a double integral parametrization) over the enclosed 3D surface.
Divergence Theorem: Relates the surface integral measuring the total flow (flux) out of a closed surface to a triple integral of the field's divergence throughout the enclosed volume.

Key Takeaways

Double integrals calculate the volume under a 3D surface $z=f(x,y)$ over a 2D region $R$ by evaluating an iterated integral (integrating one variable at a time).
Fubini's Theorem guarantees that the order of integration ( $dx\,dy$ or $dy\,dx$ ) does not matter for continuous functions over rectangular regions, but switching order on general regions requires carefully redefining the limits.
Changing to polar coordinates simplifies integrals over circular regions. Always remember to replace $dx\,dy$ with the polar area differential $r\,dr\,d\theta$ .
Triple integrals compute hypervolumes or total mass, center of mass, and moments of inertia for 3D solids. They can be evaluated in rectangular ( $dV = dx\,dy\,dz$ ), cylindrical ( $dV = r\,dz\,dr\,d\theta$ ), or spherical ( $dV = \rho^2 \sin \phi \, d\rho \, d\phi \, d\theta$ ) systems.
Parametric Surface Area: Use the cross product of partial derivative vectors ( $|\mathbf{r}_u \times \mathbf{r}_v|$ ) to find the surface area of complex geometries.
Jacobians provide the critical scaling factor for the differential integration area or volume during generalized coordinate transformations.
Multiple integrals are conceptually bridged to vector fields via Green's, Stokes', and the Divergence Theorems for calculating real-world electromagnetic and fluid behaviors.

PreviousApplications of Integration - Examples & Applications

Quiz Me

NextMultiple Integrals - Examples & Applications