We show that the saturation of resistive tearing modes in a cylindrical tokamak, as well as the corresponding island width, can be directly calculated with a magnetohydrodynamics (MHD) equilibrium code without solving the dynamics and without considering resistivity. The results are compared with initial-value resistive MHD simulations and to an analytical nonlinear theory. For small enough islands, the agreement is remarkable. For sufficiently large islands, the equilibrium calculations, which assume a flat current profile inside the island, overestimate the saturation amplitude. On the other hand, excellent agreement between nonlinear resistive MHD simulations and nonlinear theory is observed for all the considered tearing unstable equilibria.