$En=       10; 		=>  10 
$M_0=       0.511006;

$sigma_x=         0.5;
$sigma_y=         0.5;
$sigma_xpri=        0.005;
$sigma_ypri=       0.005;  =>        0.01
$Emit_x  = $sigma_x*$sigma_xpri;  =>       0.001
$Emit_y=$sigma_y*$sigma_ypri;  =>       0.001
$Beta_x=($sigma_x )*($sigma_x)/$Emit_x;  =>          10
$Beta_y=($sigma_y)*($sigma_y)/$Emit_y;  =>          10

$I=           6
$B_field=0.0195*$I;
$Grad=$B_field/2.65; 
#$Grad=           0

$Alfa_x=0.0;  =>           0
$Alfa_y=-0.00;  =>           0

#$Emit_x=270e-4 $Emit_y=270e-4 
#$Beta_x=3.333 $Beta_y=3.333

$P=sqrt($En*$En-$M_0*$M_0); 
$T=sqrt($M_0*$M_0+$P*$P)-$M_0; 
$del=          0;  =>           0
$delE=$del*$E;  =>           0

OptiM
#  Insert comment here
Energy[MeV]=$T  Mass[MeV]=$M_0
Emittance: ex[cm]=$Emit_x 	ey[cm]=$Emit_y    DP/P=0.08
Initial:        BetaX[cm]=$Beta_x       BetaY[cm]=$Beta_y
                  AlfaX=$Alfa_x               AlfaY=$Alfa_y
	DispersX[cm]=0 	DispersY[cm]=0 
	Dsp_PrimeX=0 	DspPrimeY=0 
      	X[cm]=0.000     Y[cm]=0.000     Z[cm]=0.000     S[cm]=0.000
      	tetaX[deg]=0            tetaY[deg]=0

begin lattice. Number of periods=1
 o2 Q1 o3 
end lattice
#
begin list
#Q1         	L[cm]=15         	G[kG/cm]=0.0038324 	Tilt[deg]=0
Q1         	L[cm]=15         	G[kG/cm]=$Grad	 	Tilt[deg]=0
#Q2         	L[cm]=15         	G[kG/cm]=0.0688324 	Tilt[deg]=0

 
o2 	L[cm]=41
o3 	L[cm]=56.5

end list of elements