function?out=MMSkteam（skmu1mu2T）
%多服務臺
%s——服務臺個數
%k——最大顧客等待數
%T——時間終止點
%mu1——到達時間間隔服從指數分布
%mu2——服務時間服從指數分布
%事件表：
%???arrive_time——顧客到達事件
%???leave_time——顧客離開事件
%mintime——事件表中的最近事件
%current_time——當前時間
%L——隊長
%tt——時間序列
%LL——隊長序列
%c——顧客到達時間序列
%b——服務開始時間序列
%e——顧客離開時間序列
%a_count——到達顧客數
%b_count——服務顧客數
%e_count——損失顧客數

%初始化

arrive_time=exprnd（mu1）;
leave_time=[];
current_time=0;
L=0;
LL=[L];
tt=[current_time];
c=[];
b=[];
e=[];
a_count=0;
b_count=0;
e_count=0;
%循環
while?min（[arrive_timeleave_time]）????current_time=min（[arrive_timeleave_time]）;
????tt=[ttcurrent_time];????%記錄時間序列
????if?current_time==arrive_time??????????%顧客到達子過程
????????arrive_time=arrive_time+exprnd（mu1）;??%?刷新顧客到達事件
????????a_count=a_count+1;?%累加到達顧客數
????????if??L