By a transformation of
metric using the first eigenfunction, we obtain lower bounds for all eigenvalues of
Schrödinger operator with the Neumann boundary condition. Global estimates for
first eigenfunction are needed and this is achieved by the Moser iteration
technique.