  We propose a Monte Carlo algorithm to search for the boundary points of the
orbit space which is important in determining the symmetry breaking directions
in the Higgs potential and the Landau potential. Our algorithm is robust,
efficient and generally applicable. We apply the method to the Landau potential
of the $d-$wave abnormal superconductor.
