realsets.py is a [Sage] module to support subsets of the real field consisting of intervals and isolated points and was developed to demonstrate set operations of the MGL
It is based of previous work from Interval1Sage adding integration on real sets and real intervals.
An object in this module consists of a list of disjoint open intervals plus a list of isolated points (not belonging to these intervals). Notice that
Infinite is acceptable as interval bound.