Implementing Runge–Kutta–Fehlberg.
Implementing Runge–Kutta–Fehlberg, part 2.
Basing error estimate for RKF45 on new instead of old state.
Optimising Runge–Kutta algorithms a little.
Further Optimizations.
Introducing NOGSL flag.
Making fractions in Runge–Kutta methods baseType.
changing it++ to ++it for performance