SOLVEBLOK: A Package for Solving Almost Block Diagonal Linear Systems

In spline approximation and m the numerical solution of (systems of) orchnary differential equations, one IS eventually faced with linear systems of a rather characteristic, almost block diagonal, form which reflects the underlying geometry of the problems. Whde these systems could be treated as banded systems, it is more efficient, both m the generation of the systems and m their numerical solution, to work explicitly with their block structure.