Functions that solve boundary value problems ('BVP') of systems of ordinary
differential equations ('ODE') and differential algebraic equations ('DAE').
The functions provide an interface to the FORTRAN functions
'twpbvpC', 'colnew/colsys', and an R-implementation of the shooting method.