# Renvoie sigma^2
est_variance = function(X) {
  m = mean(X)
  somme = 0
  for (x in X) {
    somme = somme + (x - m) * (x - m)
  }
  somme = somme / (count(X) - 1)
  return(somme)
}

print(sigma)
print(n)
k_alpha = function(mu_0, sigma, n, alpha) {
  return(mu_0 + sigma / sqrt(n) * qnorm(1 - alpha))
}

delta_1 = function(S, alpha, mu_0, mu_1) {
  k = k_alpha(mu_0, sigma, n, alpha)
  if (mean(S) > k)
    return(0)
  else
    return(1)
}