double point = gsl_ran_binomial_pdf (mc[i].n01, 0.5,
mc[i].n01 + mc[i].n10);
double entries[] = { n, 2.0 * sig, sig, point };
double point = gsl_ran_binomial_pdf (mc[i].n01, 0.5,
mc[i].n01 + mc[i].n10);
double entries[] = { n, 2.0 * sig, sig, point };