#5. Praktika

#1. ariketa

#banaketa exponentziala

#a) 

#EX = 1/?? VAR(X) = 1/??^2
Lambda = 1/5
EX = 1/Lambda
VARX = 1/(Lambda^2)


#b)
c <- rexp(10,Lambda)

hist(c, xlim = c(0,20), xlab = "Itxaro-denbora", col = "pink")

bb <- mean(c)

#c)
matrizea <- c
for(i in 2:500){
  c1 <- rexp(10,Lambda)
  matrizea <- rbind(matrizea,c1)
}

bbm <- rowMeans(matrizea)

hist(bbm, xlim = c(0,20), xlab = "Itxaro-denbora", col = "pink")

bb2 <- mean(bbm)

#d)
matrizea1 <- rexp(10,Lambda) 
for(i in 2:500){
  c2 <- rexp(10,Lambda)
  matrizea1 <- rbind(matrizea1,c2)
}

bbm1 <- rowMeans(matrizea1) # Bai, balioak 5en inguruan baitaude

dents1 <- dnorm(bbm1, mean(bbm1), var(bbm1))

dents2 <- dnorm(EX,VARX)

plot(dents1)

lines(dents2)