It's so quiet on the forum today, I might as well make some noice then ;) Since last time [BarrOff](https://github.com/BarrOff) has made several contributions to the ode solvers both by adding higher-order integrators like Tsit54, Vern65. The adaptive ode solver does now also use both an absolute and relative error tolerence instead of just an absolute tolerence it used before.
A new kind of Spline type has been introduced: HermiteSpline. It supports generic types in contrast to the natural spline which only support floats currently. If HermiteSpline is provided with function values and it's derivative it is third-order accurate. By suggestion of @Vindaar the adaptiveGauss integral function supports infinite integration intervals and has a much more robust global error control instead of the local error control used before. Reaching machine epsilon accuracy was impossible before but now it has been possible for "nice" integrands.
