clear set more off cap log close log using ex2, replace set obs 500 set seed 38901 * two humps, random starting point 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(3) noinit log close view ex2.smcl exit