Resolve-se a Equacao de Difusao de Neutrons a dois grupos de energia aplicando-se o Metodo dos Elementos Finitos (MEF), com elementos triangulares de primeira ordem. Para tanto foi desenvolvido o programa EFTDN (Elementos Finitos Triangulares na Difusao de Neutrons) em linguagem FORTRAN IV. A formulacao discreta da Equacao de Difusao e obtida aplicando o metodo de Galerkin. Para resolver o problema de auto-valores e utilizado o metodo da potencia e com a finalidade de acelerar a convergencia dos resultados aplica-se os polinomios de Chebshev. Na solucao dos sistemas de equacoes e aplicado o metodo de Gauss dividido em duas etapas distintas: triangularizacao da matriz dos coeficientes e a retrosubstituicao levando-se em conta a esparsidade do sistema. Alguns problemas teste sao resolvidos, dentre estes dois reatores P.W.R., o ZION-1 de 1300 MWe e o reator 2D-IAEA - Benchmark. Comparando os resultados com as solucoes padroes demonstra-se a validade da aplicacao do MEF e a precisao dos resultados. Discute-se ainda de que forma pode ser melhorado e ampliado este trabalho. (Autor)