a <- matrix(c(1, 2, 3, 4), 2, 2)
a^1
a^2