[Haskell-cafe] matrix computations based on the GSL