A<-matrix(sample(20:40,4),ncol=2)
matrix.sqrt(A)
A<-matrix(sample(20:40,16),ncol=4)
matrix.sqrt(A)
#sqrt of inverse of A, or sqrt inverse of A
matrix.sqrt(solve(A))
#non-square matrix
A<-matrix(sample(20:40,20),ncol=4)
matrix.sqrt(A)
Run the code above in your browser using DataLab