We discuss the existence and regularity of solutions to an elliptic system, whose basic example is \(\begin{aligned} {\left\{ \begin{array}{ll} \,-\Delta u =\frac{f(x)\hspace{0.55542pt}v^{\theta } }{u^{\theta }} & \text{ in }\;\; \Omega ,\\ \,-\Delta v = g(x)(1+u)^{\gamma } & \text{ in }\;\; \Omega ,\\ \,u>0,\;v>0 & \text{ in } \;\;\Omega ,\\ \,u = 0,\; v=0 & \text{ on } \;\; \partial \Omega , \end{array}\right. } \end{aligned}\) where \(0\leqslant f\leqslant g\in L^{m}(\Omega )\) are not identically zero, and \(\theta ,\gamma >0\) are fixed constants.