1 | initial version |

One possibility to implement the algebra $H_4$ from *loc. cit.*, the one where we also use the relations $x^2=0$, $g^2=1$, $gx+xg=0$ (and there is no need for an inverse, since $g^{-1}=g$) would be:

```
A.<X,G> = FreeAlgebra(QQ, 2) # or over QQbar
F = A.monoid()
X, G = F.gens()
monomials = [F(1), X, G, X*G]
MS = MatrixSpace(QQ, len(monomials)) # or respectively over QQbar
matrices = [
# matrix showing the action of the first generator, X, on the monomials
MS([0, 1, 0, 0, # X*1 = X
0, 0, 0, 0, # X*X = 0
0, 0, 0, 1, # X*G = XG
0, 0, 0, 0, # X*XG = (XX)G = 0
]),
# matrix showing the action of the second generator, G, on the monomials
MS([0, 0, 1, 0, # G*1 = G
0, 0, 0, -1, # G*X = -XG
1, 0, 0, 0, # G*G = 1
0, -1, 0, 0, # G*XG = (GX)G = (-XG)G=-X(GG)=-X
]),
]
H4.<x,g> = A.quotient(monomials, matrices)
```

Then we have for instance:

```
sage: (x+g)^2
1
sage: (1+x*g)^2
1
sage: (1+x+g)^10
512 + 512*x + 512*g
```

*Note:* It would be interesting to have a full working algebraic setting for the Hopf algebra $H_4$ (including $S$ and $\Delta$). Well, posting as anonymous is not a good idea in general, but ok, here is the state of the art so far.

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.