Boo-AI — Master Artificial Intelligence by Building from Scratch

Learning Objectives

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

Understand the slicing method for computing volumes of solids with known cross-sections
Derive the area function $A(x)$ for various cross-section shapes: squares, semicircles, triangles
Set up and evaluate volume integrals $V = \\int_a^b A(x)\\,dx$
Choose the appropriate axis of integration based on the problem geometry
Apply Cavalieri's Principle to compare volumes of different solids
Connect the cross-section method to numerical integration and 3D printing applications

The Big Picture: Building Solids Layer by Layer

"Every solid is a stack of infinitesimally thin slices, each contributing its area to the whole."

In previous sections, we computed volumes of solids of revolution using disks, washers, and shells. But what if our solid is not a surface of revolution? What if we simply know the shape of each cross-section?

This is precisely the situation in many real-world problems:

🏗️ Engineering

Bridge deck volumes with varying cross-section
Airplane fuselage design
Ship hull calculations
Structural beam analysis

🖨️ 3D Printing

Layer-by-layer construction
Material estimation
Print time calculation
Supports and infill volume

🏥 Medical Imaging

Organ volume from CT/MRI slices
Tumor size estimation
Blood vessel analysis
Prosthetic design

🎮 Computer Graphics

Constructive solid geometry
Volume rendering
Collision detection
Physics simulation

The Central Question

Given a solid where we know the cross-sectional area $A(x)$ at each position $x$ along an axis, how do we find the total volume?

The answer: integrate the area function. Each thin slice of thickness $dx$ contributes volume $A(x)\\,dx$ , and we sum all these contributions.

Historical Context: From Archimedes to Modern Slicing

Cavalieri and the Method of Indivisibles (1635)

The Italian mathematician Bonaventura Cavalieri (1598–1647) developed the method of indivisibles, a precursor to integral calculus. His fundamental insight: two solids with equal cross-sectional areas at every height have equal volumes.

This principle, now called Cavalieri's Principle, was revolutionary. It allowed comparison of volumes without computing them explicitly — if you can show two solids always have the same cross-sectional area, they must have the same volume.

Newton and Leibniz: Making It Rigorous

The cross-section method became fully rigorous with the development of integral calculus by Newton and Leibniz. What Cavalieri described intuitively — slicing and summing — became the precise operation of integration.

Modern Applications

The cross-section method is fundamental to computed tomography (CT). A CT scanner takes X-ray "slices" of the body, and software reconstructs the 3D volume by integrating these cross-sections — exactly the mathematical method we'll develop!

The Slicing Method: The General Formula

Consider a solid that extends from $x = a$ to $x = b$ along the x-axis. At each position $x$ , the cross-section perpendicular to the x-axis has area $A(x)$ .

Step 1: Slice the Solid

Divide the interval $[a, b]$ into $n$ subintervals, each of width $\\Delta x = (b-a)/n$ .

Step 2: Approximate Each Slice as a Slab

The $i$ -th slice is approximately a slab with:

Cross-sectional area $A(x_i^*)$ (evaluated at some sample point $x_i^*$ )
Thickness $\\Delta x$

The volume of this slab is approximately:

\\Delta V_i \\approx A(x_i^*) \\cdot \\Delta x

Step 3: Sum All Slabs

The total volume is approximately the Riemann sum:

V \\approx \\sum_{i=1}^{n} A(x_i^*) \\cdot \\Delta x

Step 4: Take the Limit

As $n \\to \\infty$ (equivalently, as $\\Delta x \\to 0$ ), the Riemann sum becomes a definite integral:

The Cross-Section Volume Formula

V = \\int_a^b A(x)\\,dx

The Key Insight

The formula $V = \\int_a^b A(x)\\,dx$ works for any solid where $A(x)$ is a known, integrable function. The cross-sections don't need to be circles — they can be squares, triangles, or any shape whose area you can express in terms of $x$ .

Cross-Section Area Formulas

To use the cross-section method, we need to express the area $A(x)$ as a function of position. Here are the common shapes:

