clear
set more off

cap log close
log using ex1, replace

set obs 500
set seed 60157

* two humps

g index=uniform()
g z=invnorm(uniform())+2
replace z=-z if index<0.4

kdensity z, s(.) n(100) xlab ylab gap(4)

* set more on
* more 
* set more off

cap noi denormix z, loglevel(3) icomp chi2 nmax(6)

log close
view ex1.smcl

exit