O presente trabalho teve como objetivo implementar a solução numérica da equação da difusão do calor, bidimensional, com condições de contorno de Dirichlet e em regime não permanente em malha generalizada. Foi utilizado o método de volumes finitos com uma formulação implícita, o sistema de equações lineares resultante foi resolvido pelo método de Gauss Seidel. Os resultados obtidos com as geometrias testadas indicaram que a metodologia de solução e a implementação numérica foram realizadas corretamente.