Cross-Section Shape	If Side/Diameter/Base = s	Area Formula
Square	Side = s	A = s²
Rectangle (h × w)	Height h, width w	A = hw
Circle	Diameter = s	A = π(s/2)² = πs²/4
Semicircle	Diameter = s	A = π(s/2)²/2 = πs²/8
Equilateral Triangle	Side = s	A = (√3/4)s²
Isosceles Right Triangle	Legs = s	A = s²/2
Regular Hexagon	Side = s	A = (3√3/2)s²

Finding the Side Length s(x)

The crucial step is determining $s(x)$ — the characteristic dimension of the cross-section at position $x$ .

Common Scenarios:

Case 1: Cross-sections span between two curves $y = f(x)$ and $y = g(x)$

s(x) = |f(x) - g(x)|

Case 2: Cross-sections span the width of a region in the xy-plane

For a circular base $x^2 + y^2 = r^2$ :

s(x) = 2\\sqrt{r^2 - x^2}

Case 3: Side length given directly as a function

If the problem states "square cross-sections with side $s = 1 - x$ "

A(x) = (1-x)^2

Interactive 3D Cross-Section Explorer

Use the interactive visualization below to explore how different cross-section shapes create different solids. Adjust the number of slices to see how the Riemann sum approximation converges to the exact volume.

🎯3D Cross-Section Volume Explorer

Select Solid

Square Cross-Sections (Parabola Base)

Square cross-sections perpendicular to x-axis, with side equal to vertical distance between y = √x and y = -√x

Number of Slices: 10

Highlight Slice: None

Area Function

A(x) = (2\sqrt{x})^2 = 4x

Volume Calculation

V = \int_0^4 4x\,dx = 2x^2\Big|_0^4 = 32

Metric	Value
Riemann Sum Approximation (n = 10)	32.000000
Exact Volume (by integration)	32.000000
Error	0.000000 (0.00%)

Key Insight

The volume is found by "stacking" infinitely thin cross-sections. Each slice contributes $A(x) \cdot dx$ to the total volume. As slices become infinitely thin, the Riemann sum becomes the definite integral $V = \int_a^b A(x)\,dx$ .

Try This

• Drag to rotate the 3D view and see the cross-sections from different angles
• Increase the number of slices and watch the approximation improve
• Use the highlight slider to see individual cross-sections
• Switch between different solids to see how cross-section shape affects volume

Worked Examples

Example 1: Square Cross-Sections on a Parabolic Base

Problem: Find the volume of the solid whose base is the region bounded by $y = \\sqrt{x}$ and $y = -\\sqrt{x}$ for $0 \\leq x \\leq 4$ , and whose cross-sections perpendicular to the x-axis are squares.

Step 1: Identify the side length of each square cross-section.

At position $x$ , the square's side spans from $y = -\\sqrt{x}$ to $y = \\sqrt{x}$ :

s(x) = \\sqrt{x} - (-\\sqrt{x}) = 2\\sqrt{x}

Step 2: Compute the area function.

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

Step 3: Set up and evaluate the integral.

V = \\int_0^4 4x\\,dx = 2x^2\\Big|_0^4 = 2(16) - 0 = 32

Answer: The volume is $32$ cubic units.

Example 2: Semicircular Cross-Sections on a Circular Base

Problem: The base of a solid is the disk $x^2 + y^2 \\leq 4$ . Cross-sections perpendicular to the x-axis are semicircles with diameter along the base. Find the volume.

Step 1: Find the diameter of each semicircle.

At position $x$ , the circle $x^2 + y^2 = 4$ gives $y = \\pm\\sqrt{4 - x^2}$ .

\\text{Diameter} = 2\\sqrt{4 - x^2}, \\quad \\text{Radius} = \\sqrt{4 - x^2}

Step 2: Compute the semicircle area.

A(x) = \\frac{1}{2}\\pi r^2 = \\frac{\\pi}{2}(4 - x^2)

Step 3: Set up and evaluate the integral.

V = \\int_{-2}^{2} \\frac{\\pi}{2}(4 - x^2)\\,dx = \\frac{\\pi}{2}\\left[4x - \\frac{x^3}{3}\\right]_{-2}^{2}

= \\frac{\\pi}{2}\\left[\\left(8 - \\frac{8}{3}\\right) - \\left(-8 + \\frac{8}{3}\\right)\\right] = \\frac{\\pi}{2} \\cdot \\frac{32}{3} = \\frac{16\\pi}{3}

