I'm working on a simulation program that replicates the movement of planets around a large celestial body (the sun). This is a three dimensional simulation that uses vectors. At present, I'm struggling to work out how to calculate the correct initial velocity vector so that the planet will enter into a stable orbit around the sun.
I've attempted to take the vector cross product of
v0 but to no avail. I also attempted to calculate..
$T_0 = v_1 - v_0$
$T_1 = v_0 - v_1$
$NV = T_0 * T_1$
NV is the new velocity and
* is the cross product function. This resulted in a zero vector because $T_0$ is parallel to $T_1$ so I find myself somewhat stuck. Does anyone know how to compute the correct vector?