  The Broad Histogram Method (BHM) allows one to determine the energy
degeneracy g(E), i.e. the energy spectrum of a given system, from the knowledge
of the microcanonical averages <Nup(E)> and <Ndn(E)> of two macroscopic
quantities Nup and Ndn defined within the method. The fundamental BHM equation
relating g(E) to the quoted averages is exact and completely general for any
conceivable system. Thus, the only possible source of numerical inaccuracies
resides on the measurement of the averages themselves. In this text, we
introduce a Monte Carlo recipe to measure microcanonical averages. In order to
test its performance, we applied it to the Ising ferromagnet on a 32x32 square
lattice. The exact values of g(E) are known up to this lattice size, thus it is
a good standard to compare our numerical results with. Measuring the deviations
relative to the exactly known values, we verified a decay proportional to
1/sqrt(counts), by increasing the counter (counts) of averaged samples over at
least 6 decades. That is why we believe this microcanonical simulator presents
no bias besides the normal statistical fluctuations. For counts~10**10, we
measured relative deviations near 10**(-5) for both g(E) and the specific heat
peak, obtained through BHM relation.
