i1 : R = QQ[a..e];
i2 : S = simplicialComplex {a*b*c*d*e};
i3 : faceDelete(a, S) o3 = | bcde | o3 : SimplicialComplex
i4 : faceDelete(a*b*c, S) o4 = | bcde acde abde | o4 : SimplicialComplex
i5 : faceDelete(a*b*c*d*e, S) == boundary S o5 = true