[FEATURE] Port ODE solvers from SciPy to SciJo

[shivasankarka]

What & Why

This issue is related to porting the ODE solvers available in the scipy quad module to the SciJo library. The goal is to enable SciJo to support the same ODE solving capabilities as found in scipy, including integration routines and solver interfaces.

Proposed Solution

Tasks:

• Identify all ODE solvers and related functions in scipy's quad module.
• Implement equivalent ODE solvers in SciJo.
• Ensure API compatibility with scipy's conventions where possible. API compatibility is welcome, but not a necessity if it means more performance.
• Add documentation and examples for the new solvers.
• Write unit tests to verify correctness and performance against scipy.

Extra Context

References:

• scipy.integrate.quad documentation
• SciJo repository

If anyone has experience with scipy ODE solvers or suggestions for implementation strategies, please comment below.