C <- matrix(c(1,0,1,1,2,2,-1,1,2,1,3,0,1,1,0,1), nrow=4, byrow=T) C
## [,1] [,2] [,3] [,4] ## [1,] 1 0 1 1 ## [2,] 2 2 -1 1 ## [3,] 2 1 3 0 ## [4,] 1 1 0 1
det(C)
## [1] 4