p:=5; n:=2; wp:=-7; M:=MatrixRing(Integers(),p-1); z:=M ! CompanionMatrix(CyclotomicPolynomial(p)); a:= Integers() ! PrimitiveElement(Integers(p)); Z:=z; for i in [1..p-2] do z:=z^a; Z:=DiagonalJoin(Z,z); end for; e:=M ! 1; pi:=M ! 0; for i in [1..p-2] do pi[i][i+1] := 1; end for; pi[p-1][1]:= n; P:=TensorProduct(pi,e); MM:=MatrixRing(Integers(),(p-1)^2); E:=MM ! 1; MM2:=MatrixRing(Integers(),2*(p-1)^2); UU:=sub; w:=MatrixRing(Integers(),2) ! [0,1,wp,0]; e:=MatrixRing(Integers(),2) ! 1; UUw:=sub; SchurIndexJac(UUw);