Abstract: The fundamental theory of interpolating matrix method is eatablished for solving mixed order systems of nonlinear multipoint boundary value problems of ODEs.A general purpose solver IVMMS of ODEs is written based on this know ledge,and which can support the finite element method of lines for solid mechanics.