MASS for matrix operation

library(MASS)

mat <- matrix(1:9, nrow = 3, ncol = 3)
print("Original Matrix:")

## [1] "Original Matrix:"

print(mat)

##      [,1] [,2] [,3]
## [1,]    1    4    7
## [2,]    2    5    8
## [3,]    3    6    9

inverse <- ginv(mat)
print("Inverse Matrix:")

## [1] "Inverse Matrix:"

print(inverse)

##            [,1]          [,2]       [,3]
## [1,] -0.6388889 -5.555556e-02  0.5277778
## [2,] -0.1666667 -1.387779e-17  0.1666667
## [3,]  0.3055556  5.555556e-02 -0.1944444

determinant <- det(mat)
print("Determinant of the Matrix:")

## [1] "Determinant of the Matrix:"

print(determinant)

## [1] 0

eigen_result <- eigen(mat)
eigenvalues <- eigen_result$values
eigenvectors <- eigen_result$vectors
print("Eigenvalues:")

## [1] "Eigenvalues:"

print(eigenvalues)

## [1]  1.611684e+01 -1.116844e+00 -5.700691e-16

print("Eigenvectors:")

## [1] "Eigenvectors:"

print(eigenvectors)

##            [,1]       [,2]       [,3]
## [1,] -0.4645473 -0.8829060  0.4082483
## [2,] -0.5707955 -0.2395204 -0.8164966
## [3,] -0.6770438  0.4038651  0.4082483