clear;
clc;
%M/M/1排隊系統仿真
SimTotal=input（‘請輸入仿真顧客總數SimTotal=‘）;?%仿真顧客總數；
Lambda=0.2;?????%到達率Lambda；
Mu=0.9;?????????%服務率Mu；
?
t_Arrive=zeros（1SimTotal）;??
t_Leave=zeros（1SimTotal）;?
ArriveNum=zeros（1SimTotal）;?
LeaveNum=zeros（1SimTotal）;?

Interval_Arrive=-log（rand（1SimTotal））/Lambda;%到達時間間隔
Interval_Serve=-log（rand（1SimTotal））/Mu;%服務時間
t_Arrive（1）=Interval_Arrive（1）;%顧客到達時間
?
ArriveNum（1）=1;?
for?i=2:SimTotal??
????t_Arrive（i）=t_Arrive（i-1）+Interval_Arrive（i）;?
????ArriveNum（i）=i;?
end?

t_Leave（1）=t_Arrive（1）+Interval_Serve（1）;%顧客離開時間
?
LeaveNum（1）=1;?
for?i=2:SimTotal??
????if?t_Leave（i-1）????????t_Leave（i）=t_Arrive（i）+Interval_Serve（i）;?
????else??
????????t_Leave（i）=t_Leave（i-1）+Interval_Serve（i）;?
????end??
????LeaveNum（i）=i;?
end?