bmaxcalcmj3

 E en keV 
 on se place dans la limite alfastar<<1
 on neglige la variation de la transformee rotationnelle

0001 % E en keV
0002 % on se place dans la limite alfastar<<1
0003 % on neglige la variation de la transformee rotationnelle
0004 
0005 function deltaBmax=bmaxcalcmj3(E,delta,a,R0,B0,NBOB);
0006 
0007 me=9.107e-31;
0008 qe=1.602e-19;
0009 c=3e8;
0010 
0011 Ej=1.602e-16*E;
0012 v=c*sqrt(1-1/(1+Ej/(me*c^2))^2);
0013 
0014 Talfa=12;
0015 lambdadelta=R0/(NBOB);
0016 const=Talfa*me/(2*qe*B0*NBOB*lambdadelta);
0017 deltaBmax=const*v*sqrt(delta);