The iterative method of evaluating plasma equilibria is especially simple in a magnetic coordinate representation. This method is particularly useful for clarifying the subtle constraints of three‐dimensional equilibria and studying magnetic surface breakup at high plasma beta.