Boo-AI — Master Artificial Intelligence by Building from Scratch

Learning Objectives

By the end of this section, you will be able to:

Visualize how rotating a curve around an axis creates a three-dimensional solid of revolution
Derive the disk method formula from first principles using Riemann sums
Apply the disk method formula $V = \\pi \\int_a^b [f(x)]^2 \\, dx$ to compute volumes
Set up disk method integrals for rotation around both the x-axis and y-axis
Connect the disk method to real-world applications in engineering and physics
Implement numerical volume calculations in Python

The Big Picture: From 2D Curves to 3D Solids

"The disk method transforms a one-dimensional curve into a three-dimensional object by the simple act of rotation — and integration captures its volume."

In the previous section, we learned to find the area between curves. Now we take a dramatic leap into three dimensions: what happens when we rotate a region around an axis?

Imagine taking a piece of paper with a curve drawn on it and spinning it around a fixed axis like a pottery wheel. The curve sweeps out a three-dimensional shape called a solid of revolution. This process creates many familiar shapes:

🎯 Simple Examples

Cone: Rotate a line $y = mx$ around the x-axis
Sphere: Rotate a semicircle around its diameter
Cylinder: Rotate a horizontal line around the x-axis
Paraboloid: Rotate $y = \\sqrt{x}$ around the x-axis

🏭 Engineering Applications

Vases and bowls: Pottery shapes
Rocket nozzles: Aerodynamic design
Wine glasses: Calculating liquid volume
Storage tanks: Industrial containers

The Central Question

Given a function $y = f(x)$ on $[a, b]$ , how do we calculate the volume of the solid formed when this curve is rotated around an axis?

Historical Context: From Archimedes to Cavalieri

The problem of finding volumes of curved solids has ancient origins. Around 250 BCE, Archimedes of Syracuse discovered that the volume of a sphere equals two-thirds the volume of its circumscribing cylinder — a result he was so proud of that he requested it be inscribed on his tombstone.

Archimedes used an ingenious method of "balancing" infinitesimally thin slices on a lever — an early form of integral calculus 1,900 years before Newton and Leibniz!

Cavalieri's Principle (1635)

Italian mathematician Bonaventura Cavalieri formalized the intuition behind slicing methods. His principle states:

Cavalieri's Principle: If two solids have the same cross-sectional area at every height, they have the same volume.

This principle justifies the disk method: we compute volume by summing the areas of all cross-sectional slices.

Geometric Intuition: Slicing the Solid

The key insight of the disk method is that when we slice a solid of revolution perpendicular to the axis of rotation, each slice is a circular disk.

Step-by-Step Visualization

Start with a curve: Consider $y = f(x)$ on $[a, b]$
Rotate around the x-axis: The curve sweeps out a 3D surface
Slice perpendicular to x-axis: Each slice at position $x$ is a circular disk
Measure the disk: Radius $r = f(x)$ , Area $A = \\pi r^2 = \\pi [f(x)]^2$
Sum infinitely thin disks: Integrate from $a$ to $b$

Why Circles?

When you rotate a point at height $y = f(x)$ around the x-axis, it traces a circle of radius $f(x)$ . All points at the same x-value rotate to form a circular disk!

Interactive Cross-Section Explorer

Use the interactive tool below to explore how cross-sectional disks change as you move along the x-axis. Notice how the radius of each disk equals $f(x)$ at that position.

Cross-Section Explorer

Move the slider to see how the cross-sectional area changes at different x positions. Each cross-section is a circle (disk) with radius $r = f(x)$ .

Slice Position: x = 1.00

Radius

r = f(1.00) = 1.0000

Cross-Sectional Area

A = π(1.000)² = 3.1416

Formula

A(x) = \pi [\sqrt{x}]^2

Key Insight: When we revolve y = f(x) around the x-axis, each vertical slice creates a circular cross-section. The radius of this circle equals f(x), so the area is $A(x) = \pi [f(x)]^2$ . Integrating these infinitesimally thin disks from x = a to x = b gives the total volume.

Deriving the Disk Method

Let's derive the disk method formula rigorously using our knowledge of Riemann sums and definite integrals.

Setup

Consider a continuous function $y = f(x) \\geq 0$ on $[a, b]$ . We rotate the region under this curve around the x-axis to create a solid.

Step 1: Partition the Interval

Divide $[a, b]$ into $n$ subintervals of equal width:

