A numerical study of surface- and diffuse-current, high-beta stellarator equilibria is described. Numerical results are obtained from an algorithm for the solution of the ideal magnetohydrodynamic equations in three dimensions as an initial- and boundary-value problem. Equilibria are obtained for diffuse-current plasmas just as for surface-current plasmas. The equilibrium conditions for diffuse-current plasmas are quantitatively different from those for surface-current plasmas, especially at high beta for ℓ = 0, 1 systems. In contrast, the equilibrium conditions scale according to small-δ theory even to δ ∼ 1. Results are obtained for a range of dimensionless parameters (beta, aspect ratio, etc.) and compared to an asymptotic equilibrium theory and experiments. Good agreement among the two theories and the experiments is shown.