In this work, we are concerned with the class of systems \(\begin{aligned} \left\{ \begin{array}{ll} -\Delta u=\lambda a(x)u^q+\tau c(x)u^{\alpha -1} v^\beta & \text{ in }~~\Omega \\ -\Delta v=\mu b(x)v^p+\delta c(x)u^\alpha v^{\beta -1} & \text{ in }~~\Omega \\ 0\not \equiv u\ge 0,\,\,0\not \equiv v\ge 0\,& \text{ in }~~\Omega \\ u=v=0& \text{ on }~~\partial \Omega \end{array} \right. \end{aligned}\) where \(\Omega \) is a smooth bounded domain in \({\mathbb {R}}^N\) , a, b, c belong to \(L^\infty (\Omega )\) , \(\lambda ,\,\mu ,\,\tau ,\,\delta >0\) , \(p,q\in (0,1)\) , \(\alpha ,\,\beta \ge 1\) and \(\alpha +\beta >2\) .