- Usage:
`differentials(L,j,g,s)`

- Outputs:
- a matrix

Generate the j-th differential of a Kustin-Miller resolution of length g. So, e.g., for j=1 we obtain the relations of the ring resolved and for j=2 the first syzygies of those.

We use the notation of Section 2.3 of

J. Boehm, S. Papadakis: On the structure of Stanley-Reisner rings associated to cyclic polytopes, http://arxiv.org/abs/0912.2152 [math.AC]

For any j the last entry of L should be the variable T.

For j=1 we assume L = {b_{1}, beta_{1}, a_{1}, T }.

For j=2 we assume L = {b_{2}, beta_{2}, h_{1}, a_{2}, alpha_{1}, T }.

For j=3,...,g-1 we assume L = {b_{j}, beta_{j}, h_{j-1}, a_{j}, alpha_{j-1}, b_{j-1}, T }.

For j=g-1 we assume L = {beta_{g-1}, h_{g-1}, a_{g-1}, alpha_{g-2}, b_{g-2}, T }.

For j=g we assume L = {alpha_{g-1}, a_{g}, b_{g-1}, u, T }.

Finally s equals k_{1}-k_{2}.

This is not really a user level function, however it is exported as occasionally it can be useful. The export may be removed at some point.

- kustinMillerComplex -- Compute Kustin-Miller resolution of the unprojection of I in J

- differentials(List,ZZ,ZZ,ZZ)