%《MATLAB及其在理工課程中的應(yīng)用指南》第六章例6-3-1調(diào)用的子程序
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
% 麥克斯韋分子速度分布函數(shù)
function f=mxwl(T,mu,v)
%mu	――分子量,公斤.摩爾-1（如氮為28×10-3）
%v  ――分子速度（可以是一個數(shù)組）
%T  ――氣體的絕對溫度
R=8.31	;			%氣體常數(shù)	
k=1.381*10^(-23);		%玻爾茨曼常數(shù)
NA=6.022*10^23;		%阿伏伽德羅數(shù)
m=mu/NA;					%分子質(zhì)量
%麥克斯韋分布率
f=4*pi*(m/(2*pi*k*T))^(3/2).*exp(-m*v.^2./(2*k*T)).*v.*v;