Se trata la aplicación de un método variacional de elementos finitos (Raleigh Ritz) a la resolución de la ecuación de la difusión de neutrones en un medio absorbente con fuente distribuida. Se analizan con cierto detalle los problemas de establecimiento del funcional, generación de malla y almacenamiento de la matriz. Se ofrecen los resultados obtenidos para algunos casos concretos de geometría rectangular y cilíndrica, y comparación de las distribuciones teóricas que se pretendían aproximar.The paper deals with the application of a finite element variational method to the solution of the equation of neutron diffusion in an absorbing media with distributed source. The problems of settting up the functional, lattice generation, and matrix storage, are analysed in some detail. The results obtained for some particular cases of rectangular and cylindric geometry are shown and also a comparison with the theoretical distributions that were tried to approximate.