The equilibrium magnetization of an array consisting of infinite rows containing 2N + 1 uniform Josephson junctions each is calculated from the dc Josephson equations with gauge-invariant phase differences. The calculated magnetization curves oscillate with modulated amplitudes. When the parameter α (=2Na0/λJ, where λJ is the Josephson penetration depth and a0 the lattice constant), is sufficently large, sweeping the external field up and down results in hysteresis loops. The differences from the conventional equilibrium curves for type-II superconductors are discussed. © 1994.