Add integral power routine

A real power routine depends on support for complex matrices which we don't provide yet.
Could you add a comment about using sequential squaring? It looks very different from most for loops, and seeing a comment about sequential squaring will help future readers who aren't used to bit shifting.