read RotaStanley: L:=BooleanLattice(2): UmbralRecurrenceSingle(L,1,q,[x],[a]): print(`This took`, time(), `seconds of CPU time `): quit: