构造初等矩阵

fabric_Elementary_matrix <- function(i =0,
                                     j = 0,k = 1,n)
{
 # n代表矩阵的维数
  A <- diag(1, n)
  if(i!= 0&j == 0)#i行乘k倍
  {
    A[i,] <- k*A[i,]
  }
   if(i!= 0&j != 0)
  {
     if(k!=1)#i行的k被加到j行
     {
       A[j,] <- A[j,] + k*A[i,]
     }
     else#两行互换
     {
       temp1 <- A[i,]
       A[i,] <- A[j,]
       A[j] <- temp1
     }
   }
  return(A)
}