Computation of multi-region, relaxed magnetohydrodynamic equilibria with prescribed toroidal current profile