Matrix Operations :
> sng<-matrix(1:4,2,2,byrow=F)
> sng
[,1] [,2]
[1,] 1 3
[2,] 2 4
> t(sng)
[,1] [,2]
[1,] 1 2
[2,] 3 4
> solve(sng)
[,1] [,2]
[1,] -2 1.5
[2,] 1 -0.5
> solve(t(sng))
[,1] [,2]
[1,] -2.0 1.0
[2,] 1.5 -0.5
> A<-sng
> A%*%A #A(prime).A
[,1] [,2]
[1,] 7 15
[2,] 10 22
> A%o%A #A.A(prime)
, , 1, 1
[,1] [,2]
[1,] 1 3
[2,] 2 4
, , 2, 1
[,1] [,2]
[1,] 2 6
[2,] 4 8
, , 1, 2
[,1] [,2]
[1,] 3 9
[2,] 6 12
, , 2, 2
[,1] [,2]
[1,] 4 12
[2,] 8 16