Why the random error vanishes under refinement
Since and , we have . Independence gives
The right-hand side tends to zero, proving convergence in mean square and therefore in probability at fixed . Along dyadic partitions the error probabilities are summable by Chebyshev's inequality, giving almost sure convergence for that sequence. For with constant , the squared sum is multiplied by , so .
For an uneven grid, let denote the longest interval. Each is at most , so the variance is bounded by . This tells you exactly which part of the sampling design matters. Adding a lot of points to the first half does little to control error if you leave one long interval in the second half. It also separates two different random quantities. The endpoint displacement is a single draw with variance , while the squared-step sum estimates from the whole path. On nested dyadic grids you can see the estimate improve overall without expecting each refinement to move closer to .
The grid experiment has another useful comparison. On a uniform grid with intervals, the variance of the estimate is , so doubling reduces its standard deviation by a factor of . To cut that standard deviation in half you need four times as many intervals. Move the slider and watch the measured sum fluctuate. The path and the theoretical value stay fixed. A coarser grid can land nearer by chance, so judge convergence by the shrinking error scale rather than by monotonicity.
How noisy is a finite-grid estimate?
Compare the same path across observation grids
import numpy as np
rng = np.random.default_rng(42)
horizon = 1.0
for steps in (25, 100, 400):
increments = rng.normal(0, np.sqrt(horizon / steps), (4_000, steps))
realized = np.sum(increments**2, axis=1)
print(steps, realized.mean(), realized.std(ddof=1))