To provide Jacobian in non-linear solver, time-stepping or Jacobian

It won’t if automatic differentiation is around because that will generally be accurate enough. It’s not hard to show against finite differences though, take a very stiff equation and do something like Rodas5(autodiff=false) (ROBER?) and you’ll see the difference (since Rodas methods require an accurate Jacobian)