The MRIYA program developed for calculation of transient one- ore two-dimensional temperature fields in any part of a nuclear power plant is described. The program efficiency is achieved by including the universal geometrical block and using the implicit Euler method in the thermal dynamics mathematical model differential equations system integration. The system integration is based on fast LU-decomposition algorithm for sparse system of equations solution. 9 refs.; 8 figs.; 4 tabs.