We propose a numerical method for calculating bifurcations of periodic solu
tions observed in a model equation of Hodgkin-Huxley neurons coupled by exc
itatory synapses with a time delay. To illustrate the validity of the metho
d, bifurcations in two-coupled Hodgkin-Huxley equations with Variation of a
coupling coefficient and time delay are studied.