The improved confinement of modern stellarators is largely attributed to the numerical optimization of the magnetic field. With the technique pioneered by Nührenberg and Zille [1]}, the shape of the plasma boundary, \(S_P\) , is optimized for certain properties of the magnetohydrodynamic (MHD) equilibrium,
\(\left(\nabla \times \textbf {B}\right) \times \textbf {B} &= \mu _0\nabla p &\text{in } V_P\\\textbf {B} \cdot \hat{\textbf {n}} &= 0 &\text{on } S_P,\) 
