// likelihood = P(E | H)
// prior = P(H)
// modelEvidence = P(E)
// returns P(H | E)
static double bayes(double likelihood, double prior, double modelEvidence) {
  ret likelihood*prior/modelEvidence;
}