Block-boundary value methods for the solution of ordinary differential equations