M <- matrix(c(2, 5, 6, 1, 8, 10, 9, 4, 7), nrow = 3, byrow = TRUE)

M1 <- cbind(M, c(3, 4, 1))
M2 <- cbind(c(7, 4, 5), M)
sum_matrix <- M1 + M2
prod_matrix <- M1 * M2

M
M1
M2
sum_matrix
prod_matrix