In the present work, by means of a phenomenological model, we simulate the hysteresis loop of an hexagonal array of Ni nanowires. Our model is based on the assumption that the hysteresis loop of a single wire is a rectangular box with a particular value of the coercive field, and the effect of the array is to generate a distribution of the coercive fields. Our results are in good agreement with experimental data.This work was supported by CONICYT Program Fellowship, FONDECYT-CHILE grant number 1040354, Millennium project P02-054-F.Peer reviewed