Compare fusion cross-section methods
%-------------------------------------------------------------------------- % See also MaxwellianV, SigmaV %-------------------------------------------------------------------------- %-------------------------------------------------------------------------- % Copyright (c) 2017 Princeton Satellite Systems, Inc. % All rights reserved. %-------------------------------------------------------------------------- tK = logspace(0,2.17); maxwV = MaxwellianV( {'d' 'd' 'd'}, {'t' 'he3' 'd'}, tK ); sigmaV = SigmaV( {'d' 'd' 'd'}, {'t' 'he3' 'd'}, tK ); Plot2D(tK, maxwV,'Temperature (keV)','Cross-Section','Comparison','ylog'); hold on; plot(tK,sigmaV,'--') plot([10 20 50 100],[2.3e-19 3.8e-18 5.4e-17 1.6e-16]*1e-6,'*') legend('d-t','d-3he','d-dp','d-t','d-3He','d-dp','NRL') %--------------------------------------