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