{1..d}]; 
    OrbN2:=OrbN diff {V1}; if #OrbN2 ge m then
    S2:={}; for l in {1..m-1} do
    RR:=Random(OrbN2 diff S2); S2:=S2 join {RR}; end for;
    S1:={@ q : q in S2 @}; PCM1:=P1; PCM2:=P1;
    for h in {1..N-1} do
    PCM1:=Matrix(F,1,d,[SpecialEvaluate(TT,PCM1[1][i]): i in 
    {1..d}])*CM; PCM2:=VerticalJoin(PCM2,PCM1); end for; 
    PCM3:=Transpose(PCM2); for ii in {1..m-1} do
    P:=Matrix(F,1,d,[S1[ii][j]: j in {1..d}]); PP:=P; PCM1:=P; 
    PCM2:=P; for h in {1..N-1} do
    PCM1:=Matrix(F,1,d,[SpecialEvaluate(TT,PCM1[1][i]): i in 
    {1..d}])*CM; PCM2:=VerticalJoin(PCM2,PCM1); end for; 
    PCM3:=HorizontalJoin(PCM3, Transpose(PCM2)); end for; 
    LC:=LinearCode(PCM3); print 
    "==========================================="; print " ";
    print "Skew Quasi Twisted Code of type: ", 
    N, N-d, MinimumWeight(LC); print "Polynomial generator:"; 
    print g1;  print " "; print "Companion Matrix:"; print CM; 
    print " "; print "Parity Check Matrix:"; print PCM3; print " ";
    print "==========================================="; else 
    print "There are not enough orbits of length", N ;
    end if; end if; return "end"; end function;

Consider $\F_8=\{ 0,1,w,w^2,\dots,w^6\}$. From Program \[Prog0\], by using the command +SD(7,w\^3)+ it follows that $g(x)=x^4-(w^3+x+w^2x^2)\in\F[x;\theta]$ divides $x^7-w^3$, where $\theta(x)=x^4$ for every $x\in\F_8$. Moreover, by Program \[Prog3\] with the command +SQT(2,7,\[w\^3,1,w\^2,0,1\])+ we see that the point $P_2:=[0:w^4:0:1]\in\P^3(\F_8)$ has an orbit of length $7$. So, we get a SQT $[14,4,7]_8$-code whose generator matrix is $$\begin{aligned}
 \begin{pmatrix} 
1 & 0 & 0 & 0 & w^3 & 