Предлагается новый способ вычисления вторых производных гравитационного потенциала для трехмерных тел произвольной формы и переменной плотности масс по гармоническим моментам. При вычислении учитывается необходимое число членов отрезка ряда для достижения заданной точности с учетом местоположения точки относительно аномальных масс и размеров последних. Приводятся явные выражения для всех вторых производных гравитационного потенциала, а также дается оценка погрешности предлагаемого метода.