def symmetric_closure_function(rel):
    """
    Function to determine the symmetric closure of a relation
    :param rel: A list that represents a relation
    :return: The symmetric closure of relation rel
    """
    return rel + [(y, x) for (x, y) in rel if (y, x) not in rel]