
function temperature_gain = fusion_of_T(T,N)
n=1E20; % m^-3 number density of plasma
E_G=1.18; % MeV Gammow of Energy
Q=17.5; % MeV Energy per D-T->He+n fusion
S=14; % MeV Nuclear factor 
k_b=8.617E-4; % MeV/1E7K
c=3E8;%m/s
mass_dt=2.5*938.7; % equal D-T avg 2.5*mass of proton in MeV/c^2
v=sqrt(3*k_b*T/(mass_dt))*c;
sigma=S*E_G^(-1/3)*(k_b*T)^(-2/3)*N^(-3)*exp(-(E_G/(k_b*T))); % barn


sigma_v = sigma*v; %approx D-T cross section
rate = n^2 * sigma_v * (1E-28); % 1/(s*m^3)
temperature_gain = rate*Q/k_b; % K/(s*m^3)
