clear
set more off

cap log close
log using ex4, replace

set obs 500
set seed 11405

* two humps, constrained variance

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

log close
view ex4.smcl

exit