In this article, we consider the numerical method for solving the system of inequalities under the order induced by a symmetric cone with the function involved being monotone. Based on a perturbed smoothing function, the underlying system of inequalities is reformulated as a system of smooth equations, and a smoothing-type method is proposed to solve it iteratively so that a solution of the system of inequalities is found. By means of the theory of Euclidean Jordan algebras, the algorithm is

more »

... ved to be well defined, and to be globally convergent under weak assumptions and locally quadratically convergent under suitable assumptions. Preliminary numerical results indicate that the algorithm is effective.