capture log close 
log using 771bs1, replace
set more off
set seed 123456789
local reps 1000

webuse census2, clear
summ age
ci age
bootstrap "summ age" r(mean) r(sd), reps(`reps') saving(771bs1) replace

use 771bs1,clear
graph drop _all
histogram _bs_1, normal name(g771bs1_m) 
histogram _bs_2, normal name(g771bs1_s)
graph combine g771bs1_m g771bs1_s, col(1) ysize(7) xsize(9) ///
   ti("Bootstrap distributions of mean and standard error")
graph export 771bs1.pdf, replace
log close
log2html 771bs1.smcl, replace