%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版)第六章例6－2－4程序
% 完全彈性碰撞問題
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 2007年5月
%
clear
K=logspace(-1,1,11);		% 把K設(shè)成自變量數(shù)組,從K=0.1到K=10按等比取11個點
for i=1:length(K)			% 對各個K循環(huán)計算
   ur1=roots([(1+1/K(i)),-2/K(i),(1/K(i)-1)]);	% 二次方程有兩個解
   ur(i)=ur1(abs(ur1-1)>0.001);		    % 去掉在1鄰近的庸解
end
vr=(1-ur)./K;							% 用（5）式求vr,用元素群運算
em=1-ur.*ur;						    % 主動球損失的相對能量
[K',ur',vr',em']						% 顯示輸出數(shù)據(jù)
semilogx(K',[ur',vr',em']),grid			%	繪圖