%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第六章例6－3－1程序
% 麥克斯韋氣體分子速度分布率計算
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
% 
clear
T=300;mu=28e-3;	% 給出T,mu
v=eps:1500;		% 給出自變量數組
y=mxwl(T,mu,v); 	% 調用子程序
plot(v,y),hold on % 畫出分布曲線
v1=300:500;			% 給定速度范圍
y1=mxwl(T,mu,v1);% 該范圍的分布
fill([v1,500,300],[y1,0,0],'r')
trapz(y1)		% 求該范圍概率積分
T=200;mu=28e-3;y=mxwl(T,mu,v);plot(v,y)	% 改變T，畫曲線
T=300;mu=2e-3;y=mxwl(T,mu,v);plot(v,y)	% 改變mu, 畫曲線
% gtext('T=200,\mu=28*10^{-3}'),	% 這種標注語句可生成希臘字母及上標（指數）
% gtext('T=300,\mu=28*10^{-3}'),	% 但只能在MATLAB5環境下實現
% gtext('T=300,\mu=2*10^{-3}'),	
hold off