The multistep generalized differential transform method is applied to solve the fractional-order multiple chaotic FitzHugh-Nagumo (FHN) neurons model. The algorithm is illustrated by studying the dynamics of three coupled chaotic FHN neurons equations with different gap junctions under external electrical stimulation. The fractional derivatives are described in the Caputo sense. Furthermore, we present figurative comparisons between the proposed scheme and the classical fourth-order Runge-Kutta method to demonstrate the accuracy and applicability of this method. The graphical results reveal that only few terms are required to deduce the approximate solutions which are found to be accurate and efficient
