Multiple Integration: Double, Triple Integrals & Beta-Gamma

Unit 4 — Multiple Integration (8 Hours)

Multiple integration extends the concept of integration to functions of two and three variables and provides essential tools for calculating areas, volumes, centroids, and moments of inertia in engineering and physics.

1. What Are Multiple Integrals?

A multiple integral is a generalization of the definite integral to functions of more than one variable. While single integrals compute area under a curve, double integrals compute volume under a surface, and triple integrals compute hypervolumes or accumulated quantities in 3D space.

Why They Matter

Calculate areas of complex planar regions
Compute volumes of solids and surfaces
Find centroids and moments of inertia in mechanics
Evaluate probabilities in multivariate statistics
Solve field problems in electromagnetism and fluid dynamics

2. Double Integral

2.1 Definition

Let $f(x, y)$ be defined on a closed region $R$ in the $xy$ -plane. The double integral of $f$ over $R$ is:

\iint_R f(x, y) \, dA = \lim_{\|P\| \to 0} \sum_{i=1}^{n} f(x_i^*, y_i^*) \, \Delta A_i

where $P$ is a partition of $R$ into $n$ subregions of area $\Delta A_i$ .

2.2 Evaluation Over Rectangular Regions

If $R = [a, b] \times [c, d]$ , then:

\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

Fubini's Theorem: If $f$ is continuous on $R$ , the order of integration can be interchanged.

2.3 Evaluation Over General Regions

Type I Region (bounded by vertical lines $x = a$ , $x = b$ and curves $y = g_1(x)$ , ):

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

Type II Region (bounded by horizontal lines $y = c$ , $y = d$ and curves $x = h_1(y)$ , ):

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

2.4 Applications to Area

The area of a region $R$ in the plane is:

A = \iint_R dA = \iint_R dx \, dy

If the region is bounded by $y = f(x)$ , $y = g(x)$ , $x = a$ , :

A = \int_a^b \int_{g(x)}^{f(x)} dy \, dx = \int_a^b [f(x) - g(x)] \, dx

3. Triple Integral

3.1 Definition

Let $f(x, y, z)$ be defined on a closed region $V$ in $\mathbb{R}^3$ . The triple integral is:

\iiint_V f(x, y, z) \, dV = \lim_{\|P\| \to 0} \sum_{i=1}^{n} f(x_i^*, y_i^*, z_i^*) \, \Delta V_i

3.2 Evaluation Over Rectangular Boxes

If $V = [a, b] \times [c, d] \times [p, q]$ :

\iiint_V f(x, y, z) \, dV = \int_a^b \int_c^d \int_p^q f(x, y, z) \, dz \, dy \, dx

The order of integration can be permuted (six possible orders) if $f$ is continuous.

3.3 Evaluation Over General 3D Regions

For a region bounded by surfaces $z = u_1(x, y)$ and $z = u_2(x, y)$ over a region in the -plane:

\iiint_V f(x, y, z) \, dV = \iint_R \left[ \int_{u_1(x, y)}^{u_2(x, y)} f(x, y, z) \, dz \right] dA

3.4 Applications to Volume

The volume of a solid region $V$ is:

V = \iiint_V dV = \iiint_V dx \, dy \, dz

If $V$ lies between $z = z_1(x, y)$ and $z = z_2(x, y)$ over region :

V = \iint_R [z_2(x, y) - z_1(x, y)] \, dA

4. Change of Order of Integration

4.1 Understanding Limits of Integration

When evaluating $\int_a^b \int_{g_1(x)}^{g_2(x)} f(x, y) \, dy \, dx$ , the inner limits may depend on the outer variable. Changing the order requires redescribing the region.

4.2 Reversing the Order

Procedure:

Sketch the region $R$ using the given limits
Identify the bounding curves
Express $x$ as a function of $y$ (or vice versa)
Rewrite the integral with new limits

Example Transformation:

Given:

\int_a^b \int_{y_1(x)}^{y_2(x)} f(x, y) \, dy \, dx

After changing order:

\int_c^d \int_{x_1(y)}^{x_2(y)} f(x, y) \, dx \, dy

4.3 Evaluation of Difficult Multiple Integrals

Some integrals are impossible or very difficult in one order but straightforward in the other. Always consider changing the order when:

The inner integral has no elementary antiderivative
The integrand involves functions like $e^{x^2}$ , $\frac{\sin x}{x}$ ,

5. Change of Variables

5.1 Transformation of Variables

When changing from $(x, y)$ to $(u, v)$ where $x = x(u, v)$ and $y$ :

\iint_R f(x, y) \, dx \, dy = \iint_{R'} f(x(u, v), y(u, v)) \left| \frac{\partial(x, y)}{\partial(u, v)} \right| \, du \, dv

5.2 Jacobian of Transformation

The Jacobian determinant for 2D:

J = \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}

For 3D transformations $(x, y, z) \to (u, v, w)$ :

J = \frac{\partial (x, y, z)}{\partial (u, v, w)} = ∣ \begin{array}{ccc} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} & \frac{\partial x}{\partial w} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} & \frac{\partial y}{\partial w} \\ \frac{\partial z}{\partial u} & \frac{\partial z}{\partial v} \end{array}

