Abstract. We consider the following non-local elliptic boundary value problem:where α and λ are positive constants and f is a function satisfyingThe solution of the equation represents the steady state of a thermistor device. The problem has a unique solution for a critical value λ * of the parameter λ, at least two solutions for λ < λ * and has no solution for λ > λ * . We apply a finite element and a finite volume method in order to find a numerical approximation of the solution of the problem from the space of continuous piecewise quadratic functions, for the case that λ < λ * and for the stable branch of the bifurcation diagram. A comparison of these two methods is made regarding their order of convergence for f (s) = e −s and f (s) = (1 + s) −2 . Also, for the same equation but with Dirichlet boundary conditions, a situation where the solution is unique for λ < λ * , a similar comparison of the finite element and the finite volume method is presented.