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