Численное моделирование основывается на статическом подходе, при котором планета моделируется как упругое тело с зависящими от радиуса плотностью, модулем сжатия К и модулем сдвига. Считается, что деформации и напряжения, подчиняющиеся закону Гука, обусловлены давлением на поверхность планеты топографических структур и аномалиями плотности, распределенными некоторым образом в коре и мантии. Расчеты проводятся с помощью техники функций Грина (или методе нагрузочных чисел), развитой в работах [1,2,3,4,5]. С тех пор данные о топографии и гравитационном поле Марса существенно улучшены [6,7]. В работе топография и гравитационное поле Марса определяются по отношению к референсной поверхности, за которую выбирается равновесный сфероид. [8]. В качестве тестовой модели используется модель M7 из работы [8], которая удовлетворяет всем имеющимся на сегодняшний день геофизическим и геохимическим данным. В тестовой модели M7 средняя толщина коры - 50 км и средняя плотность коры равна 3000 кг м-3. Поскольку невозможно однозначно определить плотностные аномалии по данным о гравитационном поле планеты, то для определенности задачи делалось предположение, что в Марсе существуют два уровня сосредоточения аномалий – поверхность планеты и граница кора-мантия. Аномальный слой плотности действует на планету как нагрузка, для учета подстройки коры и мантии к этим нагрузкам вводится множитель (1+kn(r)): где kn(r) – нагрузочные числа порядка n для заглубленной на глубине r аномалии [1,3,9]. Для модели М7 рассчитаны нагрузочные числа kn(r), матрица напряжений, напряжения сжатия-растяжения, а также максимальные сдвиговые напряжения.