** C31_Discussion_1_Daniel_Thonn **

Format of equation: ax + b y + c*z = d

Prepare Matrix:

M1 <- matrix(c(1,1,0,1,-1,0,-1,-1,1, -1,-1,2),3,4)
M1
##      [,1] [,2] [,3] [,4]
## [1,]    1    1   -1   -1
## [2,]    1   -1   -1   -1
## [3,]    0    0    1    2

so: z=2 ;

Subtract row 1 from row 2

M1[2,] <- M1[2,] -M1[1,] 
M1
##      [,1] [,2] [,3] [,4]
## [1,]    1    1   -1   -1
## [2,]    0   -2    0    0
## [3,]    0    0    1    2

Therefore: z = 2 ; y = 0 ;

and: x + y -z = -1 ; x + 0 -2 = -1 ;

x = 1

Conclusion: z = 2 ; y = 0 ; x = 1 ;

END