Данная работа представляет собой исследование движения в гравитационном поле с центральным пиком плотности. Рассмотрена сферическая модель гравитационного потенциала и проведено ее обобщение на трехосный случай с помощью метода эквипотенциалей. Получены ограничения на структурные параметры из условий физической корректности модели и проанализировано влияние параметров на вид орбит. Задача Коши для системы дифференциальных уравнений решается с помощью различных численные методов. Система является гамильтоновой, поэтому возможно применение симплектических методов Эйлера и Нюстрёма второго порядка. Для более точных вычислений были также использованы вложенные методы Мерсона и Дормана-Принса. По результатам численного интегрирования было построено множество орбит.