Implement LSODA

[hgrecco]

Adams/BDF method with automatic stiffness detection and switching.

See Scipy implementation

[Nicholaswogan]

I made a jit-able wrapper to LSODA: https://github.com/Nicholaswogan/NumbaLSODA/blob/main/README.md

Might be useful for implementing here. Does require a C++ compiler tho.