The statistics of stellar systems of multiplicity three and higher is reviewed. They are frequent, 0.15−0.25 of all stellar systems. Some 700 multiples are expected among the 3383 stars of spectral type F, G, and K within 50 pc, while only 76 of them are actually known. Many (if not all) close binaries have distant tertiary components, indicating that angular momentum exchange within multiple systems was probably critical in forming short-period binaries. The ratio of outer to inner periods in the best-studied nearby multiples and in low-mass pre-main sequence multiples does not exceed 104 at the formation epoch; larger ratios are produced by subsequent orbital evolution. All multiples with well-defined orbits are dynamically stable, the eccentricities of outer orbits obey the empirical stability limit Pout(1 – eout)3/Pin > 5 that is more strict than current theoretical limits. Relative orientation of orbits in triple stars shows some degree of alignment, especially in weakly-hierarchical systems. The statistics support the idea that most multiple stars originated from dynamical interactions in small clusters.