The Modified Local Green’s Function Method (MLGFM) is an integral method which uses appropriately chosen Green’s function projections obtained numerically with the aid of auxiliary finite element problems. Its applicability includes those cases for which a fundamental solution does not exist or is very cumbersome. The MLGFM was studied intensely in the 90´s with promising results, especially for tractions and heat fluxes at the boundaries. The present contribution compares this method for heat flux evaluation in anisotropic media with finite volumes and finite elements. The latter approximates heat fluxes using a superconvergent patch recovery scheme, whereas the former computes flux quantities directly at nodes. The numerical example uses linear elements and includes non-homogeneous temperature and flux boundary conditions.