\\Delta x = \\frac{b - a}{n}

Step 2: Approximate Each Slice as a Disk

At the $i$ -th subinterval, choose a sample point $x_i^*$ . The cross-section at this point is approximately a disk with:

Radius: $r_i = f(x_i^*)$
Thickness: $\\Delta x$
Volume: $\\Delta V_i = \\pi r_i^2 \\cdot \\Delta x = \\pi [f(x_i^*)]^2 \\Delta x$

Step 3: Sum All Disk Volumes

The total volume is approximately:

V \\approx \\sum_{i=1}^{n} \\pi [f(x_i^*)]^2 \\Delta x

Step 4: Take the Limit

As $n \\to \\infty$ (more, thinner disks), the approximation becomes exact:

V = \\lim_{n \\to \\infty} \\sum_{i=1}^{n} \\pi [f(x_i^*)]^2 \\Delta x = \\pi \\int_a^b [f(x)]^2 \\, dx

The Disk Method Formula

Disk Method (Rotation Around x-axis)

V = \\pi \\int_a^b [f(x)]^2 \\, dx

where

f(x) \\geq 0

[a, b]

Understanding Each Component

Symbol	Meaning	Interpretation
π	Constant pi	Appears because cross-sections are circles (A = πr²)
∫ᵃᵇ	Definite integral	Sums up infinitely many infinitesimally thin disks
[f(x)]²	Squared function	The radius of each disk is f(x), so area = πr² = π[f(x)]²
dx	Infinitesimal thickness	Each disk has thickness dx

Don't Forget to Square!

A very common mistake is to forget to square the function. The integrand is $[f(x)]^2$ , not $f(x)$ . This is because the area of a circle is $\\pi r^2$ , not $\\pi r$ .

Interactive 3D Visualizer

Explore the disk method in 3D! Rotate the view by dragging, adjust the number of disks, and watch as the approximation converges to the exact volume.

3D Disk Method Visualizer

Number of Disks: 6

Interval

[0.00, 4.00]

Δx

0.6667

Disk Approximation

V ≈ 25.1327

Sum of 6 disks

Exact Volume

V = 25.1327

\pi \int_{0}^{4} [\sqrt{x}]^2 \, dx

Error

0.00%

|Approx - Exact| / Exact

Disk Method Formula:

V = \pi \int_a^b [f(x)]^2 \, dx \approx \sum_{i=1}^{n} \pi [f(x_i^*)]^2 \Delta x

Each disk has volume $\pi r^2 \cdot \text{thickness}$ where $r = f(x_i^*)$ is the radius and thickness is $\Delta x$ .

Worked Examples

Example 1: Cone Volume

Find the volume of the cone formed by rotating $y = \\frac{R}{h}x$ around the x-axis from $x = 0$ to $x = h$ .

Solution:

Using the disk method with $f(x) = \\frac{R}{h}x$ :

V = \\pi \\int_0^h \\left[\\frac{R}{h}x\\right]^2 dx

= \\pi \\cdot \\frac{R^2}{h^2} \\int_0^h x^2 \\, dx

= \\pi \\cdot \\frac{R^2}{h^2} \\cdot \\left[\\frac{x^3}{3}\\right]_0^h

= \\pi \\cdot \\frac{R^2}{h^2} \\cdot \\frac{h^3}{3}

= \\frac{1}{3}\\pi R^2 h

Result: $V = \\frac{1}{3}\\pi R^2 h$ ✓

This confirms the familiar cone volume formula: one-third of the cylinder volume!

Example 2: Paraboloid

Find the volume of the solid formed by rotating $y = \\sqrt{x}$ around the x-axis from $x = 0$ to $x = 4$ .

Solution:

V = \\pi \\int_0^4 [\\sqrt{x}]^2 \\, dx

= \\pi \\int_0^4 x \\, dx

= \\pi \\left[\\frac{x^2}{2}\\right]_0^4

= \\pi \\left(\\frac{16}{2} - 0\\right) = 8\\pi

Result: $V = 8\\pi \\approx 25.13$

Example 3: Sphere

Derive the sphere volume formula by rotating the semicircle $y = \\sqrt{R^2 - x^2}$ around the x-axis.

Solution:

V = \\pi \\int_{-R}^{R} [\\sqrt{R^2 - x^2}]^2 \\, dx

= \\pi \\int_{-R}^{R} (R^2 - x^2) \\, dx

