The subject matter of the article is the method of determining comprehensive indicator of labor safety. The goal of the work is to develop a methodology for obtaining a comprehensive assessment of labor safety, taking into account the characteristics of harmful and dangerous production factors. The article solves the following task: to investigate the possibility of applying the method of numerical integration to obtain a comprehensive indicator of occupational safety, that is, to apply the principle of determining the area and volume under curved surfaces both on the plane and in space, which are created by combining estimates of unit indicators of harmful factors on a dimensionless scale. Methods are used: qualimetry and mathematical statistics. The following results were obtained: a number of modern scientific works were considered, substantiating the need to assess occupational safety at the workplace. A number of modern scientific works related to the quantitative assessment of the quality of qualimetric objects of different nature, which have different quality indicators and different measurement scales, have been analyzed. As a result of the analysis, the relevance of the topic has been proven and the need to develop a methodology for obtaining a comprehensive assessment of labor safety, which will be suitable for assessing working conditions at any production, has been determined. In order to obtain a comprehensive index of labor safety, it is proposed to determine an estimate for each harmful production factor, and then to determine a single estimate, taking into account all characteristics. To determine the complex indicator of a harmful production factor, it is proposed to apply the integration method, that is, to find the area under the broken surface, which resulted from the connection of points on the plane of the XOY coordinate system. For this, quadrature formulas using the method of average rectangles are used. A step-by-step algorithm for determining the comprehensive indicator of labor safety, using integration by the method of average rectangles, is proposed. The method of determining the comprehensive indicator of labor safety in production was tested. Harmful factors in production were determined, valid indicators were obtained and their estimates on a dimensionless scale were determined. A time series of changes in indicators of harmful factors over time is graphically constructed. A comprehensive indicator of occupational safety in production was determined as a basis for making decisions regarding further actions to improve working conditions. Conclusions: It is suggested to use qualitative methods for comprehensive assessment of labor safety in production. Namely, to determine an assessment for each harmful and dangerous production factor and, then, to determine a single comprehensive assessment of working conditions. The method of determining the complex indicator of labor safety, using numerical integration by the method of average rectangles, is proposed. The proposed method can be considered universal, as it can be used for any premises and enterprises.