We study positive solutions to nonlinear elliptic systems of the form: \begin{eqnarray*} -\Delta u =\lambda f(v) \mbox{ in }\Omega\\-\Delta v =\lambda g(u) \mbox{ in }\Omega\\\quad~~ u=0=v \mbox{ on }\partial\Omega \end{eqnarray*} where $\Delta u$ is the Laplacian of $u$, $\lambda$ is a positive parameter and $\Omega$ is a bounded domain in $R^n$ with smooth boundary $\partial\Omega$. In particular, we will analyze the combined effects of the nonlinearities on the existence and multiplicity of positive solutions. We also study systems with multiparameters and stronger coupling. We extend our results to $p$-$q$-Laplacian systems and to $n\times n$ systems. We mainly use sub- and super-solutions to prove our results.