# To get A261422 out to 200000 terms
N:=6: # To get all terms of A002113 of up to N digits in the list "Res"
Res:= $0..9:
for d from 2 to N do
  if d::even then
    m:= d/2;
    Res:= Res, seq(n*10^m + digrev(n), n=10^(m-1)..10^m-1);
  else
    m:= (d-1)/2;
    Res:= Res, seq(seq(n*10^(m+1)+y*10^m+digrev(n), y=0..9), n=10^(m-1)..10^m-1);
  fi
od:
Res:=[Res]:
G:=add(x^Res[i],i=1..nops(Res)):
#M:=Res[nops(Res)]+2:
M:=200003;
G1:=series(G,x,M):
G2:=series(G1^2,x,M):
G3:=series(G1*G2,x,M):
L3:=seriestolist(G3):
writeto(trans);
for n from 1 to nops(L3) do lprint(n-1,L3[n]); od: