D = pruneUnit(C, n, u, M)
Prune the unit in location {u} of the {n}-th differential in chain complex {C}. Moves the unit to the end, clears row and column, while keeping the neighboring differentials compatible.
Currently only works with 1x1 units.
The object pruneUnit is a method function with options.