请修改下面的代码使它看起来不一样,注意将变量和函数命名更换:average_m<-function(N){ set.seed(42) M<-rep(0,N) for(i in 1:N){ m<-2 U1<-runif(1,min=0,max=1) U2<-runif(1,min=0,max=1) while(U1<=U2){ m<-m+1 U1<-U2 U2<-runif(1,min=0,max=1) } M[i]=m } I<-mean(M) S<-var(M) return(list(I,S)) } N<-10000 result<-average_m(N) s<-sqrt(as.numeric(result[2])/N)

这是修改后的代码，我将变量和函数命名更换了，同时也稍微修改了一下代码风格：

compute_avg_and_var <- function(num_simulations) {
  set.seed(42)
  simulation_results <- rep(0, num_simulations)
  for (i in 1:num_simulations) {
    current_value <- 2
    U1 <- runif(1, min = 0, max = 1)
    U2 <- runif(1, min = 0, max = 1)
    while (U1 <= U2) {
      current_value <- current_value + 1
      U1 <- U2
      U2 <- runif(1, min = 0, max = 1)
    }
    simulation_results[i] <- current_value
  }
  mean_value <- mean(simulation_results)
  variance <- var(simulation_results)
  return(list(mean_value, variance))
}

num_simulations <- 10000
result <- compute_avg_and_var(num_simulations)
standard_deviation <- sqrt(as.numeric(result[2]) / num_simulations)