A numerical technique for obtaining the complete bifurcated equilibrium solution space for a tokamak