program define monte
version 2.1
set seed 31977
set more 1
log using monte
macro define _it 1
while %_it<201 {
     quietly {
          drop _all
          set obs 50
          generate X=invnorm(uniform())
          summ X, detail
          }
     display %_it
     display _result(3)
     display _result(10)
     display
     macro define _it=%_it+1
     }
log close
drop _all
infile iterate mean median using monte.log
label data "Monte Carlo experiment #1"
label variable mean "Sample mean, X normal"
label variable median "Sample median, X normal"
save monte
end