add subroutine for the calculation of statistical uncertainty of the density of states