We were using row-vector multiplication for powers (taken from the
superclass), but our vectors are column vectors. Oops. This broke
things when we assumed column vectors were being used, like when
we constructed a multiplication table. This commit fixes the powers
and adds/updates some tests.