This study concerns the existence of positive solution for the following nonlinear system $$\{-div(|x|^{-ap}|{\nabla}u|^{p-2}{\nabla}u)=|x|^{-(a+1)p+c_1}({\alpha}_1f(v)+{\beta}_1h(u)),x{\in}{\Omega},\\-div(|x|^{-bq}|{\nabla}v|q^{-2}{\nabla}v)=|x|^{-(b+1)q+c_2}({\alpha}_2g(u)+{\beta}_2k(v)),x{\in}{\Omega},\\u=v=0,x{\in}{\partial}{\Omega}$$, where ${\Omega}$ is a bounded smooth domain of $\mathbb{R}^N$ with $0{\in}{\Omega}$, 1 < $p,q$ < N, $0{{\leq}}a<\frac{N-p}{p}$, $0{{\leq}}b<\frac{N-q}{q}$ and $c_1$, $c_2$, ${\alpha}_1$, ${\alpha}_2$, ${\beta}_1$, ${\beta}_2$ are positive parameters. Here $f,g,h,k$ : $[0,{\infty}){\rightarrow}[0,{\infty})$ are nondecresing continuous functions and $$\lim_{s{\rightarrow}{\infty}}\frac{f(Ag(s)^{\frac{1}{q-1}})}{s^{p-1}}=0$$ for every A > 0. We discuss the existence of positive solution when $f,g,h$ and $k$ satisfy certain additional conditions. We use the method of sub-super solutions to establish our results.