We study the ions dynamics inside an Electrostatic Ion Beam Trap (EIBT) and show that the stability of the trapping is ruled by a Hill's equation. This unexpectedly demonstrates that an EIBT, in the reference frame of the ions works very similar to a quadrupole trap. The parallelism between these two kinds of traps is illustrated by comparing experimental and theoretical stability diagrams of the EIBT. The main difference with quadrupole traps is that the stability depends only on the ratio of the acceleration and trapping electrostatic potentials, not on the mass nor the charge of the ions. All kinds of ions can be trapped simultaneously and since parametric resonances are proportional to the square root of the charge/mass ratio the EIBT can be used as a mass spectrometer of infinite mass range. , they form a new family of traps operating at energies of a few keV. They are used for atomic and molecular metastable-states studies, molecular fragmentation and photodissociation (see, e.g., [7] for a review). Beyond providing trapping of energetic particles in a well defined direction, these traps have many interesting features : they are small, relatively inexpensive, easy to setup and operate and have a field-free region where ions move freely and where measurements can easily be performed. They can even be used as Time Of Flight (TOF) mass spectrometers [8] or cooled at cryogenic temperatures [9].Despite these interesting features, all the published theoretical models describing EIBT are based on one dimensional approximations, neglecting the radial motion. This leads to inaccurate predictions of the trap stability and operating domain, and restricts their flexibility, as finding reasonable working points requires lengthy and tedious experimental exploration. This may partly be explained by the lack of an analytical formula for the electrostatic potential inside these traps leading to the dilemma of choosing between a simplistic analytical model and a heavy numerical treatment unsuited to explore the huge space of parameters. Usual beam simulation codes fail to produce good results as the numerical inaccuracies at the ion turning points lead to energy non conservation reaching a few 100 eV over a few tens of oscillations. Here we solve the problem, using methods developed for radio-frequency quadrupole traps. We show that the radial dynamic is ruled by a Hill's equation, a particular case of the Mathieu equation that describes quadrupole traps. This model yields accurate predictions of the ions motion in the trap and of the stability region.The design and operation of the EIBT has been des- cribed previously in detail [1] and a schematic drawing of the ion trap is shown in Fig. 1. The trap consists of two sets of coaxial cylindrical electrodes roughly equivalent to two spherical mirrors, the electrostatic analog of a Fabry-Perot interferometer. The configuration of the trap is defined by the potentials of five of these electrodes {V 1 , V 2 , V 3 , V 4 , V z } on both ends of the trap, the others being grounded. ...