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 RK2 method , our ETD2RK method can enlarge time steps and raise efficiency more than one order of magnitude larger and improve computational efficiency more than ten times while excellently capturing accurate traces of membrane potentials and spike times of HH neurons.