我们由外向内依次将半径在r~r+dr内的球壳膨胀到无穷远处,计算需要的总能量就好了. 　　如果已经计算过r+dr~R部分的能量,那么意味着这部分已经飘向无穷远.现在只需要将r~r+dr与0~r两部分分离即可算出这部分的能量. 　　设0~r部分质量为m1(r),r+dr部分在方位角O~O+dO范围内的质量为m2(r,dr,O,dO),显然两者质心相距r,现在将m2移动到距球心D处,万有引力为：F=-G*m1*m2/D（G为万有引力常数）.从r移动到“无穷远”万有引力做负功：W(r,m1,m2)=∫[D,无穷]FdD=-G*m1*m2/r. 　　将m1=M*(r/R)^3,m2=[(3Mr^2)/(4*pi*R^3)]drdO带入得到W(r,dr,dO)=[-(3GM^2r^4)/(4*pi*R^6)]drdO. 　　将球壳一层层全部膨胀到无穷远处万有引力做功之和（负功）为：W=∫[0,R]drdO[-(3GM^2r^4)/(4*pi*R^6)]=[-(3G*M^2)/(5R)]. 　　所以克服万有引力需要的能量为（正）：E=+(3G*M^2)/(5R),这就是结合能. 　　特别地,如果这是地球,那么地表的重力加速度g=GM/R^2,此时E=3gMR/5（仅当此球为地球的情况下成立）.