import numpy as np
import math

# cutoff for particle number:
nmax = 50
jmax = 50
mmax = 50
kmax = 50

# system rates:
rho = 1.

# other constants:
#g = 8. #shouldn't be a constant!
qplus = 12.
qminus = 1.
n0 = 4. #in paper
nu = 4.