Answer: The volume is $\\frac{16\\pi}{3} \\approx 16.76$ cubic units.

Example 3: Equilateral Triangle Cross-Sections

Problem: Find the volume of the solid with base in the first quadrant bounded by $y = 1 - x$ , $x = 0$ , and $y = 0$ , and with equilateral triangle cross-sections perpendicular to the x-axis.

Step 1: The base of each triangle is the line segment from $y = 0$ to $y = 1 - x$ .

s(x) = 1 - x

Step 2: Area of equilateral triangle with side $s$ :

A(x) = \\frac{\\sqrt{3}}{4}(1-x)^2

Step 3: Evaluate the integral.

V = \\int_0^1 \\frac{\\sqrt{3}}{4}(1-x)^2\\,dx = \\frac{\\sqrt{3}}{4}\\left[-\\frac{(1-x)^3}{3}\\right]_0^1

= \\frac{\\sqrt{3}}{4} \\cdot \\frac{1}{3} = \\frac{\\sqrt{3}}{12} \\approx 0.144

Answer: The volume is $\\frac{\\sqrt{3}}{12} \\approx 0.144$ cubic units.

Cross-Sections Perpendicular to the y-axis

Sometimes it's more natural to take cross-sections perpendicular to the y-axis. In this case, we integrate with respect to $y$ :

Cross-Sections Perpendicular to y-axis

V = \\int_c^d A(y)\\,dy

Example 4: Square Cross-Sections Perpendicular to y-axis

Problem: The base is bounded by $y = x^2$ and $y = 1$ . Cross-sections perpendicular to the y-axis are squares. Find the volume.

Step 1: Express $x$ in terms of $y$ .

From $y = x^2$ , we get $x = \\pm\\sqrt{y}$ . The base spans from $-\\sqrt{y}$ to $\\sqrt{y}$ .

s(y) = 2\\sqrt{y}

Step 2: Area function.

A(y) = (2\\sqrt{y})^2 = 4y

Step 3: Integrate with respect to $y$ from 0 to 1.

V = \\int_0^1 4y\\,dy = 2y^2\\Big|_0^1 = 2

Answer: The volume is $2$ cubic units.

Cavalieri's Principle

Cavalieri's Principle provides a powerful way to compare volumes without computing integrals:

Cavalieri's Principle

If two solids have the same height, and at every level the cross-sectional areas are equal, then the two solids have equal volumes.

Mathematically: If $A_1(x) = A_2(x)$ for all $x \\in [a,b]$ , then $V_1 = V_2$ .

Classic Application: Sphere vs. Cylinder-Minus-Cone

Consider a hemisphere of radius $r$ sitting on a table, and a cylinder of radius $r$ and height $r$ with an inverted cone removed from it.

At height $h$ :

Hemisphere: Cross-section is a circle of radius $\\sqrt{r^2 - h^2}$ , area = $\\pi(r^2 - h^2)$
Cylinder minus cone: Outer circle has area $\\pi r^2$ , cone hole has radius $h$ , so net area = $\\pi r^2 - \\pi h^2 = \\pi(r^2 - h^2)$

Since the areas match at every height, by Cavalieri's Principle:

V_{\\text{hemisphere}} = V_{\\text{cylinder}} - V_{\\text{cone}} = \\pi r^3 - \\frac{1}{3}\\pi r^3 = \\frac{2}{3}\\pi r^3

This is exactly the known formula for a hemisphere — and we derived it without computing any integrals!

Real-World Applications

3D Printing: Layer-by-Layer Construction

3D printers build objects layer by layer, which is a direct physical implementation of the cross-section method. Each layer has a certain cross-sectional area, and the total material used is:

\\text{Material Volume} = \\sum_{i=1}^{n} A_i \\cdot \\Delta z \\approx \\int A(z)\\,dz

Slicing software uses exactly this principle to estimate print time and material consumption.

Medical Imaging: CT Scan Volume Reconstruction

A CT scanner takes hundreds of cross-sectional images (slices) of the body. To estimate the volume of an organ or tumor:

Identify the organ boundary in each slice (segmentation)
Calculate the area $A_i$ of the organ in slice $i$
Sum: $V \\approx \\sum A_i \\cdot \\Delta z$ where $\\Delta z$ is slice thickness

