Relativistic Maxwellian distribution generator normalized to the local density By Joan Decker (CEA-DRFC,joan.decker@cea.fr) and Yves Peysson (CEA-DRFC,yves.peysson@cea.fr)
0001 function [fM] = maxwellian_dke_yp(gamma,pn2,dpn,xne_norm,xTe_norm,ir) 0002 % 0003 %Relativistic Maxwellian distribution generator normalized to the local density 0004 % 0005 %By Joan Decker (CEA-DRFC,joan.decker@cea.fr) and Yves Peysson (CEA-DRFC,yves.peysson@cea.fr) 0006 % 0007 if nargin < 6, 0008 error('Wrong number of input arguments in maxwellian_dke_yp') 0009 end 0010 % 0011 fM = exp(-pn2./(1 + gamma)/xTe_norm(ir))/2/pi/xTe_norm(ir)/sqrt(2*pi*xTe_norm(ir)); 0012 normfM = 4*pi*integral_dke_jd(dpn,pn2.*fM,1); 0013 fM = fM*xne_norm(ir)/normfM; 0014