Abstract: Interpolating matrix method is a numerical method for solving multi-point boundary value problems. This paper provides the error analysis of the method, proving the same precision for the solutions y(x), y'(x), ..., y(^m)(x)obtained by the method. Author gives the convergence order estimate and proves good stability property of the method the second order scalar equations.