Test for unitary matrix
TF = isunitary(X, tol)
isunitary(X, tol) returns logical 1 (true) if the square matrix X is unitary (with transpose equal to its inverse) to within the tolerance given by the second parameter, if present. If no second parameter is specified, the tolerance defaults to eps.