* add very simple test for QR decomposition
* add function matrix_absmax that returns the larges absolute value of the
  matrix
* add function matrix_absmin that returns the smallest absolute value of the
  matrix
* add function matrix_balance that balances a matrix
* add function matrix_froebenius that calculates the froebenius norm
* fix sign error