Civil Engineering: Earthwork Calculations

When designing roads, engineers must estimate the volume of earth to be moved (cut and fill). They survey cross-sections at intervals along the road alignment, then integrate to find total volume.

Machine Learning Connections

Volume Rendering in Deep Learning

Neural Radiance Fields (NeRF) and related 3D reconstruction methods use volume rendering, which involves integrating properties along rays:

C(\\mathbf{r}) = \\int_0^{\\infty} T(t) \\cdot \\sigma(\\mathbf{r}(t)) \\cdot \\mathbf{c}(\\mathbf{r}(t))\\,dt

This is mathematically identical to cross-section integration — we're summing contributions along a path through a 3D volume.

Probability Distributions: Marginalizing Over Variables

In Bayesian inference, integrating over a variable is analogous to summing cross-sections:

p(x) = \\int p(x, y)\\,dy

Each "slice" of the joint distribution at a fixed $y$ contributes to the marginal distribution — the same conceptual framework as volume by cross-sections.

Python Implementation

Computing Volumes Numerically

Cross-Section Volume Calculator

🐍cross_section_volume.py

Explanation(5)

Code(71)

3The Core Function

This function implements the slicing method: V = ∫A(x)dx. We pass in a function A_func that computes the cross-sectional area at any position x.

15Midpoint Rule

We use the midpoint rule for numerical integration because it typically gives better accuracy than left or right Riemann sums. Each x_mid is at the center of its subinterval.

19Summing Slices

The volume is approximated as the sum of thin slabs: V ≈ Σ A(xᵢ)·Δx. NumPy's vectorized operations make this efficient.

27Square Cross-Section Area

For the parabolic base where the side length is 2√x, the square area is (2√x)² = 4x. This gives a simple linear function of x.

40Semicircle Area

For semicircular cross-sections on a circular base, the diameter at position x is 2√(4-x²). The semicircle area is (π/2)r² = (π/2)(4-x²).

66 lines without explanation

1import numpy as np
2import matplotlib.pyplot as plt
3from mpl_toolkits.mplot3d import Axes3D
4
5def cross_section_volume(A_func, a, b, n=1000):
6    """
7    Compute volume of a solid with known cross-sectional area.
8
9    Parameters:
10    - A_func: Function that returns cross-sectional area at position x
11    - a, b: Bounds of integration
12    - n: Number of subintervals for numerical integration
13
14    Returns: Approximated volume
15    """
16    # Using the midpoint rule for better accuracy
17    dx = (b - a) / n
18    x_mid = np.linspace(a + dx/2, b - dx/2, n)
19
20    # Sum up A(x) * dx for each slice
21    areas = A_func(x_mid)
22    volume = np.sum(areas * dx)
23
24    return volume
25
26# Example 1: Square cross-sections on parabolic base
27# Base: y = sqrt(x), y = -sqrt(x), 0 <= x <= 4
28# Side of square = 2*sqrt(x)
29# A(x) = (2*sqrt(x))^2 = 4x
30
31def square_area(x):
32    """Area of square cross-section at position x"""
33    return 4 * x
34
35# Calculate volume
36volume_square = cross_section_volume(square_area, 0, 4)
37print(f"Square cross-sections on parabola:")
38print(f"  Numerical: {volume_square:.6f}")
39print(f"  Exact:     {32:.6f}")
40print()
41
42# Example 2: Semicircular cross-sections on circular base
43# Base: x^2 + y^2 = 4 (circle of radius 2)
44# Diameter = 2*sqrt(4 - x^2), radius = sqrt(4 - x^2)
45# A(x) = (1/2) * pi * r^2 = (pi/2) * (4 - x^2)
46
47def semicircle_area(x):
48    """Area of semicircular cross-section at position x"""
49    return (np.pi / 2) * (4 - x**2)
50
51volume_semicircle = cross_section_volume(semicircle_area, -2, 2)
52exact_semicircle = 16 * np.pi / 3
53print(f"Semicircular cross-sections on circle:")
54print(f"  Numerical: {volume_semicircle:.6f}")
55print(f"  Exact:     {exact_semicircle:.6f}")
56print()
57
58# Example 3: Equilateral triangle cross-sections
59# Base: region bounded by y = 1 - x and axes in first quadrant
60# Side = 1 - x
61# A(x) = (sqrt(3)/4) * (1-x)^2
62
63def equilateral_triangle_area(x):
64    """Area of equilateral triangle cross-section"""
65    return (np.sqrt(3) / 4) * (1 - x)**2
66
67volume_triangle = cross_section_volume(equilateral_triangle_area, 0, 1)
68exact_triangle = np.sqrt(3) / 12
69print(f"Equilateral triangle cross-sections:")
70print(f"  Numerical: {volume_triangle:.6f}")
71print(f"  Exact:     {exact_triangle:.6f}")

