We investigate the effects of magnetic inhomogeneities and thermal fluctuations on the magnetic properties of a rare earth intermetallic compound, Nd2Fe14B. The constrained Monte Carlo method is applied to a Nd2Fe14B bulk system to realize the experimentally observed spin reorientation and magnetic anisotropy constants K A m (m = 1, 2, 4) at finite temperatures. Subsequently, it is found that the temperature dependence of K A 1 deviates from the Callen-Callen law,3 , even above room temperature, TR ∼ 300 K, when the Fe (Nd) anisotropy terms are removed to leave only the Nd (Fe) anisotropy terms. This is because the exchange couplings between Nd moments and Fe spins are much smaller than those between Fe spins. It is also found that the exponent n in the external magnetic field Hext response of barrier height FB = F 0 B (1 − Hext/H0) n is less than 2 in the low-temperature region below TR, whereas n approaches 2 when T > TR, indicating the presence of Stoner-Wohlfarth-type magnetization rotation. This reflects the fact that the magnetic anisotropy is mainly governed by the K A 1 term in the T > TR region.