> Zeck(30); rz(30); compact(rz(30));
                                         [8, 6, 2]
                                         [0, 2, 0]
                                            0 2 0
> rz(30);
                                        [0, 2, 0]
> for i from 1 to 10 do print(i,rz(F(7*i)/F(7))) od;
                                           1, [0]
                                         2, [0, 4]
                                  3, [1, 0, 0, 0, 0, 0, 0]
                                4, [1, 0, 0, 0, 0, 1, 1, 3]
                               5, [1, 0, 0, 0, 0, 1, 1, 8, 0]
                             6, [1, 0, 0, 0, 0, 1, 1, 9, 0, 4]
                      7, [1, 0, 0, 0, 0, 1, 1, 9, 1, 0, 0, 0, 0, 0, 0]
                    8, [1, 0, 0, 0, 0, 1, 1, 9, 1, 0, 0, 0, 0, 1, 1,
3]
                   9, [1, 0, 0, 0, 0, 1, 1, 9, 1, 0, 0, 0, 0, 1, 1, 8,
0]
                 10, [1, 0, 0, 0, 0, 1, 1, 9, 1, 0, 0, 0, 0, 1, 1, 9,
0, 4]
> for i from 1 to 30 do print(compact(rz(F(i)^3)))od;
                                              0
                                              0
                                              4
                                            1 1 0
                                            0 4 1
                                         1 0 1 0 0 0
                                         1 0 0 0 2 2
                                         1 0 0 1 6 0
                                       1 0 0 1 1 0 3 1
                                     1 0 0 1 2 0 2 0 0 0
                                      1 0 0 1 1 3 3 0 2
                                      1 0 0 1 1 6 1 4 0
                                   1 0 0 1 1 5 0 0 0 1 1 1
                                  1 0 0 1 1 5 1 2 1 0 0 0 0
                                    1 0 0 1 1 5 0 4 5 0 2
                                 1 0 0 1 1 5 1 0 1 0 1 1 4 0
                               1 0 0 1 1 5 1 0 0 0 3 0 0 1 1 1
                               1 0 0 1 1 5 1 0 0 1 7 1 0 0 0 0
                               1 0 0 1 1 5 1 0 0 1 1 0 3 5 0 2
                             1 0 0 1 1 5 1 0 0 1 2 0 2 0 1 1 4 0
                            1 0 0 1 1 5 1 0 0 1 1 3 3 1 0 0 1 1 1
                            1 0 0 1 1 5 1 0 0 1 1 6 1 5 1 0 0 0 0
                           1 0 0 1 1 5 1 0 0 1 1 5 0 0 0 1 1 5 0 2
                          1 0 0 1 1 5 1 0 0 1 1 5 1 2 1 0 0 1 1 4 0
                          1 0 0 1 1 5 1 0 0 1 1 5 0 4 5 1 0 0 1 1 1
                       1 0 0 1 1 5 1 0 0 1 1 5 1 0 1 0 1 1 5 1 0 0 0 0
                       1 0 0 1 1 5 1 0 0 1 1 5 1 0 0 0 3 0 0 1 1 5 0 2
                       1 0 0 1 1 5 1 0 0 1 1 5 1 0 0 1 7 1 0 0 1 1 4 0
                     1 0 0 1 1 5 1 0 0 1 1 5 1 0 0 1 1 0 3 5 1 0 0 1 1
1
                   1 0 0 1 1 5 1 0 0 1 1 5 1 0 0 1 2 0 2 0 1 1 5 1 0 0
0 0
> seq(F(n) mod 13, n=0..30);
0, 1, 1, 2, 3, 5, 8, 0, 8, 8, 3, 11, 1, 12, 0, 12, 12, 11, 10, 8, 5,
0, 5, 5, 10, 2, 12, 1,

    0, 1, 1
> period(13);
                                             28
> for m from 1 to 20 do print(m, period(m)) od;
                                            1, 1
                                            2, 3
                                            3, 8
                                            4, 6
                                           5, 20
                                           6, 24
                                           7, 16
                                           8, 12
                                           9, 24
                                           10, 60
                                           11, 10
                                           12, 24
                                           13, 28
                                           14, 48
                                           15, 40
                                           16, 24
                                           17, 36
                                           18, 24
                                           19, 18
                                           20, 60
> for k from 1 to 20 do 
>   print(ithprime(k), period(ithprime(k))) 
> od;

                                 2, 3


                                 3, 8


                                5, 20


                                7, 16


                                11, 10


                                13, 28


                                17, 36


                                19, 18


                                23, 48


                                29, 14


                                31, 30


                                37, 76


                                41, 40


                                43, 88


                                47, 32


                               53, 108


                                59, 58


                                61, 60


                               67, 136


                                71, 70

> 

                             3, 6, 12, 24

> Zp(0,2);
                                            [2]
> Zp(1,2);
                                            [2]
> Zp(2,2);
                                            [3]
> Zp(0,3);
                                           [5, 3]
> Zp(1,3);
                                           [6, 4]
> Zp(2,3);
                                           [7, 5]
> Zp(0,4);
                                            [3]
> Zp(1,4);
                                            [4]
> Zp(0,5);
                                   [16, 13, 11, 9, 6, 3]
> Zp(0,6);
                                 [20, 15, 13, 11, 8, 6, 3]
> for i from 3 to 10 do print(`m=`.i,`period=`.(period(i)),Zp(0,i)) od;
                                   m=3, period=8, [5, 3]
                                     m=4, period=6, [3]
                           m=5, period=20, [16, 13, 11, 9, 6, 3]
                         m=6, period=24, [20, 15, 13, 11, 8, 6, 3]
                               m=7, period=16, [11, 9, 7, 5]
                                   m=8, period=12, [7, 5]
                           m=9, period=24, [19, 15, 13, 11, 9, 5]
      m=10, period=60, [55, 50, 46, 44, 39, 37, 33, 31, 29, 27, 23,
20, 18, 12, 9, 5]
>