Abstract: In this work, the magnetic properties of superparamagnetic iron oxide nanoparticles (SPIONs) submitted to an external magnetic field are studied using a Metropolis algorithm. The influence on the M(B) curves of the size distribution of the nanoparticles, of uniaxial anisotropy, and of dipolar interaction between the cores are examined, as well as the influence of drying the samples under a zero or non-zero magnetic field. It is shown that the anisotropy impacts the shape of the magnetization curves, which then deviate from a pure Langevin behaviour, whereas the dipolar interaction has no influence on the curves at 300 K for small particles (with a radius of 3 nm). The fitting of the magnetization curves of particles with magnetic anisotropy to a Langevin model (including a size distribution of the particles) can then lead to erroneous values of the distribution parameters. The simulation results are qualitatively compared to experimental results obtained for iron oxide nanoparticles (with a 3.21 nm median radius).