Background. Mantle density models are key tools for understanding fundamental geological and physical processes occurring within the Earth. Many parameters used in mantle density models remain poorly understood and undefined. Among others, these include data on the composition and rheology of the mantle, which can vary significantly. Methods. The method of creating density models (density) significantly influences the final result. Modeling with one-dimensional models simplifies the calculation process but generalizes the distribution of mantle density, assuming it is homogeneous in the horizontal direction. This limitation does not allow for the consideration of lateral variations in mantle density, which can be important at the regional level. In this study, we present a quasi-three-dimensional model of mantle density beneath the Ukrainian Shield, obtained on the basis of a set of one-dimensional density curves. Polynomial corrections for heterogeneity were applied during the calculations, compensating for the shortcomings of one-dimensional models. This three-dimensional model was derived by recalculating one-dimensional velocity curves obtained by seismic tomography for 21 mantle domains in the depth range of 50 to 2600 km. The process of transforming P-wave velocity curves into a density model includes the following steps: determining seismic boundaries in the mantle as points of inflection of the first derivative of P-wave velocity curves for each mantle domain; creating a synthetic S-wave mantle model beneath the Ukrainian Shield by recalculating P-wave velocity curves; solving the Adams-Williamson equation using seismic velocities (P, S) for each domain with subsequent polynomial correction to account for heterogeneity; selecting a reference mantle model that would serve as the basis for converting velocity curves into density through the comparison of gravitational potential on the Earth's surface and calculated values from existing reference mantle models (PREM, PREMA, PREMC, IASP91 AK135). The AK135 model was chosen as the reference model based on the comparison of calculated and observed gravitational potential at the central point of the Ukrainian Shield. This study focuses on the final stages of constructing the mantle density model, taking into account mass balancing of the upper and lower mantle for each domain when determining density using the Adams-Williamson equation and introducing polynomial corrections relative to the AK135 reference model; calculating densities for each of the 21 mantle domains and their three-dimensional integration. Results. In this study, we present a quasi-three-dimensional model of mantle density beneath the Ukrainian Shield, obtained on the basis of a set of one-dimensional density curves, with polynomial corrections for inhomogeneity incorporated into the calculations, compensating for the shortcomings of one-dimensional model calculations. This three-dimensional model was obtained by recalculating one-dimensional velocity curves obtained by the seismic tomography method for P-waves, calculated for 21 mantle domains in the depth range from 50 to 2600 km. Conclusions. This study focuses on the final stages of constructing the mantle density model, considering balancing the mass of the upper and lower mantle for each domain in determining density using the Adams-Williamson equation and introducing polynomial corrections relative to the AK135 reference model; calculating densities for each of the 21 mantle domains and their three-dimensional integration. The obtained mantle density model of the Ukrainian Shield is well aligned with the division of the mantle into three main layers: lithosphere, upper mantle, and lower mantle. Each of the structural layers has its own visual pattern of density heterogeneity.