I'm pretty new to SageMath and I'm currently trying to work on solving geodesic equations in well-defined spacetime manifolds. I found the documentation of SageMath pretty extensive and broad and it gave me a lot of confidence in the software. But I'm still hesitating to use Sage ODE solvers because I know the capabilities of RK45 from scipy.integrate as I have use that many times. I just want to know on a broad scale, are SageMath solvers faster than Scipy solvers and how reliable and stable are the solutions in general ?
Thanks in advance for the replies!
