Gradient estimates of mean curvature equation with Neumann boundary condition in domains of Riemannian manifold

We study the prescribed mean curvature equation with Neumann boundary conditions in domains of Riemannian manifold. The main goal is to establish the gradient estimates for solutions by the maximum principle. As a consequence, we obtain an existence result.