Я же сказал: Эта формула всем хороша, кроме того, что она записана не через ту

. В эту величину входит масса также и самого гравитационного поля. А нам нужны только камни (составляющие сферу).
Честно говоря, уже перестал понимать как решение, так и , главное, условие задачи.
Хотелось бы еще раз понять, что вычисляется. Есть тонкостенная сфера , массой

(для определенности радиуса

).

- интеграл движения, входящий в шварцшильдовское решение. Можно записать такое выражение:


- это масса вещества внутри оболочки (которая входит в ваши вычисления),

- энергия гравитационного поля внутри оболочки (где вещество), третий член соответственно энергия поля вне сферы.
Далее, я так понял, вы берете слой

из данной оболочки и сферически симметрично поднимаете на величину

. Теперь у нас 4 области:
![$Mc^2=[\bar{W}_m(r=a)+\bar{W}_g(r=a)]+\bar{W}_g(a<r<a+\Delta{l})+[\bar{W}_m(r=a)+\bar{W}_g(r=a)]+\bar{W}_g(r>a+\Delta{l})\quad(b)$ $Mc^2=[\bar{W}_m(r=a)+\bar{W}_g(r=a)]+\bar{W}_g(a<r<a+\Delta{l})+[\bar{W}_m(r=a)+\bar{W}_g(r=a)]+\bar{W}_g(r>a+\Delta{l})\quad(b)$](https://dxdy-01.korotkov.co.uk/f/8/6/3/863bccec483c43b835dd3fc35e91ed7a82.png)
Первая скобка - полная энергия внутри первой оболочки, второй член - энергия поля между оболочками, третий - полная энергия внутри второй оболочки, третий - энергия поля вне вещества до бесконечности.
Насколько я понял Вы хотите определить полную энергию поля между оболочками (?). Как определить перераспределение энергии поля между участками при данной процедуре не очень понятно.