fixed errors

1. In order to have reproducible results using multicores, I need to
load the library `parallel` first, and use L'Ecuyer-CMRG random number
generation.
2. I added a function to calculate a bootstrap confidence interval
3. Changed the names of the generated files
4. Increased the number of bootstrapped iterations