“…The numerous underlying problems were a source of motivation for many researchers. In order to solve (1.1) in the standard case, ν is a function, many aspects have been developed : Hardy inequality, see the references [5,15,19,21,28,29], variational methods [16,20] for example, the singularities of semilinear Hardy elliptic equation by [7,9,14,22,23,24,26,30] and the references therein. When N ≥ 3, µ 0 ≤ µ < 0, Dupaigne in [15] (also see [4,5]) studied the weak solutions for the equation L µ u = u p + tf in Ω, subject to the zero Dirichlet boundary condition, with t > 0, p > 1 and f ∈ L 1 (Ω), in the distributional sense that u ∈ L 1 (Ω), u p ∈ L 1 (Ω, ρdx) and…”