This is a Mathematica code for computing higher cup products by Steenrod. Our notation is as follows.
A^{(p)} ⌣_i B^{(q)} ⌣_j C^{(r)}
= Cup[A[p], B[q], C[r]][{i, j}]
(Error for wrong numbers of forms and cups.)
Associativity is supported in the appropriate way; an overall factor will be factorized.
The exterior derivative is implemented as Del
as the usual Mathematica code;
here Cup
, which mixes some other forms and higher cups appropriately;
it is surely nilpotent and satisfies homomorphism.
The Leibniz rule holds.