a <- matrix(1:4, nrow=2)
b <- c(1, 2)
r <- solve(a,b)
a %*% as.matrix(r)
a_inv <- solve(a)
a_inv %*% a

####################################################

det(a)

####################################################

cor(1:10,2:11)