Hello, I'm currently studying general relatively and am trying to plot orbits of planets around the sun using the schwarzchild metric. I've worked out the geodesic equations, working with c=1 to simplify things, and written a MATLAB script to plot trajectories, but I'm struggling to work out...