Before revision:  Thus, we adaptively assign the ETD2RK and RK2 schemes depending on the stiff and flat regimes to achieve optimal accuracy and efficiency . Compared with the regular RK2 method , our ETD2RK method can enlarge time steps and raise efficiency more than one order of magnitude while capturing excellent accurate traces of membrane potentials and spike times of HH neurons.
After revision:  Thus, we adaptively assign the ETD2RK and RK2 schemes depending on the stiff and flat regimes to achieve optimal accuracy and efficiency . Compared with the regular second-order Runge-Kutta method (RK2), our AETD2 method can use time steps one order of magnitude while capturing excellent accurate traces of membrane potentials and spike times of HH neurons.