sigmacalcmj3

 on se place dans la limite alfastar<<1 
 constantes :

0001 function sigma=sigmacalcmj3(E,delta,ne,zeff,a,R0,B0,NBOB);
0002 
0003 % on se place dans la limite alfastar<<1
0004 % constantes :
0005 
0006 logcou=17;               % logarithme coulombien
0007 me=9.107e-31;            % masse de l'electron (kg)
0008 qe=1.602e-19;            % charge de l'electron (C)
0009 epsi0=8.854e-12;         % permittivite dans le vide (F.m-1)
0010 c=3e8;                   % vitesse de la lumiere (m.s-1)
0011 
0012 Ej=1.602e-16*E;         %energie en joules
0013 v=c*sqrt(1-1/(1+Ej/(me*c^2))^2);
0014 frecoll=(zeff+1)*ne*logcou*qe^4/(4*pi*epsi0^2*me^2*v^3);
0015 Lalfa=16/pi;
0016 deltastar=2*delta;
0017 sigma2=R0*frecoll*Lalfa/(v*NBOB*sqrt(deltastar));
0018 sigma=sqrt(sigma2);