Correction of an error in the matrix elements used by Salwen & Grosch (1972) has brought the results of the matrix-eigenvalue calculation of the linear stability of Hagen–Poiseuille flow into complete agreement with the numerical integration results of Lessen, Sadler & Liu (1968) for azimuthal index n = 1. The n = 0 results were unaffected by the error and the effect of the error for n > 1 is smaller than for n = 1. The new calculations confirm the conclusion that the flow is stable to infinitesimal disturbances.

Further calculations have led to the discovery of a degeneracy at Reynolds number R = 61·452 ± 0·003 and wavenumber α = 0·9874 ± 0·0001, where the second and third eigenmodes have equal complex wave speeds. The variation of wave speed for these two modes has been studied in the vicinity of the degeneracy and shows similarities to the behaviour near the degeneracies found by Cotton and Salwen (see Cotton 1977) for rotating Hagen-Poiseuille flow. Finally, new results are given for n = 10 and 30; the n = 1 results are extended to R = 106; and new results are presented for the variation of the wave speed with αR at high Reynolds number. The high-R results confirm both Burridge & Drazin's (1969) slow-mode approximation and more recent fast-mode results of Burridge.