#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 term x11
set key off
#set xrange [0:755.8698]
set yrange [-2.:3.]
set xlabel 's [m]'
set ylabel 'x [cm]'
set label 'GeV' at screen 0.26, 0.15
unset arrow
set arrow from 531.344, 2.7 to 531.344, 2.9 lw 2 nohead 
set label 'septum' at 550.,2.7
#set arrow from 0.03,-0.00004 to 0.03,0.00004 lw 2 nohead
omega=2*pi*60
delta = 0.15/5.3
t0 = acos(1.-2*delta)/omega
do for [j=1:400:1]{
set label 99 'turn '.int(j*10) at screen 0.2,0.2
energy = 5.300/2 * (1.-cos(j*10 * omega *2.52131e-6 +omega*t0 ))
set label 98 sprintf('%01.3f',energy) at screen 0.2, 0.15
#outfile = sprintf('animation_4000/orbitdistortion%04.0f.png',j)
#set output outfile
#print outfile
plot 'full_turn_orbit.dat' u 2:($3*100) i j:j w l
#pause 0.1
reread
}