User Rating 0.0 ★★★★★

Total Usage 0 times

Category Physics & Science

Initial Band Length (m) L₀ — length of the rubber band at t=0

Stretch Rate (m/s) v_s — speed at which the far end recedes

Bug Speed (m/s) v_b — constant walking speed of the bug

Time Compression 1× Simulation seconds per real second

Max Sim Time (s) Stop simulation after this many sim-seconds

Presets:

Time 0.00 s

Bug Position 0.00 m

Band Length 10.00 m

Fraction Complete 0.000%

Analytical Arrival —

Is this tool helpful?

Your feedback helps us improve.

★ ★ ★ ★ ★

About

The Bug-Rivet paradox (also known as the Ant on a Rubber Band) asks whether a bug walking at constant speed v_bug along a rubber band that stretches uniformly at rate v_stretch can ever reach the far end. Intuition suggests no - the band grows faster than the bug walks. Yet the bug's fractional progress per unit time forms terms of the harmonic series 1k, and since ∞∑k=1 1k diverges, the fraction eventually reaches 1. The bug always arrives, though the required time grows exponentially - for a 1 m band stretching at 1 m/s with a bug at 1 cm/s, arrival takes roughly e¹⁰⁰ − 1 seconds, vastly exceeding the age of the universe.

This simulator lets you configure all parameters and watch the paradox unfold in real time or accelerated mode. The continuous-time analytical solution is computed alongside a step-by-step numerical integration so you can verify convergence. Note: the model assumes perfectly uniform stretching and zero band mass. Real rubber bands break. The paradox is a mathematical statement about the harmonic series, not a materials science claim.

Formulas

The band has initial length L₀ and stretches uniformly at constant rate v_s. At time t, the band length is:

L(t) = L₀ + v_s ⋅ t

The bug at position x is carried along by the stretch and also walks at speed v_b. Its governing ODE:

dxdt = v_b + xL(t) ⋅ v_s

Defining the fractional position f(t) = x(t) ÷ L(t), the analytical solution is:

f(t) = v_bv_s ⋅ ln(1 + v_s ⋅ tL₀)

The bug reaches the end when f(t) = 1. Solving for arrival time T:

T = L₀v_s ⋅ (e^{v_s ÷ v_b} − 1)

In the discrete model (bug walks one step per second, band stretches after each step), the fractional progress after n steps equals:

f_n = v_bv_s ⋅ n∑k=1 1k

Since the harmonic series diverges, f_n eventually exceeds 1 for any positive v_b and v_s. This is the core of the paradox.

Where: L₀ = initial band length, v_s = stretch rate (speed at which the far end recedes), v_b = bug walking speed, x(t) = bug absolute position, f(t) = bug fractional position (0 to 1), T = arrival time.

Reference Data

Scenario	Band Length L₀	Stretch Rate v_s	Bug Speed v_b	Speed Ratio	Arrival Time t	Final Band Length
Classroom demo	1 m	1 m/s	0.5 m/s	2	6.39 s	7.39 m
Equal speeds	1 m	1 m/s	1 m/s	1	1.72 s	2.72 m
Slow bug	1 m	1 m/s	0.1 m/s	10	2.20 × 10⁴ s	2.20 × 10⁴ m
Very slow bug	1 m	1 m/s	0.01 m/s	100	e¹⁰⁰ − 1 s	≈ 2.69 × 10⁴³ m
Fast bug	10 m	2 m/s	3 m/s	0.67	4.87 s	19.74 m
Ant on a rope	1 km	1 km/s	1 cm/s	10⁵	e^10⁵ s	Astronomical
Snail on cosmos	1 ly	c	0.001 m/s	≈ 3 × 10¹¹	Finite but unimaginable	Beyond observable universe
No stretch	5 m	0 m/s	1 m/s	0	5 s	5 m
Micro scale	100 μm	10 μm/s	5 μm/s	2	63.9 s	739 μm
Ratio = 3	1 m	3 m/s	1 m/s	3	5.70 s	18.09 m
Ratio = 5	1 m	5 m/s	1 m/s	5	29.48 s	148.41 m
Ratio = 20	1 m	1 m/s	0.05 m/s	20	4.85 × 10⁸ s	4.85 × 10⁸ m
Harmonic series: 10 terms	10∑k=1 1k ≈ 2.9290					Partial sum
Harmonic series: 100 terms	100∑k=1 1k ≈ 5.1874					Partial sum

Frequently Asked Questions

The key insight is that stretching is uniform - every point on the band moves proportionally to its distance from the origin. The bug, already partway along the band, gets "carried forward" by the stretch. Its fractional progress each step is v_b / L(t), and since L(t) grows linearly, these fractions form terms proportional to 1/k. The harmonic series ∑(1/k) diverges, so the cumulative fraction inevitably reaches 1. The bug doesn't outrun the stretching - it accumulates fractional progress that sums to completion.

The arrival time is T = (L₀ / v_s) × (e^(v_s/v_b) − 1). This grows exponentially with the ratio v_s / v_b. At ratio 1, T ≈ 1.72 × L₀/v_s. At ratio 10, T ≈ 22,026 × L₀/v_s. At ratio 100, T involves e^100 ≈ 2.69 × 10^43. The function is mathematically finite for all finite ratios, but physically meaningless for large ones - the band would need to stretch beyond the observable universe.

Yes. The paradox is directly analogous to photons traveling through an expanding universe. A photon (the bug) moves at speed c through space that expands (the rubber band). In a universe with Hubble-like linear expansion (v_s proportional to distance), distant light can still reach us because the photon's fractional progress accumulates. However, if expansion accelerates (dark energy), the analogy breaks - the band equivalent would stretch exponentially, and the integral converges. In that regime, there exist galaxies whose light will never reach us, defining the cosmological event horizon.

Uniform stretching is essential to the paradox. If only the far end moves (like pulling a rope from one end), the bug gets no "free ride" - the portion behind the bug doesn't expand. In that case, the bug's progress is simply v_b / (L₀ + v_s × t), which integrates to a logarithm that diverges, so the bug still arrives (same math, different physics). However, if the band stretches exponentially - L(t) = L₀ × e^(αt) - the fractional progress integral converges, and the bug may never arrive if v_b < α × L₀.

The discrete model (stretch after each step) and the continuous model (simultaneous walking and stretching) give slightly different arrival times, but both guarantee arrival. The discrete harmonic sum H_n ≈ ln(n) + γ (where γ ≈ 0.5772 is the Euler-Mascheroni constant) approximates the continuous solution's ln(1 + v_s × t / L₀). For large speed ratios, the difference becomes negligible relative to the exponential arrival time. This simulator uses fourth-order Runge-Kutta integration for the continuous model, giving accuracy to within 10⁻⁸ of the analytical solution.

For ratios above about 20, the arrival time exceeds 10^8 in normalized units. The simulator caps visualization at a configurable maximum time to remain interactive. For large ratios, it computes and displays the analytical arrival time T = (L₀/v_s)(e^(v_s/v_b) − 1) directly, along with a log-scale progress chart showing f(t) approaching 1. The animation uses time compression - each frame can represent thousands of time units - letting you see the characteristic logarithmic growth curve even for extreme parameters.