gap> P:=PoincareSeries(R,20);
(1)/(-x_1^3+3*x_1^2-3*x_1+1)

gap> ExpansionOfRationalFunction(P,2000)[2000];
2001000
