Abstract We will introduce numerical methods for systems of ODEs by considering the celebrated FitzHugh-Nagumo model published FitzHugh [1] in 1961 and, independently, Nagumo et. al. [2] 1962. The is a system ordinary differential equations with two unknowns, and commonly used as simple action potentials excitable pacemaker cells.