library(matlib)
A <- matrix(c(1,2,1,1), nrow=2, ncol=2)
A
## [,1] [,2]
## [1,] 1 1
## [2,] 2 1
b <- c(4,5)
b
## [1] 4 5
Solve(A,b)
## x1 = 1
## x2 = 3
plotEqn(A,b)
## x[1] + x[2] = 4
## 2*x[1] + x[2] = 5