20 years after the discovery of the first planets outside our solar system,the current exoplanetary population includes more than 700 confirmed planets around main sequence stars.Approximately 50％ belong to multiple-planet systems in very diverse dynamical configurations,from two-planet hierarchical systems to multiple resonances that could only have been attained as the consequence of a smooth largescale orbital migration.The first part of this paper reviews the main detection techniques employed for the detection and orbital characterization of multiple-planet systems,from the (now) classical radial velocity (RV) method to the use of transit time variations (TTV) for the identification of additional planetary bodies orbiting the same star.In the second part we discuss the dynamical evolution of multi-planet systems due to their mutual gravitational interactions.We analyze possible modes of motion for hierarchical,secular or resonant configurations,and what stability criteria can be defined in each case.In some cases,the dynamics can be well approximated by simple analytical expressions for the Hamiltonian function,while other configurations can only be studied with semi-analytical or numerical tools.In particular,we show how meanmotion resonances can generate complex structures in the phase space where different libration islands and circulation domains are separated by chaotic layers.In all cases we use real exoplanetary systems as working examples.