Benchmark suite

[gabrielgellner]

It would be great to have a set of benchmarks that output figures like in Hairer and Wanner showing the tolerance vs the number of function calls, also we could add tolerance vs cpu time.

It would also be good to see this for increasing the size of the problem using a general code for linear systems. That way we could see how the solvers perform for growing problem size.

[gabrielgellner]

As this package is doing reimplementation of ODE solvers in pure Julia it is essential that it has a good test suite. Thing to look at:

• The DETEST paper and its examples
• The Julia package IVPTestSuite.jl I will want to have this level of functionality at a minimum. Some other ideas is to have a general sweet of linear problems where we can compare the analytic to the approximate solutions. As well as a set of problems that we can compare to solvers from things like Matlab and Mathematica which are known good solvers.