  Standard algorithms for the numerical integration of the Langevin equation
require that interactions are slowly varying during to the integration
timestep. This in not the case for hard-body systems, where there is no
clearcut between the correlation time of the noise and the timescale of the
interactions. Starting from a short time approximation of the Smoluchowsky
equation, we introduce an algorithm for the simulation of the overdamped
Brownian dynamics of polydisperse hard-spheres in absence of hydrodynamics
interactions and briefly discuss the extension to the case of external drifts.
