Третья вертикальная производная аномалии силы тяжести находит применение в гравиразведке для оценки параметров аномальных тел, в физической геодезии при расчете третьего поправочного члена в рамках решения задачи М. С. Молоденского методом аналитического продолжения. Основная цель исследования заключалась в разработке и анализе алгоритмов вычисления данной производной двумя способами: с использованием численного дифференцирования на основе пятиточечной конечноразностной схемы, когда аномалии заданы вдоль вертикальной оси Oz относительно расчетной точки с постоянным шагом по высоте, и аналитическим методом, основанным на разложении в усеченный ряд Фурье по сферическим гармоникам с применением коэффициентов глобальной модели EIGEN-6C4. По итогам реализации созданных алгоритмов и программного обеспечения построены две глобальные цифровые модели третьей вертикальной производной аномалии силы тяжести. Анализ показал, что при решении задачи М. С. Молоденского с использованием аналитического приближения третьего порядка вычисление этой производной предпочтительнее осуществлять аналитическим путем – через разложение в усеченный ряд сферических функций, что обеспечивает большую экономическую эффективность по сравнению с численными методами