# S 5
sample(1:6, 1, replace=T)

sample(1:6, 10, replace=T)
mean(sample(1:6, 10, replace=T))

wuerfel <- NULL
for(j in 1:5000){
ergebnis = mean(sample(1:6, 10, replace=T))
wuerfel = c(wuerfel, ergebnis)
}

# S 6
mean(wuerfel)

# S 8
proben <- function(unten=1, oben = 6, k = 10, N = 50)
{
# default: wir werfen 10 Wuerfel 50 Mal
alle <- NULL
for(j in 1:N){
ergebnis = mean(sample(unten:oben, k, replace=T))
alle = c(alle, ergebnis)
}
alle
}


# S 9
o = proben(0, 99, 10, 50)
hist(o, col=3)

# S 10
hist(o, col=3, freq=F)


# S 12
osehrviele = proben(0, 99, 10, 50000)
h4 = hist(osehrviele, col=3, freq=F, breaks=200)


# S 15
sigma <- function(unten=1, oben = 6)
{
x = unten:oben
n = length(x)
m = mean(x)
sqrt((sum(x^2)/n - m^2))

}
sigma()

# S. 19
o = proben(0, 99, 10, 50)
hist(o, col=3, freq=F)
mu = mean(0:99)
sig = sigma(0, 99)/sqrt(10)

curve(dnorm(x, mu, sig), 30, 80, add=T)

# S. 23
curve(dnorm(x, 0, 1), -3, 3)

# S. 29
mu = mean(1:6)
SE = sigma(1, 6)/sqrt(9)
pnorm(2.9, mu, SE)
qnorm(0.1459479)
mu + qnorm(0.1459479) * SE

# S. 30
mu = mean(1:6)
SE = sigma(1, 6)/sqrt(9)
pnorm(2.9, mu, SE)