[adgjlsfhk1]

You might want to try OrdinaryDiffEq.jl. As of last summer, we have really fancy radau methods that are arbitrary order, arbitrary tolerance (see https://arxiv.org/abs/2412.14362 for the paper)