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