An expansion method is developed for long-wavelength heliac and stellarator equilibrium. Time-dependent equations of motion are solved numerically to give equilibria which can include magnetic islands and stochastic field regions. Analytic low-beta equilibria are found which are in good agreement with computations.