Results from classical molecular dynamics simulations of infinite nuclear systems with varying density, temperature and isospin content are used to calculate the symmetry energy at low densities. The results show an excellent agreement with the experimental data and corroborate the claim that the formation of clusters has a strong influence on the symmetry energy in the liquid–gas coexistence region.