A <- matrix(c(2, 1, 1, -1), nrow=2, ncol=2, byrow=TRUE)
b <- c(5, 1)
Ab <- cbind(A, b)
Ab
##           b
## [1,] 2  1 5
## [2,] 1 -1 1
m <- nrow(Ab)
n <- ncol(Ab)
print(m,n)
## [1] 2
cat("Matriks koefisien A:\n")
## Matriks koefisien A:
print(A)
##      [,1] [,2]
## [1,]    2    1
## [2,]    1   -1
cat("\nVektor konstanta b:\n")
## 
## Vektor konstanta b:
print(b)
## [1] 5 1
x <- solve(A, b)
cat("\nSolusi dari sistem persamaan linear Ax=b adalah:\n")
## 
## Solusi dari sistem persamaan linear Ax=b adalah:
print(x)
## [1] 2 1

Pada contoh kode di atas, kita menginput matriks koefisien A dan vektor hasil b, kemudian kita menggabungkan keduanya menjadi matriks augmented Ab.