  Numerical calculations have been made on the spike-train response of a pair
of Hodgkin-Huxley (HH) neurons coupled by synapses and axons with time delay.
The recurrent excitatory-excitatory, inhibitory-inhibitory,
excitatory-inhibitory, and inhibitory-excitatory couplings are adopted. The
coupled, excitable HH neurons are assumed to receive the two kinds of
spike-train inputs: the transient input consisting of M impulses for the finite
duration (M: integer) and the sequential input with the constant interspike
interval (ISI). The HH neurons in all the kinds of couplings are found to play
a role of memory storage with on-off switching. When the coupling strength and
the time delay are changed, the distribution of the output ISI $T_{\rm o}$
shows bifurcation (multifurcation), metastability and chaotic behavior. The
calculation of the time correlation shows that the synchronization between the
two HH neurons is well preserved even when the distribution of their $T_{\rm
o}$ is chaotic. The correlation dimension of the cycles of $T_{\rm o}$ is shown
to depend not only on the model parameters but also on the type of input ISIs.
