def size(N, nb):
"""
Return the size of the boson many-body basis.
Parameters
----------
N : int
Number of sites
nb : int
Number of bosons
"""
return int(binom(nb + N - 1, nb))