To further investigate the behavior of the chaotic boundary we set up three new set of simulations.  In the first simulation we investigate the behavior of whole libration island with a low-resolution grid of particles.   The initial conditions were chosen with similar criterions for the low-resolution survey shown in S2000/S5: limits of the libration island and an application of the frequency analysis method,  only this time the initial values of the inclination was of 36.83^o,(with a step of 0.6^o) while we started with \omega =40.146^o (and a step of 2.5^o),to include a larger region of circulating orbits.   In the second survey we concentrate on the lower boundary of the separatrix, with initial values of inclination of 41.2^o (the step is of 0.02^o in this case) and an initial value of \omega of 60.14^o.  In the third survey we study the region where S/2000S5 reach the maximum of the eccentricity.   The initial value of \omega is still 60.14^o, but the initial inclination is 41.675^o in this case.
 

For the filtering process we choose these parameters:

3
BAB
2 10 3
80 0.024d0 10.0d0
80 0.10d0 20.0d0
182.625d0
bin.filter.dat
F ! write real*8 bin.dat output of osculating elements? [T|F] (useful for precise restart)

This set of parameters give a sampling time of 30 yr (which corresponds to a resolution of 1.32 "/yr), a passband of 200 yr (=1/2 of the smallest precession period of \omega), and a stop-band of 66.6 yrs, more than enough to stop frequencies associated with the orbital period of Saturns (29.6 yr).