%opony
dane = [32, 33, 28, 37, 29, 30, 25, 27, 39, 40, 26, 26, 27, 30, 25, 30, 31, 29, 24, 36, 25, 37, 37, 20, 22, 35, 23, 28, 30, 36, 40, 41];
mu=mean(dane)
sigma=std(dane)
wynik=[(mu-tinv(0.005,31)*sigma/sqrt(32)) (mu +tinv(0.005,31)*sigma/sqrt(32)) ]