Gillespie (1977): The exact stochastic simulation of master equations. Draw time τ from exponential distribution with rate Σaᵢ, choose reaction μ with prob aμ/Σaᵢ. No approximations — captures intrinsic molecular noise that ODEs miss. Crucial for gene expression, cell signaling with low copy numbers.