# Daten
d <- c(163, 185, 174, 180, 175, 179)

# Angabe
md <- mean(d)
sdd <- sd(d)
sapply(c(0.5,0.75,0.9,0.95,0.975),
       qnorm,mean=md,sd=sdd)

# Loesung
# b)
mean(d)
var(d)*5/6

# c)
pdf("figA20lsg.pdf")
plot(ecdf(d),ylab="Fn(z)",xlab="z")
lines(c(161,164,166,171,176,181,186,188,191),
      c(0.025,0.05,0.1,0.25,0.5,0.75,0.9,0.95,0.975))
dev.off()