1gap> #Constructing form: HermitianFormByPolynomial
2gap> r := PolynomialRing( GF(9), 4);
3GF(3^2)[x_1,x_2,x_3,x_4]
4gap> vars := IndeterminatesOfPolynomialRing( r );
5[ x_1, x_2, x_3, x_4 ]
6gap> poly := vars[1]*vars[2]^3+vars[1]^3*vars[2]+
7>              vars[3]*vars[4]^3+vars[3]^3*vars[4];
8x_1^3*x_2+x_1*x_2^3+x_3^3*x_4+x_3*x_4^3
9gap> form := HermitianFormByPolynomial(poly,r);
10< hermitian form >
11gap> Display(form);
12Hermitian form
13Gram Matrix:
14 . 1 . .
15 1 . . .
16 . . . 1
17 . . 1 .
18Polynomial: x_1^3*x_2+x_1*x_2^3+x_3^3*x_4+x_3*x_4^3
19
20gap> quit;
21