We present a finite volume discretization of the nonlinear elliptic problems. The discretization results in a nonlinear algebraic system of equations. A Newton-Krylov algorithm is also presented for solving the system of nonlinear algebraic equations. Numerically solving nonlinear partial differential equations consists of discretizing the nonlinear partial differential equation and then solving the formed nonlinear system of equations. We demonstrate the convergence of the discretization scheme and also the convergence of the Newton solver through a variety of practical numerical examples.