I think you'd also want your matrix to be positive semi-definite. Unless you're ok with complex numbers in your answer. What you'd normally do then is diagonalize your matrix thus: , where P is (obviously) invertible, and D is diagonal. Finally, you'd take the square root thus: However, this may not be the most general method available.