“…In designing a poloidal-field system, an expeditious approach is to use the approximate Shafranov formulae to compute the poloidal field at the shell for a desired plasma equilibrium. From this field, the required winding distribution is determined [1]. Since, however, the Shafranov formulae use expansions in e, the inverse toroidal aspect ratio, questions arise regarding their accuracy as a function of plasma beta, plasma inductance, and toroidal aspect ratio.…”