quiet
do i = 1, 11
 change [[i]]@beam_start pz @(sqrt((1e6*([[i]]+9))^2-m_electron^2)-ele::0[p0c])/ele::0[p0c]
enddo
set default universe=6
place r13 none
place r23 none
place r33 none
place r14 beta
place r24 dispersion
place r34 all_orbits
place r44 tune
set plot beta autoscale_y=t
set plot dispersion autoscale_y=t
set plot all_orbits autoscale_y=t
set plot tune autoscale_y=t
set dat *@*|invalid = 1
veto dat 1@fit.orbit[2]
veto dat 11@fit.orbit[1]
set dat 1@fit.tune[1]|meas  = 2*cos(0.38*twopi)
set dat 11@fit.tune[2]|meas = 2*cos(0.05*twopi)