Pada program ini menggunakan library(matlib). Perhitungan matriks kali ini dibantu dengan library tersebut. 

library(matlib)
## Warning: package 'matlib' was built under R version 4.2.2

Pada kode di bawah ini adalah contoh dari pengkodean matriks dengan baris 3 dan kolom 3. Dan untuk variabel b yaitu hasil dari setiap perhitungan baris dari matriks. 

A <- matrix(c(0, -2, -1, 2, 3, 2, 3, -2, 3), nrow = 3, ncol = 3)
b <- c(6, -1, 2)
A
##      [,1] [,2] [,3]
## [1,]    0    2    3
## [2,]   -2    3   -2
## [3,]   -1    2    3
b
## [1]  6 -1  2

Jika ingin megetahui nilai x dari setiap variabel x, kita dapat menghitungnya dengan Solve(A, b). Ketika di-run, maka akan ada nilai x yang dimana nilai tersebut untuk menyelesaikan persamaan linear. 

Solve(A, b)
## x1      =           4 
##   x2    =  2.53846154 
##     x3  =  0.30769231

Berikut ini adalah contoh lainnya: 

A <- matrix(c(1,0,1,0,1,0,0,1,0,1,1,0,0,1,0,1), nrow = 4, ncol = 4)
b <- c(475, 489, 542,422)
A
##      [,1] [,2] [,3] [,4]
## [1,]    1    1    0    0
## [2,]    0    0    1    1
## [3,]    1    0    1    0
## [4,]    0    1    0    1
b
## [1] 475 489 542 422
Solve(A, b)
## x1     - 1*x4  =   53 
##   x2     + x4  =  422 
##     x3   + x4  =  489 
##             0  =    0

Sumber Referensi: Ruriko Yoshida - Linear Algebra and its Application with R