In this work, we introduce multipoint flux (MF) approximation method to the problem of conduction heat transfer in anisotropic media. In such media, the heat flux vector is no longer coincident with the temperature gradient vector. In this case, thermal conductivity is described as a second order tensor that usually requires, at least, six quantities to be fully deflned in general three-dimensional problems. The tvt'o-point flux ßnite differences approximation may not handle such anisotropy and essentially more points need to be involved to describe the heatflta vector. In the framework of mixed flnite element method (MFE), the MFMFE methods are locally conservative with continuous normal fluxes. We consider the lowest order Brezzi-Douglas-Marini (BDM) mixed flnite element method with a special quadrature rule that allows for nodal velocity elimination resulting in a cell-centered system for the temperature. We show comparisons with some analytical solution of the problem of conduction heat transfer in anisotropic long strip. We also consider the problem of heat conduction in a bounded, rectangular domain with different anisotropy scenarios. It is noticed that the temperature fleld is signiflcantly affected by such anisotropy scenarios. Also, the technique used in this work has shown that it is possible to use the flnite difference settings to handle heat transfer in anisotropic media. In this case, heatflta vectors, for the case of rectangular mesh, generally require six points to be described.