Multivariable Integration- Techniques and Applications
What Multivariable Integration Actually Is
Single-variable integration was simple: find the area under a curve. Multivariable integration does the same thing but in higher dimensions. You're summing up tiny pieces across a region, not just along a line.
The math gets messier. The concepts stay grounded. You still build up quantities by adding infinitely small pieces—you're just doing it over surfaces and volumes instead of intervals.
This isn't advanced math for its own sake. Engineers, physicists, and anyone modeling the real world need these tools constantly.
Double Integrals: Volume Under a Surface
The double integral is your starting point. You're integrating a function f(x, y) over a region R in the xy-plane.
Think of it as stacking infinitely thin volume elements: dA = dx dy. The result is the volume trapped between the surface and the plane below it.
Order of Integration
You can integrate x first, then y. Or flip it. The result stays the same if you're careful with your limits.
For rectangular regions, the math is straightforward. For irregular regions, you need to set up your limits properly—which is where most people mess up.
Sketch the region first. Always. Trying to set limits without visualization is guessing.
Iterated Integrals
You write double integrals as iterated single-variable integrals:
∬R f(x,y) dA = ∫ab ∫cd f(x,y) dy dx
The inner integral gets evaluated first, treating the outer variable as constant. Work from the inside out.
Triple Integrals: Adding a Third Dimension
Triple integrals extend the same logic to volumes. You're integrating f(x, y, z) over a three-dimensional region.
The volume element becomes dV = dx dy dz. The result could represent mass, charge, probability—whatever quantity you're accumulating.
Setting up limits in 3D is harder. You need to understand the shape of your region and project it onto coordinate planes.
When the Order Matters
Sometimes one order makes the problem solvable. Sometimes another order saves you hours of grinding through ugly integrals.
If one approach looks brutal, switch. Fubini's Theorem says the answer is the same regardless—just expect different levels of suffering getting there.
The Jacobian: Change of Variables
Sometimes the region or integrand screams for a different coordinate system. Spherical coordinates for spheres. Polar for circles. Cylindrical for pipes.
When you switch variables, you need a correction factor. That's the Jacobian determinant.
For a transformation (u, v) → (x, y):
J = |∂(x,y)/∂(u,v)| = xuyv - xvyu
The Jacobian tells you how much areas or volumes stretch when you deform your coordinate system. You multiply your integrand by |J| when making the switch.
Polar Coordinates: The Common Case
Going from (r, θ) to (x, y) gives J = r. That's why you get that extra r in polar integrals:
∬R f(x,y) dA = ∫∫ f(r cos θ, r sin θ) · r dr dθ
Spherical coordinates have their own Jacobian: ρ2 sin φ.
Line Integrals: Integration Along Paths
Line integrals break from the pattern. Instead of integrating over a region, you integrate along a curve.
Line Integral of a Scalar Field
For a scalar function f along curve C:
∫C f ds
Here ds is the differential of arc length. You parameterize the curve, find ds = |r'(t)| dt, and plug through.
Line Integral of a Vector Field
This is where physics kicks in. The line integral of a vector field F along curve C:
∫C F · dr
This gives work done by a force field along a path, or circulation, or flux depending on how you set it up.
Parameterize your curve, compute dr, take the dot product, integrate. That's the whole procedure.
Surface Integrals: Integrating Over Surfaces
Surface integrals extend integration to two-dimensional surfaces in 3D space.
Scalar Surface Integrals
For a scalar function over surface S:
∬S f dS
The element dS is surface area. You parameterize the surface, compute the magnitude of the cross product of partial derivatives, and integrate.
Vector Surface Integrals (Flux)
For a vector field flowing through a surface:
∬S F · dS = ∬S F · n dS
This is flux. How much of the field passes through the surface. Critical in electromagnetism and fluid dynamics.
Green's, Stokes', and the Divergence Theorem
These are the big theorems that tie everything together. They let you convert between integrals over regions and integrals over boundaries.
Green's Theorem
In the plane, for a closed curve C bounding region R:
∮C F · dr = ∬R (∂Q/∂x - ∂P/∂y) dA
Line integral around a closed curve equals double integral over the region inside. Useful for computing circulation or converting ugly line integrals into manageable double integrals.
Stokes' Theorem
The 3D version. Surface integral of curl F over S equals line integral of F around the boundary curve:
∬S (∇ × F) · dS = ∮∂S F · dr
Divergence Theorem
Flux through a closed surface equals triple integral of divergence in the volume inside:
∮S F · dS = ∭V (∇ · F) dV
This is Gauss's Law in disguise. Electric flux through a closed surface relates to enclosed charge.
Practical Applications
You won't find this in textbooks, but here's where multivariable integration actually shows up:
- Calculating mass from density functions over irregular shapes
- Moments of inertia for engineering design—how resistant a shape is to rotation
- Electric and gravitational potential from distributed charges or masses
- Fluid flow through curved pipes or around obstacles
- Heat transfer in multidimensional systems
- Probability in multiple dimensions—joint distributions, expectations
- Computer graphics for rendering curved surfaces and calculating lighting
Common Mistakes to Avoid
- Forgetting the Jacobian when changing variables—guaranteed wrong answer
- Setting wrong limits on iterated integrals—sketch the region first
- Confusing dS and dV—different things, different uses
- Losing track of which coordinate system you're in
- Forgetting that surface orientation matters for flux calculations
- Treating line integrals and surface integrals the same way—they're not
Getting Started: A Worked Example
Evaluate ∬R (x² + y²) dA where R is the region in the first quadrant bounded by x² + y² = 4 and x² + y² = 9.
Step 1: Recognize this screams polar coordinates. The region is an annular sector in the first quadrant.
Step 2: Set up in polar.
- r goes from 2 to 3
- θ goes from 0 to π/2
- x² + y² = r², so the integrand becomes r²
- The Jacobian is r
Step 3: Write the integral.
∫0π/2 ∫23 r² · r dr dθ = ∫0π/2 ∫23 r³ dr dθ
Step 4: Evaluate.
∫23 r³ dr = [r⁴/4]23 = (81/4) - (16/4) = 65/4
∫0π/2 (65/4) dθ = (65/4)(π/2) = 65π/8
The answer is 65π/8.
Quick Reference: Coordinate Systems and Jacobians
| System | Variables | Volume Element | Best For |
|---|---|---|---|
| Cartesian | x, y, z | dx dy dz | Rectangular everything |
| Cylindrical | r, θ, z | r dr dθ dz | Pipes, symmetry about axis |
| Spherical | ρ, φ, θ | ρ² sin φ dρ dφ dθ | Spheres, point sources |
When to Use What
Choosing the right coordinate system saves time. Here's the real-world breakdown:
- Circular symmetry in 2D → polar
- Cylindrical symmetry in 3D → cylindrical
- Spherical symmetry in 3D → spherical
- No obvious symmetry → Cartesian and suffer through it
Sometimes you'll need to split a region into pieces and integrate separately. That's fine. Reality doesn't give you clean boundaries.
The Bottom Line
Multivariable integration is single-variable integration with extra dimensions. The core idea—summing small pieces to get a whole—doesn't change. What changes is the geometry you have to handle.
Master the coordinate systems. Learn when to use the big theorems. Practice setting up integrals correctly. The evaluation is often just arithmetic—the setup is where the thinking happens.