p = 5 Here the gamma-values gamma[1], gamma[2], gamma[3], gamma[4], `...` are 4, 4, 4, 4, () .. () 

 

 

 

>    M := 2:

   p := 5:

bound := 6:

 for alpha to bound do

          GF||alpha := PI(p^alpha, M, 1):
         ord||alpha := order(GF||alpha, p^alpha):
      
od: print(``);

print(array([

['p', ``, 'alpha', ``, ``, 'p^alpha', ``, ``, '(({p^alpha-1}/2))[p]!', ORDER],

seq([p, ``, alpha, ``, ``, ifactor(p^alpha), ``, ``, GF||alpha, ord||alpha],

alpha = 1..bound)])): print(``);

lprint(`See how the order value is FIXED at 4.`);
 

 

 

 

array( 1 .. 7, 1 .. 10, [( 3, 2 ) = (``), ( 5, 2 ) = (``), ( 4, 6 ) = (`*`(`^`(``(5), 3))), ( 7, 9 ) = (-1068), ( 5, 3 ) = (4), ( 7, 4 ) = (``), ( 5, 9 ) = (182), ( 5, 10 ) = (4), ( 4, 1 ) = (5), ( 7,...
`See how the order value is FIXED at 4.`
 

>