Let Omega be a bounded domain with a smooth C(2) boundary in R(N)(N >= 3), 0 is an element of (Omega) over bar, and n denote the unit outward normal to partial derivative Omega. We are concerned with the Neumann boundary problems: -div(vertical bar x vertical bar(alpha)vertical bar del u vertical bar(p-2)del u) = vertical bar x vertical bar(beta)u(p(alpha,beta)-1) - lambda vertical bar x vertical bar(gamma)u(p-1), u(x) > 0, x is an element of Omega, partial derivative u/partial derivative n = 0 on partial derivative Omega, where 1 < p < N and alpha < 0, beta < 0 such that p(alpha, beta) (sic) ...