In this paper, we consider the following elliptic systems involving critical Sobolev growth and Hardy potential: -Δu-λ1u|x|2= a1|u|2*-2u+bh(x)αα+β|u|α-2u|v|β,x∈RN, -Δv-λ2v|x|2=a2|v|2*-2v+bh(x)βα+β|u|α|v|β-2v,x∈RN, where N ≥3,λ1,λ2 ∈[0,ΛN), ΛN:=N-222 is the best Hardy constant. 2*=2NN-2 is the critical Sobolev exponent. a1,a2, b are positive parameters, α,β>0 and 1 <α+ β: = q <2 <2*. h(x)∈Lq′(RN) with q′=2*2*-q. By mean...