5.3 Common Coordinate Transformations

Polar Coordinates: $x = r\cos\theta$ , $y = r\sin\theta$

\frac{\partial(x, y)}{\partial(r, \theta)} = r \quad \Rightarrow \quad dx \, dy = r \, dr \, d\theta

Cylindrical Coordinates: $x = r\cos\theta$ , $y = r\sin\theta$ , $z = z$

\frac{\partial(x, y, z)}{\partial(r, \theta, z)} = r \quad \Rightarrow \quad dx \, dy \, dz = r \, dr \, d\theta \, dz

Spherical Coordinates: $x = \rho\sin\phi\cos\theta$ , $y = \rho\sin\phi\sin\theta$ , $z = \rho\cos\phi$

\frac{\partial(x, y, z)}{\partial(\rho, \phi, \theta)} = \rho^2 \sin\phi \quad \Rightarrow \quad dx \, dy \, dz = \rho^2 \sin\phi \, d\rho \, d\phi \, d\theta

5.4 Applications

Simplifying circular or spherical domains
Evaluating Gaussian integrals: $\int_{-\infty}^{\infty} e^{-x^2} dx = \sqrt{\pi}$

6. Beta and Gamma Functions

6.1 Definition of Beta Function

The Beta function is defined as:

B(m, n) = \int_0^1 x^{m-1} (1-x)^{n-1} \, dx, \quad m > 0, \, n > 0

Alternative Forms:

B(m, n) = 2 \int_0^{\frac{\pi}{2}} \sin^{2m-1}\theta \cos^{2n-1}\theta \, d\theta

B(m, n) = \int_0^{\infty} \frac{x^{m-1}}{(1+x)^{m+n}} \, dx

6.2 Definition of Gamma Function

The Gamma function is defined as:

\Gamma(n) = \int_0^{\infty} e^{-x} x^{n-1} \, dx, \quad n > 0

Important Values:

$n$	$\Gamma(n)$
$1$	$1$
$\frac{1}{2}$

6.3 Relationship Between Beta and Gamma Functions

B(m, n) = \frac{\Gamma(m)\Gamma(n)}{\Gamma(m+n)}

6.4 Important Properties

Gamma Function Properties:

$\Gamma(n+1) = n\Gamma(n)$ (Recursion relation)
$\Gamma(n+1) = n!$ for positive integers

Beta Function Properties:

$B(m, n) = B(n, m)$ (Symmetry)
$B(m, n) = \frac{(m-1)!(n-1)!}{(m+n-1)!}$ for positive integers

6.5 Wallis' Integrals

\int_0^{\frac{\pi}{2}} \sin^n x \, dx = \int_0^{\frac{\pi}{2}} \cos^n x \, dx = \begin{cases} \frac{(n-1)!!}{n!!} \cdot \frac{\pi}{2} & \text{if } n \text{ is even} \\ \frac{(n-1)!!}{n!!} & \text{if } n \text{ is odd} \end{cases}

where $n!!$ is the double factorial.

Using Gamma functions:

\int_0^{\frac{\pi}{2}} \sin^p x \cos^q x \, dx = \frac{1}{2} B\left(\frac{p+1}{2}, \frac{q+1}{2}\right) = \frac{\Gamma\left(\frac{p+1}{2}\right)\Gamma\left(\frac{q+1}{2}\right)}{2\Gamma\left(\frac{p+q+2}{2}\right)}

7. Dirichlet's Integral

7.1 Definition and Evaluation

Dirichlet's Integral generalizes the Beta function to multiple variables. For positive constants $l, m, n$ :

\iiint_V x^{l-1} y^{m-1} z^{n-1} \, dx \, dy \, dz = \frac{\Gamma(l)\Gamma(m)\Gamma(n)}{\Gamma(l+m+n+1)}

where $V$ is the region bounded by $x \geq 0$ , $y \geq 0$ , $z \geq 0$ and .

General Form:

\idotsint_V x_1^{l_1-1} x_2^{l_2-1} \cdots x_n^{l_n-1} \, dx_1 \cdots dx_n = \frac{\Gamma(l_1)\Gamma(l_2)\cdots\Gamma(l_n)}{\Gamma(l_1 + l_2 + \cdots + l_n + 1)}

over the region $x_i \geq 0$ and $\sum x_i \leq 1$ .

7.2 Applications to Area and Volume

Area of a 2D Simplex (Triangle): Setting $n = 2$ with $x + y \leq 1$ :

\iint_R dx \, dy = \frac{\Gamma(1)\Gamma(1)}{\Gamma(3)} = \frac{1}{2}

Volume of a 3D Simplex (Tetrahedron): Setting $n = 3$ with $x + y + z \leq 1$ :

\iiint_V dx \, dy \, dz = \frac{\Gamma(1)\Gamma(1)\Gamma(1)}{\Gamma(4)} = \frac{1}{6}

7.3 Generalization with Linear Constraint

For the region $x + y + z \leq h$ with $x, y, z \geq 0$ :

