Numerical computation of toroidal equilibria with anisotropic pressure in the presence of a dominant vacuum toroidal field