Numerical calculation of stable and unstable periodic solutions to the Hodgkin-Huxley equations