Abstract
Based on recent observational data on transverse and radial velocities of the bright Orion Trapezium stars we study the dynamical evolution of ensembles of multiple systems mimicking the Trapezium. To this end we conduct numerical N-body integrations using the observed planar positions and velocities, the radial velocities, and random z-positions for all components. We include perturbations in these quantities compatible with the observational errors. We briefly discuss the properties of the evolved systems and of the resulting binaries.