Add support for finding the power set

This patch adds powerset() method to sympy.core.sets:Set to find the
power set of a set.