sqrtmat function

Square root of nonsymmetric matrix