  We use the Made-to-Measure method to construct N-body realizations of
self-similar, triaxial ellipsoidal halos having cosmologically realistic
density profiles. Our implementation parallels previous work with a few
numerical refinements, but we show that orbital averaging is an intrinsic
feature of the force of change equation and argue that additional averaging or
smoothing schemes are redundant. We present models having the Einasto radial
mass profile that range from prolate to strongly triaxial. We use a
least-squares polynomial fit to the expansion coefficients to obtain an
analytical representation of the particle density from which we derive density
contours and eccentricity profiles more efficiently than by the usual particle
smoothing techniques. We show that our N-body realizations both retain their
shape in unconstrained evolution and recover it after large amplitude
perturbations.
