Fig_15

pdf("fig_15.pdf", width = 5, height = 5)
par(bty = "n", las = 0, mar = c(4, 4, 2, 1) + 0.2)

# función suma armónicos
sum_arm <- function(n) {
  a <- numeric(n)
  for (i in 1:n) {
    for (j in 1:i) {
      a[j] <- a[j] + 1/i
    }
  }
  a/n
}

longitud <- 10

plot(NA, NA, xlab = "CV corregido", ylab = "Densidad de probabilidad", xlim = c(0, 2), ylim = c(0, 3))

nn <- c(5, 10, 20, 40)
#colores <- rainbow(4)
colfunc <- colorRampPalette(c("orange", "royalblue"))
colores <- colfunc(4)
tipo = 0
for (n in nn) {

      tipo <- tipo + 1
      repetir <- 100000
      cvS <- numeric(repetir)
      trozosS <- numeric(n)
      for (j in 1:repetir){
        p_azar <- longitud*runif(n - 1)    # elegir los puntos de corte al azar
        p_azar <- sort(c(0, p_azar, longitud))               # ordenar los cortes
        for (i in 1:n){
          trozosS[i] <- p_azar[i + 1] - p_azar[i]
        }
        trozosS <- sort(trozosS, decreasing = TRUE)
        cvS[j] <- sd(trozosS)/mean(trozosS)*n/(n - 1)
      }

      # histograma
      h <- density(cvS)
      lines(h, lwd = 3, col = colores[tipo])

}

legend("topright", c("n = 5", "n = 10", "n = 20", "n = 40"), bty = "n", lwd = 2, col = colores, cex = 0.8)

dev.off()