Ways have been found to significantly upgrade the resolution of an equilibrium and stability code that is based on the variational principle of ideal magnetohydrodynamics. Nonlinear saturation of ballooning modes for tokamaks has been demonstrated using a revised version of the code. Second stability regions are calculated from the nonlinear dependence of δW on β. An l=2,3 torsatron is described that appears to have no ideal magnetohydrodynamic instabilities for β as high as 5%.