A Collocation Solver for Mixed Order Systems of Boundary Value Problems

Implementation of a spline collocation method for solving boundary value problems for mixed order systems of ordinary differential equations is discussed. The aspects of this method considered include error estimation, adaptive mesh selection, B-spline basis function evaluation, linear system solution and nonlinear problem solution. The resulting general purpose code, COLSYS, is tested on a number of examples to demonstrate its stability, efficiency and flexibility.