>gauss<-rnorm(100,2,6)

>boot<-0

>for(i in 1:1000){

+boot=mean(sample(gauss,replace=T))

+}

>summary(boot)

   Min. 1st Qu. Median    Mean 3rd Qu.    Max.

0.3345 1.9540  2.3350  2.3230 2.7020  4.2330

>summary(gauss)

   Min. 1st Qu. Median    Mean 3rd Qu.    Max.

-13.380  -2.238  2.570   2.296   6.861 16.230

>sd(boot)

[1]0.599087

>sd(gauss)/sqrt(100)

[1]0.5906275