We consider the Neumann problem for the Schrödinger equations
$-\Delta u\,+\,Vu\,=\,0$ , with singular nonnegative potentials
$V$ belonging to the reverse Hölder class
${{\mathcal{B}}_{n}}$ , in a connected Lipschitz domain
$\Omega \,\subset \,{{\text{R}}^{n}}$ . Given boundary data
$g$ in
${{H}^{p}}\text{or}\,{{L}^{p}}\,\text{for}\,\text{1}-\in \,<\,p\,\le \,2,\text{where}\,\text{0}<\in <\frac{1}{n}$ , it is shown that there is a unique solution,
$u$ , that solves the Neumann problem for the given data and such that the nontangential maximal function of
$\nabla u$ is in
${{L}^{p}}(\partial \Omega )$ . Moreover, the uniform estimates are found.