Mathematica finite element codes for the paper "Giant non-linear susceptibility of hydrogenic donors in silicon and germanium" by Nguyen H. Le et al.