A general precise integration method for twopoint boundary value problems of ordinary differential equations is presented. Using the transfer matrix

interval algebraic equations are established and the recursive formulas for interval combination and condensation are derived for various boundary conditions. By utilizing the results of the transfer matrix directly

the process of interval combination and condensation has high efficiency. In addition

the present method has wide application for which can deal with various boundary value problems. Numerical examples show the validity of the present method.