Mir ist noch nicht klar, was genau du machen willst. Keplers Gesetze modellieren ja ein
ideales Universum. Dafür bräuchtest du nur die exakten Daten der Ellipsen und könntest über die Gesetze die jeweilige Position der Planeten auf diesen zu einem gegebenen Zeitpunkt berechnen.
Oder willst du eine
reale Körpersimulation bauen? Dafür brauchst du dann wirklich Beschleunigung, Geschwindigkeit, Position. Zentrifugalkraft allerdings nicht, das ist eine Scheinkraft
.
In diesem Fall könnten - im Gegensatz zu Keplers statischem Modell - sich allerdings Ungenauigkeiten hochschaukeln, was neben dem Zeit-Delta vor allem von den Startparametern abhängt: Anfangsposition und -geschwindigkeit jedes einzelnen Körpers.
Genauer gehst du so vor: Jeder Körper besitzt eine Position p und einen Geschwindigkeitsvektor v. Wie du schon beschrieben hast, summierst du die Kräfte zu allen anderen Körpern und kommst so auf einen Beschleunigungsvektor a.
Die neue Geschwindigkeit ist dann v' = v + dt * a und die Position p' = p + dt * v.