I would like to write a program that solves the trajectories of objects (think rockets) that are influenced by mass of other objects (think planets). I saw that I can do this using Newton's laws, but I would like to confirm that another way is also possible.
If I could find a metric of a spacetime distorted by masses of some planets, then I could find Christoffel symbols and plug them into the geodesic equation to find ODEs describing the geodesics in this spacetime (for rockets). Than I could solve the regular N-body problem and recompute the metric. Is this right? If so, how can I build a metric of such a spacetime?
I've seen people find trajectories of photons in the Schwarzschild metric, and I consider it similar to my problem (except I want to have more than one object distorting the spacetime). I went through all questions with [differential-geometry] tag, but couldn't find anything related to my problem.
I'm very sorry about this question being not very specific. I know very little about GR and differential geometry, but I don't think it prevents me from using it as a tool.
Found a chapter "Multiple Black Hole Solutions" in http://www.aei.mpg.de/~gielen/black.pdf. It shows how to build a Reissner–Nordström metric describing a system of few black holes. I think, that this could be simplified into Schwarzschild metric by neglecting the charge of the black holes. If so, how can one obtain the metric for an object which is not a black hole?