In this paper, we consider second order neutrons diffusion problem withcoefficients in L ∞(Ω). Nodal method of the lowest order is applied to approximate the problem's solution. The approximation uses special basis functions [1] in which the coefficientsappear. The rate of convergence obtained is O(h2) in L 2(Ω), with a free rectangular triangulation.