# Modeling acceleration due to gravity for large bodies as a function of time

 P: 343 You could use the differential form of Gauss' Law of Gravity: $$\nabla \cdot g = -4 \pi G \rho$$ You would have to specify the mass density at each point on the masses, but then you can obtain the gravitational vector field for each mass. That might help you get started.