%《MATLAB及其在理工課程中的應(yīng)用指南》第六章例6－3－1程序
% 麥克斯韋氣體分子速度分布率計算
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
% 氣體分子速度的分布率
clear
T=300;mu=28e-3;	                            % 給出T,mu
v=eps:1500;		                            % 給出自變量數(shù)組
y=mxwl(T,mu,v); 	                        % 調(diào)用子程序
plot(v,y),hold on                           % 畫出分布曲線
v1=300:500;			                        % 給定速度范圍
y1=mxwl(T,mu,v1);                           % 該范圍的分布
fill([v1,500,300],[y1,0,0],'r')
set(gcf,'color','w')                        % 置圖形背景色為白色
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}'),	            % 這種標注語句可生成希臘字母及上標（指數(shù)）
gtext('T=300,\mu=28*10^{-3}'),	            % 但只能在MATLAB5環(huán)境下實現(xiàn)
gtext('T=300,\mu=2*10^{-3}'),	
hold off