The equations of MHD equilibrium are solved for a large class of low-β stellarators with arbitrarily prescribed pressure and current density profiles by means of a large-aspect-ratio ordering scheme. The fields and flux surfaces are found to sufficiently high-order that surface shaping and axis shift are recovered; the relation between the latter and the applied vertical field is given. These results are applied to a discussion of several special cases.