Symbolic powers is a classical commutative algebra topic that relates to primary
decomposition, consisting, in some circumstances, of the functions that vanish up to
a certain order on a given variety. However, these are notoriously difficult to
compute, and there are seemingly simple questions related to symbolic powers that
remain open even over polynomial rings. In this paper, we describe a Macaulay2
software package that allows for computations of symbolic powers of ideals and
which can be used to study the equality and containment problems, among
others.