The example tracks the equilibrium beam sizes computed in section 6.2. The routine
extract_beam_sizes(xs,nf%s_ij0,sizes) \end{verbatim}
uses \eq{benv} where
%
%]|Expr|[#b @`b___}(#" Helvetica|: ;bP8&c0!*,D$^"!Symbol^:!&c0 S|
%|_^""*|:"&c0!*0: ,D}]|[
${\Sigma}^{0}$
is {\tt nf\%s\_ij0},
%
%]|Expr|[#b @`b___})## b'4" Helvetica|: ;bP8&c0!*,D$^"!Symbol^:!&c0 S|
%|_^""*|:"&c0!*f: ,D}& b!( b"0 b#8 b$@ b%H b&P!WW}]|[
${\Sigma}^{f}$ is {\tt sizes} and the $E$ is contained in the probe\_8 {\tt xs}.
\begin{verbatim}
!!!!!!!!!! Tracking a probe_8 to beam sizes around the ring !!!!!!!!!!
call kanalnummer(mf,"beam_xx_inf.dat")
state=default0+radiation0+envelope0
xs0=x ! closed orbit intializing a probe
m=1 ! damapspin set to identity
xs=xs0+m ! Probe_8 = closed orbit probe + Identity
p=>als%start
call extract_beam_sizes(xs,nf%s_ij0,sizes)
write(mf,*) p%mag%name, sizes(1,1)
do i=1,als%n
call track_probe(als,xs,state,fibre1=i,fibre2=i+1)
call extract_beam_sizes(xs,nf%s_ij0,sizes)
p=>p%next
write(mf,*) p%mag%name, sizes(1,1)
enddo
close(mf)
!!!!!!!!!! Tracking a probe_8 to beam sizes around the ring using emittances !!!!!!!!!!