$En= 10; beam energy $M= 0.511; particle mass $gamma= 1+$En/$M; $beta=sqrt(1-1/$gamma/$gamma); =>0.00163231933 $R=25; => 25 $L1=$R*$PI*160/180; => 69.8131701 $E1=$M*$gamma*$beta*$beta/$R*1000; => 1.99999867 $G1=$E1/$R; =>0.0799999467 $L2=10; => 10 $E2=$M*$gamma*$beta*$beta/$L2*$PI*10/180*1000; => 0.872664045 OptiM # Insert comment here Energy[MeV]=$En Mass[MeV]=$M Emittance: ex[cm]=1 ey[cm]=1 DP/P=1e-4. Initial: BetaX[cm]=90.0 BetaY[cm]=90.0 AlfaX=1.35069e-15 AlfaY=-2.0564e-16 DispersX[cm]=122.967 DispersY[cm]=0 Dsp_PrimeX=2.78194e-16 DspPrimeY=0 X=0. Y=0. Z=0. tetaX=0. tetaY=0. begin lattice. Number of periods=1 o1 Q1 o2 Q2 o3 end lattice # begin list Q1 L[cm] = 15 G[kG/cm]=0.04 Tilt[deg]=0 Q2 L[cm]=13 G[kG/cm]=0.04 Tilt[deg]=0 o1 L[cm]=4.0329349 o2 L[cm]=4.0329349 o3 L[cm]=4.0329349 end list of elements