He estado jugando con algunos simuladores orbitales -muy sencillos-, sobre todo utilizando librerías de física preexistentes (también intenté hacerlo con vectores). Lo que me confunde es que mis órbitas no se comportan como en la realidad: el primario siempre está en el centro de la elipse y no en uno de los focos. Obtengo el mismo resultado independientemente del motor o la biblioteca que utilice.
Me he limitado a poner un primario y un orbitador, con el primario en el centro de la distribución. Introduzco la fórmula $$F = G\frac{m_1m_2}{r_1^2}$$ con la fuerza dirigida hacia el primario. He intentado ajustar el paso de tiempo, pero obtengo el mismo resultado. Simplemente estoy confundido en cuanto a lo que podría causar esto.
Actualización
Tienes razón, estaba multiplicando, sólo que no lo sabía (algo sobre la forma en que el Construct2 interpretó mis comandos). Lo tenía configurado para aplicar la fuerza que era esencialmente G(arbitrario) / distancia * distancia (no me gustó cuando probé ^2). Cuando puse la distancia al cuadrado en una variable y luego dividido por ella las cosas funcionaron.
Ahora ya había probado esto en Panda3d con el mismo problema, así que tendré que volver a mirar eso.