p_union function

Set operations