In 1978, Ermak and McCammon suggested an algorithm for efficiently computing Brownian dynamics with hydrodynamic interactions. Hydrodynamic interactions occur when the particles interact indirectly by generating and reacting to local velocities in the solvent. For a system of N three-dimensional particle diffusing subject to a force vector F(X), the derived Brownian dynamics scheme becomes: