reset
#set terminal pngcairo size 350,262 enhanced font 'Verdano,6'

#set style line 1 lc rgb '#0060ad' lt 1 lw 1 pt 7 ps 1.5 # --- blue

set xrange [-40.:40.]
set yrange [-0.8:0.8]
set xlabel 'x[mm]'
set ylabel "x' [mrad]"
set label 'septum' at 27,0.1
set label 'GeV' at screen 0.26,0.15
unset arrow
set key off
set arrow from 31. ,0.2 to 31.,.6 lw .5 nohead
omega=2*pi*60
delta = 0.15/5.3
t0 = acos(1.-2*delta)/omega
do for [j=1:4000:1]{
set label 98 'turn '.int(j) at screen 0.2,0.2
energy = 5.300/2 * (1.-cos(j * omega *2.52131e-6 +omega*t0 ))
set label 99 sprintf('%01.3f',energy) at screen 0.2, 0.15
outfile = sprintf('animation_4000/resonantextraction%04.0f.png',j)
#set output outfile
#print outfile
plot dir.'phasespace.dat' u ($1*1000):($2*1000/(1+$6)) i j:j pt 6 ps 1.
#pause 0.001
reread
}