log using cauchy_log ,replace local i = 10 while `i' < 10001 { capture drop_all set obs `i' gen normal`i' = invnorm(uniform()) gen cauchy`i' = normal`i' / invnorm(uniform()) ci normal`i' cauchy`i' local i = `i' * 10 } log close