In this manuscript, based on the most widespread fixed point theories in literature. The existence of solutions to the system of nonlinear fractional differential equations with Caputo Hadmard fractional operator in a bounded domain is verified by using Mönoch’s fixed point theorem, The stability of the coupled system is also investigated via Ulam-Hyer technique. Finally, an applied numerical example is presented to illustrate the theoretical results obtained.