This report is devoted to the study of temperature dependence of bulk viscosity in SU(3)-gluodynamics. To calculate bulk viscosity wemeasured the correlation function of the trace anomaly for a set of temperatures in the region T/T_c \in (0.9, 1.5). We used multilevel algorithm which allowed us to improve the accuracy of the data. To extract the values of bulk viscosity we used two approaches: fitting of the data by physically motivated ansatz and the Backus-Gilbert method.

Nonzero Temperature and Density