3D Visualization

3D Cross-Section Visualization

🐍visualize_cross_sections.py

Explanation(4)

Code(65)

8Two-Panel Visualization

We create two subplots: a 3D view of the actual solid and a 2D plot of the area function A(x). This helps students connect the geometry to the integration.

23Building Each Slice

For each x position, we calculate the cross-section dimensions and create a polygon in the y-z plane. The collection of these polygons forms the solid.

303D Collection

Matplotlib's Poly3DCollection allows us to add 3D polygons to the plot. Setting alpha=0.3 makes them semi-transparent so we can see through the solid.

45Area Function Plot

The right panel shows A(x) vs x. The area under this curve equals the volume of the solid — a perfect visual connection to V = ∫A(x)dx.

61 lines without explanation

1import numpy as np
2import matplotlib.pyplot as plt
3from mpl_toolkits.mplot3d import Axes3D
4from mpl_toolkits.mplot3d.art3d import Poly3DCollection
5
6def visualize_cross_sections(shape='square', n_slices=20):
7    """
8    Create 3D visualization of a solid with known cross-sections.
9    """
10    fig = plt.figure(figsize=(12, 5))
11
12    # Left: 3D view of the solid
13    ax1 = fig.add_subplot(121, projection='3d')
14
15    # Right: Cross-sectional area function
16    ax2 = fig.add_subplot(122)
17
18    if shape == 'square':
19        # Square cross-sections on y = sqrt(x), y = -sqrt(x)
20        x_vals = np.linspace(0.01, 4, n_slices)
21
22        for x in x_vals:
23            # Side length at this x
24            s = 2 * np.sqrt(x)
25            half_s = s / 2
26
27            # Square vertices (in y-z plane at position x)
28            verts = [
29                [x, -half_s, 0],
30                [x, half_s, 0],
31                [x, half_s, s],
32                [x, -half_s, s]
33            ]
34            ax1.add_collection3d(Poly3DCollection(
35                [verts], alpha=0.3,
36                facecolor='blue', edgecolor='darkblue'
37            ))
38
39        # Plot the base curves
40        x_curve = np.linspace(0, 4, 100)
41        ax1.plot(x_curve, np.sqrt(x_curve), np.zeros_like(x_curve),
42                 'r-', linewidth=2, label='y = √x')
43        ax1.plot(x_curve, -np.sqrt(x_curve), np.zeros_like(x_curve),
44                 'r-', linewidth=2, label='y = -√x')
45
46        # Area function
47        x_area = np.linspace(0, 4, 100)
48        A_x = 4 * x_area
49        ax2.fill_between(x_area, A_x, alpha=0.3, color='blue')
50        ax2.plot(x_area, A_x, 'b-', linewidth=2)
51        ax2.set_xlabel('x')
52        ax2.set_ylabel('A(x) = 4x')
53        ax2.set_title('Cross-Sectional Area Function')
54
55    ax1.set_xlabel('X')
56    ax1.set_ylabel('Y')
57    ax1.set_zlabel('Z')
58    ax1.set_title(f'Solid with {shape.title()} Cross-Sections')
59
60    plt.tight_layout()
61    plt.savefig(f'cross_section_{shape}.png', dpi=150)
62    plt.show()
63
64# Generate visualizations
65visualize_cross_sections('square', n_slices=15)

Common Pitfalls

Pitfall 1: Choosing the Wrong Axis

Make sure your cross-sections are perpendicular to the axis of integration. If cross-sections are perpendicular to the x-axis, integrate with respect to $x$ . If perpendicular to the y-axis, integrate with respect to $y$ .

