clear
set more off

cap log close
log using ex5, replace

set obs 1500
set seed 20097

* lognormal

g z=exp(invnorm(uniform())/2)

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(5)

log close
view ex5.smcl

exit