powerset_zones function

Creates a set of all non-empty subsets of the integers from 1 to $n$.