Very nice piece of work. As a pointer to the extra documentation I would suggest to add the reference /Sch86a/ in full. A minor
detail: the variables (e.g. x, y, z) are eliminated left to right, leaving
constant bounds for the inner (last) variable, say min<= z <= max. I can
think of applications where it is interesting to obtain constant
bounds for the outer (first) variable, i.e. min<= x <= max. This can easily done by reordering the columns of matrix A.
