c-d.r

mu <- 2.5
tau <- 0.5
xn <- rnorm(50, mu, 1/tau)

# test for convergence
library(purrr)

mu_i = 0
alpha = 0.01

for(i in 1:10000000){
  mu_i <- mu_i - alpha * (tau *(mu_i - mean(xn)))
  print(mu_i)
}