A <- matrix(c(3,9,1,4), nrow=2, ncol=2, byrow = TRUE)
b <- c(6,1)
Ab <- cbind(A,b)
Ab
##          b
## [1,] 3 9 6
## [2,] 1 4 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,]    3    9
## [2,]    1    4
cat("\nVektor konstanta b:\n")
## 
## Vektor konstanta b:
print(b)
## [1] 6 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]  5 -1

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