Three-dimensional, pseudo-spectral computation is used to follow the evolution of a resistive, incompressible magnetofluid. The magnetofluid is confined by rigid, free-slip, perfectly-conducting square boundaries in the x, y directions (‘poloidal’ boundaries), and periodic boundary conditions are assumed in the z direction (‘toroidal’ direction). A constant, uniform d.c. magnetic field B0 is assumed in the z direction and a non-uniform current density j flows along it initially. Starting from a non-equilibrium hollow current profile, the evolution is followed for several tens of Alfvén transit times. Considerable small-scale turbulence develops, which causes energy to decay more rapidly than magnetic helicity. The average toroidal magnetic field at the (x, y) boundary reverses sign spontaneously. The near spatial constancy of the ratio jB/(jB) ≡ cos θ, in the relaxed state at late times, suggests that the state is nearly force-free. However, the ratio j. B/B2 ≡ α is considerably less uniform than is cos θ suggesting more residual disorder than a pure minimum-energy state would display.