added adaptive runge-kutta order 5/6. Try to understand this method. Just as accurate as CVode but much slower. However, using this integration scheme, we can handle spike-mediated synaptic events. CVode needs to know the times of these discontinuities ahead of time.
Merge branch 'master' of https://github.com/dmf6/neuralsim2