sincos[k_, ls_, fi_] := Sin[fi] + (Sqrt[k]) ls Cos[fi]

intinter[fi_, k_, ls_] := (∫_0^fi) 
        sincos[k, ls, x] \[DifferentialD]x

intest[k_, l_, 
        ls_] := (∫_0^((Sqrt[k]) l)) ((Abs[
                    sincos[k, ls, y]])^3 (intinter[y, k, 
                    ls])^2) \[DifferentialD]y

intest[0.1158931845985401, 2.74, 3.5]

re = 2.817940325 10^(-15)

lame = 3.861592678 10^(-13)

emi =  0.2 10^-7

sigma min

(7/
            5)^0.5 (275/3/(6 Pi)^0.5 re lame intest[
              0.1158931845985401, 2.74, 3.5]) ^(1/7) (emi)^(5/
            7)

emi1 =  0.1 10^-7

(7/
            5)^0.5 (275/3/(6 Pi)^0.5 re lame intest[
              0.1158931845985401, 2.74, 3.5]) ^(1/7) (emi1)^(5/
            7)

sigma

betas = 17.92472388 10^(-6)

17.92472388`

(betas emi 0.000511/
              1500 + 110/
                3/(6 Pi)^0.5 re lame  (1500/
                  0.000511)^5 intest[0.1158931845985401, 2.74, 
              3.5] (emi 0.000511/1500/betas )^(5/
                  2))^0.5

cambio l

(betas emi 0.000511/
              1500 + 110/
                3/(6 Pi)^0.5 re lame  (1500/
                  0.000511)^5 intest[0.1158931845985401, 3, 
              3.5] (emi 0.000511/1500/betas )^(5/
                  2))^0.5

(betas emi 0.000511/
              1500 + 110/
                3/(6 Pi)^0.5 re lame  (1500/
                  0.000511)^5 intest[0.1158931845985401, 1.37, 
              3.5] (emi 0.000511/1500/betas )^(5/
                  2))^0.5

Plot[(betas emi 0.000511/
                1500 + 110/
                  3/(6 Pi)^0.5 re lame  (1500/
                    0.000511)^5 intest[0.1158931845985401, l, 
                3.5] (emi 0.000511/1500/betas )^(5/
                    2))^0.5, {l, 1, 5}]

betay*

(275/3/(6 Pi)^0.5 re lame intest[0.1158931845985401, 
              2.74, 3.5]) ^(2/7)  1500/
        0.000511 (emi)^(3/7)

sigma* 
      vs betay*  table plot and creazione del  . 
        dat  (CLIC - BDS  l*  = 3.5 )

tdati = 
      Table[{b, (b 10^(-6) emi 0.000511/
                      1500 + 110/
                        3/(6 Pi)^0.5 re lame  (1500/
                          0.000511)^5 intest[0.1158931845985401, 2.74, 
                      3.5] (emi 0.000511/1500/b/ 
                            10^(-6))^(5/
                          2))^0.5 10^(9)}, {b, 5, 190, 
          2}]

ListPlot[tdati, PlotJoined -> True,  
      PlotRange -> {0, 12}]

Export["OideEffect.dat", tdati]