Fix reference to matrix.alloc in matrix implementation
Implement matrix.zero method
Use native GSL functions to set matrix to zero or identity
Implement matrix powers with integer exponents