可以按http://sobereva.com/131的方法。
注意轨道贡献的定义不是唯一的，所以计算方法也不唯一，以上链接只是可能的计算方法之一。如果要严格重复文献的做法，还是得看文献的计算细节部分的描述才行。