A numerical procedure for determining the hydromagnetic stability of axisymmetric systems is described. The growth rates and eigenfunctions are determined by means of a modified initial value method. A number of tokamak configurations have been studied and stability boundaries and growth rates are given. With a conducting wall on the plasma the fundamental toroidal mode is stable and the m = 1 mode is stable for q ⪆ 1/2.