mat_a <- a <- matrix(c(1,2,3,-2,0,1,-2,-1,1,1,3,1), nrow = 3, byrow = TRUE)
mat_a
## [,1] [,2] [,3] [,4]
## [1,] 1 2 3 -2
## [2,] 0 1 -2 -1
## [3,] 1 1 3 1
vec_b <- matrix(c(-7,3,1,1))
vec_b
## [,1]
## [1,] -7
## [2,] 3
## [3,] 1
## [4,] 1
mat_a %*% vec_b
## [,1]
## [1,] 0
## [2,] 0
## [3,] 0