Abstract. Numerical integration formulas of interpolatory type are generated by the integration of g-splines. These formulas, which are best in the sense of Sard, are used to construct predictor-corrector and block implicit schemes. The schemes are then compared with Adams-Bashforth-Adams-Moulton and Rosser schemes for a particular set of prototype problems. Moreover, an improved error bound for linear multistep formulas based on g-splines and a comparison of L2 norms of Peano kernels for Adams and natural spline formulas are given.