We consider the solution of ordinary differential equations (ODEs) using boundary value methods. These methods require the solution of one or more unsymmetric, large and sparse linear systems. The GMRES method with the Strang-type block-circulant preconditioner is proposed for solving these linear systems. We show that if an A-stable boundary value method is used for an m-by-m system of ODEs, then our preconditioners are invertible and all the eigenvalues of the preconditioned systems are 1 except for at most 2m(k + k) outliers. It follows that when the GMRES method is applied to solving the preconditioned systems, the method will converge in at most 2m(k + k) + 1 iterations. Numerical results are given to illustrate the effectiveness of our methods.