Summary_2014_10_09

1. The rates of the acceptance and rejection are only estimated after
the equilibrium.
2. To accelerate the computation, we introduce the global variant
"former_energy". The value will be changed only if the new configuration
is accepted. So as to the energy.
3. The radius distribution function is rewritten. And we add the other
head files. The G[] is also defined as the global array with
initialization, statistic and writing to files overall 3 modes.