A rapid, stable, and simple triangulation method is presented for approximating equi-valued surfaces by polyhedra. It is applicable to various scientific and engineering fields, because equi-valued surfaces can be closed or open, and because spatial functions can be given by a parameter set of defining functions or by the numerical values of functions at grid points. The method is tested numerically to compare its accuracy, computation time, and data volume with those of other methods.