A direct method for the numerical computation of bifurcation points underlying symmetries