Suppose u is a harmonic function on a domain D and x, x' are in D. We estimate |u(x) — u(x')| using two Brownian motions started at x and x' and killed on exiting a cube Q ⊂ D. By selecting appropriate versions of the two Brownian motions, a classical gradient estimate for u is easily derived.