  We describe method for calculating the density of states by combining several
canonical monte carlo runs. We discuss how critical properties reveal
themselves in $g(\epsilon)$ and demonstrate this by applying the method several
different phase transitions. We also demonstrate how this can used to calculate
the conformal charge, where the dominating numerical method has traditionally
been transfer matrix.
