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