A program for the solution of the two-dimensional diffusion equation by the finite element method was developed and applied to realistic reactor calculations. Comparisons were made with well-known diffusion theory codes and with nodal programs. Results of these comparisons and first experiences with the finite element method and its success in our program FEM 2D are reported.