On the numerical determination of ideal MHD limits of stability of axisymmetric toroidal configurations