= \\pi \\left[R^2 x - \\frac{x^3}{3}\\right]_{-R}^{R}

= \\pi \\left[\\left(R^3 - \\frac{R^3}{3}\\right) - \\left(-R^3 + \\frac{R^3}{3}\\right)\\right]

= \\pi \\left(\\frac{2R^3}{3} + \\frac{2R^3}{3}\\right) = \\frac{4}{3}\\pi R^3

Result: $V = \\frac{4}{3}\\pi R^3$ ✓

This confirms the famous sphere volume formula that Archimedes discovered!

Rotating Around the y-axis

When we rotate a curve around the y-axis instead of the x-axis, we can still use the disk method — but we need to set up the integral differently.

Key Changes

Express $x$ as a function of $y$ : $x = g(y)$
The radius of each disk is now $r = g(y)$ (horizontal distance)
Integrate with respect to $y$ from $c$ to $d$

Disk Method (Rotation Around y-axis)

V = \\pi \\int_c^d [g(y)]^2 \\, dy

where

x = g(y) \\geq 0

[c, d]

Example: Paraboloid (y-axis rotation)

Find the volume of the solid formed by rotating $x = \\sqrt{y}$ around the y-axis from $y = 0$ to $y = 4$ .

Solution:

V = \\pi \\int_0^4 [\\sqrt{y}]^2 \\, dy

= \\pi \\int_0^4 y \\, dy

= \\pi \\left[\\frac{y^2}{2}\\right]_0^4 = 8\\pi

Result: $V = 8\\pi$

The Washer Method (Preview)

What if there's a hole in the middle of your solid? This happens when you rotate a region between two curves, creating a hollow center.

In this case, each cross-section is not a disk but a washer (or annulus) — a disk with a circular hole cut out.

Washer Method Formula (Preview):

V = \\pi \\int_a^b \\left([R(x)]^2 - [r(x)]^2\\right) dx

where $R(x)$ = outer radius and $r(x)$ = inner radius

We will explore the washer method in detail in the next section!

Real-World Applications

🚀 Aerospace Engineering

Rocket nozzles are designed as solids of revolution. Engineers use the disk method to calculate nozzle volumes and optimize fuel efficiency.

🏭 Manufacturing

Lathes create solids of revolution by spinning material. Calculating volumes helps estimate material costs and weight.

🍷 Food & Beverage

Wine glass and bottle designs are solids of revolution. Volume calculations ensure accurate labeling and consistent serving sizes.

🏗️ Civil Engineering

Domes, silos, and water towers are often designed as surfaces of revolution. Volume calculations are essential for capacity planning.

Machine Learning Connections

The disk method and volumes of revolution appear in machine learning in several surprising contexts:

1. Gaussian Distributions in High Dimensions

The volume of a high-dimensional Gaussian "ball" is related to integrals of the form $\\int e^{-r^2} r^{n-1} dr$ . Understanding these volumes explains the "curse of dimensionality" — why most of a high-dimensional Gaussian's probability mass concentrates in a thin shell.

2. Kernel Density Estimation

When using radially symmetric kernels (like Gaussian or Epanechnikov), the normalization constants involve computing volumes of solids of revolution.

3. Monte Carlo Integration

Estimating volumes of complex solids using random sampling is a foundational technique in computational statistics. The disk method provides exact solutions for comparison when validating Monte Carlo methods.

4. 3D Computer Vision

Reconstructing 3D objects from 2D silhouettes often involves assuming rotational symmetry. The disk method helps compute volumes of these reconstructed shapes.

Python Implementation

Computing Volumes Numerically

Here's a complete Python implementation for computing volumes using the disk method:

Disk Method Implementation

🐍disk_method.py

Explanation(5)

Code(70)

3Disk Method Volume Function

