# Computing the matrix-square-root of a positive definite square matrix:
## Not run:
# # generate random data
# rand_data(500,5000)
#
# W <- ccross(M)
#
# # this is the implementation of the matrix power-operator '%**%'
# W_sqrt <- with(eigen(W), ccross(vectors,values**0.5))
# ## End(Not run)
Run the code above in your browser using DataLab