Langevin Dynamics — Brownian Harmonic Oscillator

Overdamped: γẋ = −kx + √(2γkT)·η(t). Boltzmann distribution emerges.

⟨x⟩: 0.00 ⟨x²⟩: - kT/k theory: - MSD(t): - t: 0