Function is_closed_under 

pub fn is_closed_under<T, F>(elements: &[T], operation: F) -> bool
where
    T: Copy + PartialEq,
    F: Fn(T, T) -> T + Copy,

Returns whether applying operation to any pair of values in elements yields another member of elements.