In this article, we consider the following p‐q‐Laplacian system with singular and critical nonlinearity
−normalΔpu−normalΔqu=h1false(xfalse)ur+λαα+βuα−1vβ0.1em0.1emin0.4emnormalΩ,−normalΔpv−normalΔqv=h2false(xfalse)vr+λβα+βuαvβ−10.1em0.1emin0.4emnormalΩ,u,v>00.1em0.1em0.1em0.1em0.1em0.1emin0.4emnormalΩ,0.1em0.1em0.1em0.1em0.1emu=v=00.1em0.1em0.1em0.1em0.1em0.1em0.1emon0.4em∂normalΩ,
where Ω is a bounded domain in
ℝn with smooth boundary ∂Ω. 1 < q < p < α + β = p∗, 0 < r < 1, α, β > 1, λ ∈ (0, Λ∗) is parameter and h1(x), h2(x) ∈ L∞, h1(x), h2(x) > 0. We show the existence and multiplicity of weak solution of equation above for suitable range of λ.