This paper deals with the Neumann problem for an elliptic equation {-△u-μu/|x|^2=|u|^2*(s)-2u/|x|^s+λ|u|^q-2u,x∈Ω, Dγu+a(x)u=0, x∈偏dΩ/{0}, where Ω is a bounded domain in R^N with C^1 boundary, 0 ∈ 偏dΩ, N 〉 5. 2^*(s) = 2(N-s)/N-2 (0 ≤ s ≤ 2) is the critical Sobolev-Hardy exponent, 1 〈 q 〈 2, 0 〈 μ 〈 μ^*, γ denotes the unit outward normal to boundary 偏dΩ. By variational method and the dual fountain theorem, the existence of infinitely many solutions with nega...