I have differential equations of geodesics. I need to solve them numerically and the solution must be on a sphere.This is what the equations looks like
The geodetic on the sphere should look like this depending on the parameters
How do I numerically solve the geodesic differential equations and then plot the solution on the sphere, which would turn out as in the image above?




