smpl 1 100 genr white=nrnd smpl 1 2 genr archmy=0 genr archmerr=0 genr archmh=0 genr mean=0 smpl 3 100 genr archmh =.8+.5*(archmerr(-1))^2+.3*(archmerr(-2))^2 genr mean=1+.4*archmh genr archmerr=white*((archmh)^.5) genr archmy=mean+archmerr