IDL does not have a specific routine dedicated to find the square root(s) of a matrix. You would need to write the code for that yourself using the theory behind the calculation.

However, I do note that light online searching turned up a few existing third -artly IDL routines that would do the trick!