<div>In this article, we consider the following p-q-Laplacian system with
singular and critical nonlinearity \begin{equation*}
\left \{
\begin{array}{lllll}
-\Delta_{p}u-\Delta_{q}u=\frac{h_{1}(x)}{u^{r}}+\lambda\frac{\alpha}{\alpha+\beta}u^{\alpha-1}v^{\beta}
\ \ in\
\Omega ,\\
-\Delta_{p}v-\Delta_{q}v=\frac{h_{2}(x)}{v^{r}}+\lambda\frac{\beta}{\alpha+\beta}u^{\alpha}v^{\beta-1}
\ \ in\
\Omega, \\
u,v&gt;0 \ \ \
\ \ \ in \
\Omega, \ \
\ \ \ u=v=0
\ \ \ \
\ \ \ on \
\partial\Omega,
\end{array} \right.
\end{equation*} where Ω is a bounded domain in
$\mathbb {R}^{n}$ with smooth boundary
$\partial\Omega$.
$11,\lambda\in(0,\Lambda_{*})$
is parameter with $\Lambda _{*}$ is a positive
constant and $h_{1}(x),h_{2}(x)\in
L^{\infty},h_{1}(x),h_{2}(x)&gt;0$.
We show the existence and multiplicity of weak solution of equation
above for suitable range of $\lambda$.</div>