I come from the field of accelerator physics, where we study the dynamics of electrons traversing a circular ring guided by magnetic fields. Suppose one knows the one turn map which takes an electron from one phase space ($x_0$,$p_{x0}$ in one-D) point back to another ($x_1$,$p_{x1}$) and that one can express this map as a power series which terminates at some order $N$:
$x_1 = a_{10} x_0 + a_{01}p_0 + a_{20} x_0^2 + a_{11} x_0 p_0 +\dots + a_{nm}x_0^n p_0^m+\dots$ $p_1 = b_{10} x_0 + b_{01}p_0 + b_{20} x_0^2 + b_{11} x_0 p_0 +\dots + b_{nm}x_0^n p_0^m+\dots$
and $a_{nm}=b_{nm}=0$ for $n+m>N$. A typical case would be
$a_{10}=b_{01}=\cos\theta$, $a_{01}=-b_{10}=\sin\theta$,
for some real angle $\theta$. One may also impose additional requirements (in the case of non-radiating particles) that the map be symplectic, but the case with some damping is important as well. And this problem is the 1-D (2-D phase space) case, and the real problem is 3-D. The question is about the long term stability of this system. Can one formulate the stability region within the $x-p$ plane in terms of this map without actually computing the iterated map for each point? My sense is that a solution to this problem is still unknown to accelerator physicists (e.g. those working on the LHC and those working on 3rd generation light sources), but I wonder if it has been solved by mathematicians. Certainly for special cases, such as the Henon map and related maps much is known, but I don't know if a general solution exists. The method of normal form is the standard approach within these fields, but does not supply a solution. There are also methods for judging chaoticity of an orbit, but these methods are also inadequate from my understadning. For the case with damping (for light sources) the medium term stability (say 5000 terms) is of interest as well. Apologies if the problem is not well enough defined, but I wonder if someone has some knowledge of recent (or old!) solutions, or approaches to analysis.
Perhaps another useful thing to say is that it is assumed that the point $(0,0)$ is a fixed point of the map. So one may phrase this as finding the stability region expressed as a function of the polynomial map, in the vicinity of a fixed point.