I should add some day a trivial error checking, etc...
Also I had to write a series of functions in order to solve a large
number of small linear systems. For instance, say 10000 systems of
3x3. This arises in finite element analysis in order to avoid loops over
elements or gauss points.