An iterative method for computing square roots of a class of matrices