(wip) add frequentist intervals

Add 1D frequentist intervals, calculated by running many MCs at each point
in parameter space to decide if it is in the acceptance region. To avoid
grid scanning, we perform a root-finding binary search; still tweaking this
to improve convergence.
merge bonventre/threef 162854afd6