\iiint_V x^{l-1} y^{m-1} z^{n-1} \, dx \, dy \, dz = \frac{\Gamma(l)\Gamma(m)\Gamma(n)}{\Gamma(l+m+n+1)} \cdot h^{l+m+n}

8. Liouville's Extensions of Dirichlet's Integral

8.1 Extended Form

Liouville extended Dirichlet's integral to more general regions. If $f$ is continuous and $x_i \geq 0$ with $\sum x_i \leq 1$ :

\idotsint_V f(x_1 + x_2 + \cdots + x_n) \, x_1^{l_1-1} \cdots x_n^{l_n-1} \, dx_1 \cdots dx_n = \frac{\Gamma(l_1)\cdots\Gamma(l_n)}{\Gamma(l_1 + \cdots + l_n)} \int_0^1 f(u) \, u^{l_1 + \cdots + l_n - 1} \, du

8.2 First Extension

For $x, y, z \geq 0$ and $x + y + z \leq 1$ :

\iiint_V f(x + y + z) \, x^{l-1} y^{m-1} z^{n-1} \, dx \, dy \, dz = \frac{\Gamma(l)\Gamma(m)\Gamma(n)}{\Gamma(l+m+n)} \int_0^1 f(u) u^{l+m+n-1} \, du

8.3 Second Extension (Over Infinite Domain)

For $x, y, z \geq 0$ :

\iiint_V f(x + y + z) \, e^{-(x+y+z)} x^{l-1} y^{m-1} z^{n-1} \, dx \, dy \, dz = \frac{\Gamma(l)\Gamma(m)\Gamma(n)}{\Gamma(l+m+n)} \int_0^{\infty} f(u) e^{-u} u^{l+m+n-1} \, du

8.4 Applications in Multiple Integration

Reducing multiple integrals to single integrals
Evaluating integrals of the form $\idotsint f(\sum x_i) \, dV$
Computing statistical moments over simplex regions
Simplifying integrals in probability theory (Dirichlet distributions)

9. Engineering Applications

Domain	Application
Mechanical Engineering	Moments of inertia: $I = \iiint_V r^2 \rho \, dV$
Civil Engineering	Centroid and center of mass calculations

10. Common Mistakes to Avoid

Mistake	Correction
Forgetting the Jacobian when changing variables	Always multiply by $\\|J\\|$
Using wrong limits after changing integration order	Always sketch the region first
Assuming Fubini's theorem applies to discontinuous functions	Verify continuity before interchanging
Confusing Beta and Gamma arguments	$B(m,n)$ integrates $x^{m-1}(1-x)^{n-1}$ ; integrates

11. Best Practices & Architecture

Always sketch the region before setting up limits of integration
Choose coordinates wisely — polar for circles, cylindrical for cylinders, spherical for spheres
Verify Jacobian sign — use absolute value for geometric applications
Exploit symmetry — even/odd properties can reduce computation by half
Use Gamma/Beta functions for integrals involving powers of sine, cosine, or exponential decay
Apply Dirichlet/Liouville when integrating over simplex regions with power-law integrands

12. Practice Problems

Problem 1: Evaluate $\displaystyle \int_0^1 \int_0^x (x + y) \, dy \, dx$ and verify by changing the order of integration.

Problem 2: Find the area enclosed by the parabola $y = x^2$ and the line $y = 2x$ using a double integral.

Problem 3: Evaluate $\displaystyle \iint_R e^{x^2 + y^2} \, dA$ where is the unit disk using polar coordinates.

Problem 4: Find the volume of the solid bounded by the paraboloid $z = x^2 + y^2$ and the plane $z = 4$ .

Problem 5: Prove that $\displaystyle B\left(\frac{1}{2}, \frac{1}{2}\right) = \pi$ using the relationship with Gamma functions.

Problem 6: Evaluate $\displaystyle \int_0^{\infty} x^4 e^{-2x} \, dx$ using Gamma functions.

Problem 7: Evaluate the Dirichlet integral $\displaystyle \iiint_V x y z \, dx \, dy \, dz$ over the region $x + y + z \leq 1$ in the first octant.

Problem 8: Use Liouville's extension to evaluate $\displaystyle \iiint_V (x + y + z) \, dx \, dy \, dz$ over , .

13. Mini Project

Title: Volume and Center of Mass of a Complex 3D Component

Objective: Use triple integrals, coordinate transformations, and special functions to analyze an engineering component.

Given:

A machine component occupies the region bounded by the cylinder $x^2 + y^2 = 4$ , the plane $z = 0$ , and the paraboloid

Tasks:

Set up and evaluate the triple integral for the volume using cylindrical coordinates
Compute the total mass $M = \iiint_V \rho \, dV$
Find the center of mass coordinates $(\bar{x}, \bar{y}, \bar{z})$ :

Deliverable: A computational notebook (Python/MATLAB) with symbolic integration, 3D visualization of the region, and validation of results using alternative coordinate systems.

Quick Reference Formula Sheet

Concept	Formula
Double integral (Type I)	$\int_a^b \int_{g_1(x)}^{g_2(x)} f(x,y) \, dy \, dx$

End of Unit 4 — Multiple Integration