Dear all, I am trying to plot the cumulative number counts of galaxy cluster (N>S) as a function of X-ray flux. The method I am following is using the mass function and scaling relation to calculate dN/dVdz. Then I do a double integrate over the redshift range and flux range. But the plot that I get has no connection with observational results. I am wondering if you could have a look at the codes and help me wo solve the problem.