Hello there, I've been considering the geodesic equations of motion for a test particle in Schwarzschild geometry for some time now. Similar to what we can do with the Kepler problem I would like to be able to numerically integrate the equations of motion. I'm quite interested to see how coordinate time and proper time change over a simulation and also the trajectory itself. But the problem is.. I have no idea how to do that. Any suggestions?