HenonStokesConst-1-3
A set of scripts that compute the Stokes constant for 2 instances of the Henon map. I assume that someone who is interested in this already knows what is the point of the Stokes constant, so I will not explain anything here. There will be a paper that explains this and the procedure soon, unfortunately it is not available yet.
The equations were embedded as images since there is no native support for them. They are ugly and misalligned, but it will make do.
The folder Henon
has the scripts that compute the Stokes constant of the map
with
and the folder RotHenon
has the scripts that compute the Stokes constant of the map
with T_{2π/3} the 2π/3 rotation matrix.
For both maps the origin is a fixed point at 1:3 resonance and there exists a formal series that satisfies formally the equation
with
Let , t>0 and, then the Stokes constant is the double limit
The script computes an approximation of this limit for fixed t, n and m.
The script (Rot
)HenonCoeffsComp.wl
computes the constants w_n. It is written in Mathematica and it is not as optimized as it should be. Each time is being run it computes 10 more coefficients and stores them in the file (Rot
)HenonCoeffs.dat
. Then the sript (Rot
)HenonCoeffsToJulia.wl
writes them in the file (Rot
)HenonCoeffsList.jl
in a form that JuliaLanguage can understand.
The script (Rot
)HenonConst.jl
computes the Stokes constant of the map. There are some comments that state which variables can be redefined safely. The result is written in the file \texttt{(Rot
)HenonConst.dat}. Each line has the format:
[n] [real part] [error of the real part] [imaginary part] [error of the imaginary part]