read RotaStanley:

for i from 1 to 3 do
print(`The generating function for the`, i, `-dim. Boolean latticee in q is:`):
lu:=GFchainQ(BooleanChain(i),q):
lprint(lu):
gu:=taylor(lu,q=0,41):
gu:=[seq(coeff(gu,q,i1),i1=0..40)]:
print(`The first 40 terms are `):
print(gu):
od:

time();


quit: