Add integral power routine

A real power routine depends on support for complex matrices which we don't provide yet.
Given these fast paths are here, should we also have 2, which does m.Mul(a,a)? It's probably a common case.