GAP 4.8.9 installation with standard packages -- copy to your CoCalc project to get it

1
\\ compute basis of maximal order
2
b = nfbasis( f );
3

4
p2v(n,b)=vector(n,j,polcoeff(b,j-1));
5

6
\\ print result 
7
print("[ ");
8
for(i=1,#b, print(p2v(#b,b[i]),","));
9
print("];");
10

11

12