library(MASS)
A <- matrix(c(1, 2, 3, 4, 5, 6, 7, 8, 9), nrow = 3, ncol = 3)
B <- matrix(c(5, 6, 7), nrow = 3, ncol = 1)
print (A,B)
## [,1] [,2] [,3]
## [1,] 1 4 7
## [2,] 2 5 8
## [3,] 3 6 9
A_inv <- ginv(A)
print(A_inv)
## [,1] [,2] [,3]
## [1,] -0.6388889 -5.555556e-02 0.5277778
## [2,] -0.1666667 -5.551115e-17 0.1666667
## [3,] 0.3055556 5.555556e-02 -0.1944444
X <- A_inv %*% B
print(X)
## [,1]
## [1,] 0.1666667
## [2,] 0.3333333
## [3,] 0.5000000
Penjelasan: