Added functions for pseudo-inverse.

I added the pinv and pinvTol functions.
The types are essentially the same as for inv, except without the requirement that the input matrix be square.
The pinvTol is conditionally compiled because the implementation in hmatrix 13 that one of my computers is stuck on doesn't have a sufficiently general type signature.