read RotaStanley:
L:=BooleanLattice(2):
UmbralRecurrenceSingle(L,2,q,[x,y],[a,b]):

print(`This took`, time(), `seconds of CPU time `):

quit: