The ordinary differential equations of Newtonian dynamics are used in atomic simulations with the method of molecular dynamics. The basic issues are surveyed and standard algorithms are described. Several algorithmic variants are discussed. Some advanced ideas relating to parallel computation are considered.