The poliastro.threebody contains a set of modules related to this physical problem:
poliastro.threebody