Статья посвящена методике вычисления изостатических аномалий силы тяжести на основе гипотезы региональной компенсации. В качестве исходных данных используются аномалии Буге. Предполагаемая мощность H земной коры вычисляется по средневзвешенным высотам рельефа; за весовую функцию принят направляющий косинус давления масс рельефа на плоскость z=H0. Интегралы, входящие в расчетные формулы, вычисляются методом Люстерника-Диткина. Приводятся численные значения коэффициентов кубатурных формул и операторная схема программы, составленной для ЭВМ БЭСМ-6