Application of the Parker-Sochacki numerical integration method to neuron models