Handle Bessel correction properly

- distinguish between N and N-1 and do not approximate N-1 with N
Change interface of MC helper functions

- cumulative_result and chi_squred_dof no longer need to be called with
  the numeric type