An exact analytic formula for the magnetostatic interaction energy of a three-dimensional array of ferromagnetic cubes is presented. The magnetization within each cube is assumed constant but varies ...