Rodas4 (crashes) versus Tsit5 (works)

[Oscar_Smith]

Implicit solvers use autdiff by default to compute the jacobian. For functions that don’t work with autodiff, you should use Rodas4(autodiff=false)