We explain how periodic boundary conditions are applied to molecular dynamics simulations in order to mimic the behaviour of macroscopic systems. We pay special attention to the computation of the long-range electrostatic interactions by introducing the Ewald summation method in a simple way. Thermostats and barostats aimed to control the temperature and pressure are shortly described.