Calcul de la fraction restant dans les miroirs en fonction de l'energie : E = energie initiale en keV de la particule piegee
0001 % Calcul de la fraction restant dans les miroirs en fonction de l'energie : 0002 % E = energie initiale en keV de la particule piegee 0003 0004 function f=fracmircalcmj3(E,delta,ne,zeff,a,R0,B0,NBOB); 0005 0006 n=5000; 0007 rho1=rhocalcmj3(E,delta,ne,zeff,a,R0,B0,NBOB); 0008 nume=(rho1/2)*(1+erf(sqrt(n/2)*rho1))+(1/sqrt(2*pi*n))*exp((-1)*n*rho1^2/2); 0009 deno=(rho1/2)*(1+erf(sqrt(1/2)*rho1))+(1/sqrt(2*pi))*exp((-1)*rho1^2/2); 0010 f=nume/deno;