Hi there guys, Currently writing and comparing two separate Mathematica scripts which can be found here and also here. The first one I've slightly modified to suit my needs and the second one is meant to reproduce the same results. Both scripts are attempting to simulate the trajectory of a timelike geodesic using the geodesic equations of motion. The Hamiltonian formalism for equations of motion. Now, the results should be equivelent. However, my problem is as follows: The Hamiltonian produces a timelike geodesic that falls radially inward directly to the event horizon, where the geodesic equation approach produces a beautiful precessing orbit. Now, the approach using the geodesic equation is correct, however, can anyone see why this is different in the second link above? I have no idea why they are not producing the same results!!