  We have developed an orbit-based method for constructing triaxial models of
elliptical galaxies, which fit their observed surface brightness and
kinematics. We have tested the method against analytical models with general
distribution functions. Here, we present models that fit integral-field SAURON
observations of NGC3379 and NGC4365.
