The action for four-dimensional Einstein-power-Yang-Mills (EPYM) gravity with a cosmological constant \(\Lambda \)  was given by [1]}, [2]}, [3]}, [4]}
\(I=\frac{1}{2}\int d^4x\sqrt{g}\left(R-2\Lambda -[Tr(F^{(a)}_{{\mu \nu }}F^{{(a)\mu \nu }})]^\gamma \right)\) 
