Find all Solutions to this linear system:
x + y - z = -5 x - y - z = -3 x + y + z = 0
M <- matrix(c(1,1,-1,-5,1,-1,-1,-3,1,1,1,0),nrow = 3, byrow = TRUE)
M
## [,1] [,2] [,3] [,4]
## [1,] 1 1 -1 -5
## [2,] 1 -1 -1 -3
## [3,] 1 1 1 0
M[3,] <- M[2,]+ M[3,]
M
## [,1] [,2] [,3] [,4]
## [1,] 1 1 -1 -5
## [2,] 1 -1 -1 -3
## [3,] 2 0 0 -3
M[2,] <- M[2,]-M[1,]
M
## [,1] [,2] [,3] [,4]
## [1,] 1 1 -1 -5
## [2,] 0 -2 0 2
## [3,] 2 0 0 -3
X <- M[3,4]/M[3,1]
X
## [1] -1.5
Y <- M[2,4]/M[2,2]
Y
## [1] -1
Z <- (M[1,4]-M[1,2]*Y - M[1,1]*X)/M[1,3]
Z
## [1] 2.5