Построен эффективный алгоритм решения обратной задачи гравиметрии о восстановлении переменной плотности в горизонтальном слое по гравитационным данным. После дискретизации и аппроксимации задача сводится к решению системы линейных алгебраических уравнений. Идея алгоритма основана на использовании блочно-тёплицевой структуры матрицы коэффициентов. Использование данного алгоритма значительно снижает расходы памяти и уменьшает время решения задачи. Проведено распараллеливание алгоритма и его реализация на суперкомпьютере "Уран". Решена модельная задача о восстановлении плотности по синтетическим гравитационным данным. An effective algorithm for an inverse gravity problem of finding a variable density in a horizontal layer using gravitational data is constructed. After the discretization and approximation, the problem reduces to solving a system of linear algebraic equations. The idea of this algorithm is based on exploiting block-Toeplitz structure of coefficients matrix. Utilizing of this algorithm drastically reduces memory usage and also reduces computation time. The algorithm was parallelized and implemented using the Uran supercomputer. A model problem with synthetic gravitational data was solved.