bruteForcePsi.MLE = function() {

  z = 1
  n = 5
  eps = 1e-6
  psi = seq(eps,1-eps, by=eps)
  logLike = z*log(psi) + (n-z)*log(1-psi)
  psi[logLike == max(logLike)]
}