This function computes the exact volume using numerical integration (Simpson's rule). The disk method formula V = π∫[f(x)]² dx requires integrating the squared function.

17Simpson's Rule Setup

Simpson's rule provides high accuracy by using parabolic interpolation. The weights alternate: 1, 4, 2, 4, 2, ..., 4, 1. This achieves O(h⁴) error convergence.

EXAMPLE

For n=4: weights = [1, 4, 2, 4, 1]

24Disk Approximation Function

This function approximates the volume using n discrete disks, similar to Riemann sums. Each disk has volume πr²Δx where r = f(x) at the midpoint.

36Computing Disk Radii

At each x position, the radius equals f(x). We use midpoint sampling for better accuracy, similar to the midpoint Riemann sum.

44Example: Rotating √x

The function y = √x rotated around the x-axis creates a paraboloid. The exact volume V = π∫x dx = 8π can be verified analytically.

65 lines without explanation

1import numpy as np
2import matplotlib.pyplot as plt
3from mpl_toolkits.mplot3d import Axes3D
4
5def disk_method_volume(f, a, b, n=1000):
6    """
7    Compute volume of solid of revolution using the disk method.
8
9    The disk method formula: V = π ∫_a^b [f(x)]² dx
10
11    Parameters:
12    - f: Function y = f(x) to rotate around x-axis
13    - a, b: Interval endpoints
14    - n: Number of subintervals for numerical integration
15
16    Returns: Approximate volume
17    """
18    # Simpson's rule for numerical integration
19    x = np.linspace(a, b, n + 1)
20    y_squared = f(x) ** 2
21
22    h = (b - a) / n
23    # Simpson's rule: (h/3)[f(x_0) + 4f(x_1) + 2f(x_2) + 4f(x_3) + ... + f(x_n)]
24    weights = np.ones(n + 1)
25    weights[1:-1:2] = 4  # Odd indices get weight 4
26    weights[2:-1:2] = 2  # Even indices (except first/last) get weight 2
27
28    integral = (h / 3) * np.sum(weights * y_squared)
29    volume = np.pi * integral
30
31    return volume
32
33def disk_approximation(f, a, b, n_disks):
34    """
35    Approximate volume using n discrete disks (Riemann sum).
36
37    Each disk has:
38    - Radius: r_i = f(x_i)
39    - Thickness: Δx = (b-a)/n
40    - Volume: π * r_i² * Δx
41
42    Total volume ≈ Σ π * [f(x_i)]² * Δx
43    """
44    dx = (b - a) / n_disks
45    x_midpoints = np.linspace(a + dx/2, b - dx/2, n_disks)
46    radii = f(x_midpoints)
47    disk_volumes = np.pi * radii**2 * dx
48
49    return np.sum(disk_volumes), x_midpoints, radii, dx
50
51# Example: Rotate y = √x from x = 0 to x = 4
52f = lambda x: np.sqrt(np.maximum(x, 0))
53
54# Exact volume: V = π ∫_0^4 x dx = π[x²/2]_0^4 = 8π
55exact_volume = 8 * np.pi
56
57print("Volume of y = √x rotated around x-axis from x = 0 to x = 4")
58print(f"Exact volume: V = 8π ≈ {exact_volume:.6f}")
59print()
60print(f"{'Disks':>8} {'Approx Volume':>16} {'Error':>12} {'% Error':>10}")
61print("-" * 50)
62
63for n in [4, 8, 16, 32, 64, 128, 256]:
64    approx, _, _, _ = disk_approximation(f, 0, 4, n)
65    error = abs(approx - exact_volume)
66    pct_error = 100 * error / exact_volume
67    print(f"{n:>8} {approx:>16.6f} {error:>12.6f} {pct_error:>9.4f}%")
68
69print()
70print(f"Simpson's rule (n=1000): {disk_method_volume(f, 0, 4, 1000):.10f}")

3D Visualization

Create stunning 3D visualizations of solids of revolution:

3D Visualization

🐍visualize_solid.py

Explanation(3)

Code(67)

33D Visualization Setup

Using matplotlib's 3D toolkit to visualize solids of revolution. We create two side-by-side plots: disk approximation and smooth surface.

15Disk Construction Loop

For each disk, we compute the position and radius. The disk is drawn as a filled circle at position x with radius r = f(x).

40Parametric Surface

The smooth solid is created using parametric equations: (x, r(x)cos(θ), r(x)sin(θ)) where r(x) = f(x). This creates the surface of revolution.

64 lines without explanation

1import numpy as np
2import matplotlib.pyplot as plt
3from mpl_toolkits.mplot3d import Axes3D
4
5def visualize_solid_of_revolution(f, a, b, n_disks=20):
6    """
7    Create a 3D visualization of a solid of revolution.
8    Shows both the approximating disks and the smooth surface.
9    """
10    fig = plt.figure(figsize=(14, 6))
11
12    # Left plot: Disk approximation
13    ax1 = fig.add_subplot(121, projection='3d')
14
15    dx = (b - a) / n_disks
16    theta = np.linspace(0, 2*np.pi, 50)
17
18    for i in range(n_disks):
19        x_left = a + i * dx
20        x_right = x_left + dx
21        x_mid = (x_left + x_right) / 2
22        r = f(x_mid)
23
24        # Create disk surface
25        X = np.array([[x_left, x_right], [x_left, x_right]])
26        for t in theta:
27            y = r * np.cos(t)
28            z = r * np.sin(t)
29            ax1.plot([x_left, x_right], [y, y], [z, z],
30                    color='green', alpha=0.3, linewidth=0.5)
31
32        # Draw disk edges
33        y_circle = r * np.cos(theta)
34        z_circle = r * np.sin(theta)
35        ax1.plot(np.full_like(theta, x_right), y_circle, z_circle,
36                'g-', linewidth=1.5, alpha=0.8)
37
38    ax1.set_title(f'Disk Approximation (n={n_disks})')
39    ax1.set_xlabel('x')
40    ax1.set_ylabel('y')
41    ax1.set_zlabel('z')
42
43    # Right plot: Smooth solid
44    ax2 = fig.add_subplot(122, projection='3d')
45
46    x = np.linspace(a, b, 100)
47    theta = np.linspace(0, 2*np.pi, 60)
48    X, T = np.meshgrid(x, theta)
49    R = f(X)
50    Y = R * np.cos(T)
51    Z = R * np.sin(T)
52
53    ax2.plot_surface(X, Y, Z, cmap='viridis', alpha=0.8,
54                     edgecolor='none', antialiased=True)
55
56    ax2.set_title('Smooth Solid of Revolution')
57    ax2.set_xlabel('x')
58    ax2.set_ylabel('y')
59    ax2.set_zlabel('z')
60
61    plt.tight_layout()
62    plt.savefig('solid_of_revolution.png', dpi=150)
63    plt.show()
64
65# Visualize y = √x rotated around x-axis
66f = lambda x: np.sqrt(np.maximum(x, 0))
67visualize_solid_of_revolution(f, 0, 4, n_disks=12)

Common Mistakes

Mistake 1: Forgetting to Square the Function

The most common error is writing $V = \\pi \\int f(x) \\, dx$ instead of $V = \\pi \\int [f(x)]^2 \\, dx$ . Remember: disk area is $\\pi r^2$ , so we need $[f(x)]^2$ !

Mistake 2: Wrong Variable of Integration

For x-axis rotation, integrate with respect to $dx$ . For y-axis rotation, express the curve as $x = g(y)$ and integrate with respect to $dy$ .

Mistake 3: Incorrect Bounds

The bounds of integration depend on which axis you're rotating around. For x-axis rotation, use x-values; for y-axis rotation, use y-values.

Mistake 4: Negative Radii

The radius must be non-negative. If $f(x) < 0$ on part of the interval, use $|f(x)|$ or $[f(x)]^2$ (which is always positive).

Test Your Understanding

Question 1 of 8Score: 0

When using the disk method to find the volume of a solid of revolution, what shape are the cross-sections perpendicular to the axis of rotation?

Summary

The disk method transforms the problem of finding 3D volumes into a 1D integration problem by slicing the solid into infinitesimally thin circular disks.

Key Formulas

Axis of Rotation	Formula	Notes
x-axis	V = π∫ₐᵇ [f(x)]² dx	y = f(x), integrate along x
y-axis	V = π∫ₓᵈ [g(y)]² dy	x = g(y), integrate along y

Key Concepts

Cross-sections are circles: When slicing perpendicular to the axis of rotation, each slice is a disk with radius $r = f(x)$
Disk volume is πr²dx: Each infinitesimally thin disk has volume $\\pi [f(x)]^2 dx$
Integration sums all disks: The definite integral adds up all disk volumes from $a$ to $b$
Square the function: Never forget the exponent 2 in $[f(x)]^2$
Match variables to axes: Use $dx$ for x-axis rotation, $dy$ for y-axis rotation

The Essence of the Disk Method:

"To find the volume of a solid of revolution, slice it into circular disks perpendicular to the axis, compute each disk's area (πr²), and integrate."

Coming Next: In the next section, we'll explore the Shell Method — an alternative approach that slices the solid into cylindrical shells instead of disks. This method is often more convenient when the axis of rotation is parallel to the axis of the function!