Abstract. We examine stable solutions of the following symmetric system on a complete, connected, smooth Riemannian manifold M without boundary,when ∆g stands for the Laplace-Beltrami operator,This system is called symmetric if the matrix of partial derivatives of all components of H, that is H(u) = (∂ j H i (u)) m i,j=1 , is symmetric. We prove a stability inequality and a Poincaré type inequality for stable solutions using the Bochner-Weitzenböck formula. Then, we apply these inequalities to establish Liouville theorems and flatness of level sets for stable solutions of the above symmetric system, under certain assumptions on the manifold and on solutions.