Computation of zero β three-dimensional equilibria with magnetic islands