Pitfall 2: Confusing Side Length with Area

Don't forget to square (or apply the correct area formula) to convert from side length to area. If $s(x) = 2\\sqrt{x}$ , then for squares, $A(x) = 4x$ , not $2\\sqrt{x}$ .

Pitfall 3: Incorrect Bounds of Integration

The bounds must correspond to where the solid starts and ends along the axis of integration. For a circular base $x^2 + y^2 = 4$ with cross-sections perpendicular to x-axis, the bounds are $x = -2$ to $x = 2$ .

Check Your Setup

Before integrating, verify:

Which axis are cross-sections perpendicular to?
What is the shape of each cross-section?
What dimension varies with position (side, diameter, base)?
What are the correct bounds of integration?

Test Your Understanding

❓Test Your Understanding

Question 1 of 8

What is the general formula for finding volume using cross-sections perpendicular to the x-axis?

Summary

The cross-section method provides a general framework for computing volumes when we know the shape of each cross-section. This approach is more flexible than disks/washers/shells because it applies to solids that are not surfaces of revolution.

Key Concepts

Concept	Description
Slicing Method	Divide solid into thin slabs, sum contributions A(x)·dx
Volume Formula	V = ∫ A(x) dx where A(x) is cross-sectional area
Area Functions	Square: s², Circle: πs²/4, Semicircle: πs²/8, Eq. Triangle: (√3/4)s²
Cavalieri's Principle	Equal cross-sections at all heights ⟹ equal volumes
Integration Axis	Cross-sections perpendicular to x-axis ⟹ integrate dx

Problem-Solving Strategy

Identify the base region and the shape of cross-sections
Determine the dimension that varies (side length, diameter, etc.)
Express the area $A(x)$ as a function of position
Set up the integral with correct bounds
Evaluate using integration techniques or numerical methods

The Essence of Cross-Sections:

"Any solid is just a stack of infinitely thin slices. Know the slice shape, integrate the area, and you have the volume."

Coming Next: In the next section, we'll explore Arc Length — using integration to measure the length of curves, another fundamental application of "slicing and summing."

Learning Objectives

The Big Picture: Building Solids Layer by Layer

🏗️ Engineering

🖨️ 3D Printing

🏥 Medical Imaging

🎮 Computer Graphics

The Central Question

Historical Context: From Archimedes to Modern Slicing

Cavalieri and the Method of Indivisibles (1635)

Newton and Leibniz: Making It Rigorous

Modern Applications

The Slicing Method: The General Formula

Step 1: Slice the Solid

Step 2: Approximate Each Slice as a Slab

Step 3: Sum All Slabs

Step 4: Take the Limit

The Key Insight

Cross-Section Area Formulas

Finding the Side Length s(x)

Case 1: Cross-sections span between two curves y=f(x)y = f(x)y=f(x) and y=g(x)y = g(x)y=g(x)

Case 2: Cross-sections span the width of a region in the xy-plane

Case 3: Side length given directly as a function

Interactive 3D Cross-Section Explorer

Square Cross-Sections (Parabola Base)

Try This

Worked Examples

Example 1: Square Cross-Sections on a Parabolic Base

Example 2: Semicircular Cross-Sections on a Circular Base

Example 3: Equilateral Triangle Cross-Sections

Cross-Sections Perpendicular to the y-axis

Example 4: Square Cross-Sections Perpendicular to y-axis

Cavalieri's Principle

Classic Application: Sphere vs. Cylinder-Minus-Cone

Real-World Applications

3D Printing: Layer-by-Layer Construction

Medical Imaging: CT Scan Volume Reconstruction

Civil Engineering: Earthwork Calculations

Machine Learning Connections

Volume Rendering in Deep Learning

Probability Distributions: Marginalizing Over Variables

Python Implementation

Computing Volumes Numerically

3D Visualization

Common Pitfalls

Pitfall 1: Choosing the Wrong Axis

Pitfall 2: Confusing Side Length with Area

Pitfall 3: Incorrect Bounds of Integration

Check Your Setup

Test Your Understanding

What is the general formula for finding volume using cross-sections perpendicular to the x-axis?

Summary

Key Concepts

Problem-Solving Strategy

Case 1: Cross-sections span between two curves $y = f(x)$ and $y = g(x)$