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).

Double Integral Formula

The basic formula for evaluating a double integral over a region R.

\iint_R f(x, y) \, dA

Variables

Symbol	Description	Unit
$R$	Region of integration in the xy-plane	-
$f (x, y)$	Integrand function (height of the surface)	-
$d A$	Differential area element	-

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

If $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.

Iterated Integrals

Evaluation of a double integral by integrating one variable at a time.

\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

Variables

Symbol	Description	Unit
$a, b$	Limits of integration for x	-
$c, d$	Limits of integration for y	-

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.

Type I Region Integral

Iterated double integral integrating with respect to y first.

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

Variables

Symbol	Description	Unit
$g_{1} (x)$	Lower bounding curve	-
$g_{2} (x)$	Upper bounding curve	-
$a, b$	Constant limits for x	-

Type II Region Integral

Iterated double integral integrating with respect to x first.

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

Variables

Symbol	Description	Unit
$h_{1} (y)$	Left bounding curve	-
$h_{2} (y)$	Right bounding curve	-
$c, d$	Constant limits for y	-

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$ !

Interact with the simulation below to explore double integrals in polar coordinates, visualizing area elements and polar wedges.

Polar Coordinates: Double Integrals Visualizer

Understand why the polar area element is $dA = r \, dr \, d\theta$ . Sweeping further from the origin enlarges the wedge area, proving that the extra $r$ factor represents coordinate scaling.

Select Integrand Function $f(r, \theta)$

Radius Limit (

R_{max}

): 3 m

Angle Limit (

\theta_{max}

): 180°

Wedge Element

dA

radius

r

: 2.0 m

Double Integral Formulation

\iint_R f(r, \theta) \, dA = \int_0^{\theta_{max}} \int_0^{R_{max}} r \, dr \, d\theta

Wedge Area dA:0.1200 m²

Integral Value:14.1372 m²

Polar Sweeping Grid

The Polar Area Element

In Cartesian coordinates, the differential area is a constant rectangle: $dA = dx \\, dy$ .

In Polar coordinates, the area element is a curved wedge of radial thickness $dr$ and arc length $r \\, d\\theta$ .

Therefore, the area of the wedge is:

dA = (r \\, d\\theta) \\cdot dr = r \\, dr \\, d\\theta

Drag the slider and watch how $dA$ (the red wedge) grows physically larger as you slide $r$ outwards. Closer to the origin, wedges are tightly squeezed; further out, they expand. The extra $r$ factor mathematically offsets this geometrical widening!

Polar Double Integral

The transformation introduces the scaling factor $r$ to the differential area element.

Polar Double Integral Formula

Transforming a double integral to polar coordinates.

\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

Variables

Symbol	Description	Unit
$r$	Radial distance (Jacobian scaling factor)	-
$\theta$	Angular coordinate	-
$\alpha, \beta$	Angular limits of integration	-

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.

Area of a Planar Region

Calculating area using a double integral.

A = \iint_R 1 \, dA

Variables

Symbol	Description	Unit
$A$	Area of region R	-

Mass of a Lamina

Calculating mass of a 2D object with variable density.

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

Variables

Symbol	Description	Unit
$m$	Total mass	-
$\rho(x,y)$	Density function	-

Center of Mass (X)

The x-coordinate of the center of mass.

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

Variables

Symbol	Description	Unit
$\bar{x}$	x-coordinate of the center of mass	-
$m$	Total mass	-
$\rho(x,y)$	Density function	-

Center of Mass (Y)

The y-coordinate of the center of mass.

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

Variables

Symbol	Description	Unit
$\bar{y}$	y-coordinate of the center of mass	-
$m$	Total mass	-
$\rho(x,y)$	Density function	-

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).

Triple Integral Formula

General notation for a triple integral over a 3D solid E.

\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

Variables

Symbol	Description	Unit
$E$	3D solid region of integration	-
$d V$	Differential volume element (e.g., dx dy dz)	-

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.

Total Mass

Calculating total mass in a 3D solid.

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

Variables

Symbol	Description	Unit
$m$	Total mass	-
$\rho(x,y,z)$	Density function in 3D	-

Center of Mass (X, 3D)

The x-coordinate of the center of mass in 3D.

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

Variables

Symbol	Description	Unit
$\bar{x}$	x-coordinate of the center of mass	-
$m$	Total mass	-

Moment of Inertia (Z-axis)

The moment of inertia around the z-axis.

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

Variables

Symbol	Description	Unit
$I_{z}$	Moment of inertia around the z-axis	-
$x^{2} + y^{2}$	Square of distance from z-axis	-

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.

Surface Area Formula

Parametric surface area formula using cross products.

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

Variables

Symbol	Description	Unit
$A (S)$	Surface area	-
$\mathbf{r}_u, \mathbf{r}_v$	Partial derivative vectors	-
$\|\mathbf{r}_u \times \mathbf{r}_v\|$	Magnitude of the normal vector	-

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.

Surface Area Function Formula

Simplified surface area formula for explicitly defined functions.

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

Variables

Symbol	Description	Unit
$S$	Surface area	-
$\frac{\partial f}{\partial x}, \frac{\partial f}{\partial y}$	Partial derivatives of f	-

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:

Jacobian Determinant

Calculating the scaling factor for a change of variables in double integrals.

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}

Variables

Symbol	Description	Unit
$J (u, v)$	Jacobian determinant	-
$\frac{\partial(x, y)}{\partial(u, v)}$	Notation for the Jacobian	-

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