x <- matrix(c(4,5,6,7), nrow=2, ncol=2, byrow = TRUE)
y <- c(1,2)
xy <- cbind(x,y)
xy
##          y
## [1,] 4 5 1
## [2,] 6 7 2
m <- nrow(xy)
n<- ncol(xy)
print(x,y)
##      [,1] [,2]
## [1,]    4    5
## [2,]    6    7
cat("Matriks Koefisien x:\n")
## Matriks Koefisien x:
print(x)
##      [,1] [,2]
## [1,]    4    5
## [2,]    6    7
cat("\nVektor konstanta y:\n")
## 
## Vektor konstanta y:
print(y)
## [1] 1 2
x<- solve(x,y)
cat("\nSolusi dari sistem persamaan linear Xx=y adalah:\n")
## 
## Solusi dari sistem persamaan linear Xx=y adalah:
print(x)
## [1]  1.5 -1.0

Pada contoh kode si atas, kita menginput matriks koefisien A dan vektor hasil b, kemudian kita menggabungkn keduanya menjadi matriks augmented Ab