We propose the new system of equations for magnetodynamics. It is based on BBGKY hierarchy and uses a new approximations of mulit-particles distribution functions, which takes nearest-neighbour correlations into account. Thus we call the new system of equations correlational magnetodynamics system (CMD). It consists of well known Landau–Lifshitz–Bloch(LLB) equation complemented with an equation for two-particles correlations. Compared to traditional LLB equation, numerical simulations with CMD produces results that are closer to atomistic simulations.