Compute probability of observations across many potential vectors This will try (1/stepsize)^nflips possible vectors, computing the probability of the observation for each

try_many_vectors(
nflips,
number_samples = 1000,
stopval = 1e-05,
log = FALSE,
)