function [Ls Lq Ws Wq Pd]=queue1(labda,mu,c)
%输入:labda-到达率,mu-服务率,c-服务台数
%系统容量和顾客源都当做无限
%输出:Ls-平均逗留的顾客人数,Lq-平均队长,Ws-平均等待时间,Wq-平均逗留时间
%%
if c==1
rou=labda/mu;
Wq=rou/(mu-labda);
Ws=Wq+1/mu;
Ls=labda*Ws;
Lq=labda*Wq;
Pd=rou;
if Pd>1
Pd=1;
end
elseif c>1
rou=labda/(c*mu);%计算服务强度
P0=(labda/mu)^c/(factorial(c)*(1-rou));
for k=0:(c-1)
P0=P0+(labda/mu)^k/factorial(k);
end
P0=1/P0;
Lq=((c*rou)^c*rou)*P0/(factorial(c)*(1-rou)^2);
Ls=Lq+labda/mu;
Wq=Lq/labda;
Ws=Ls/labda;
Pd=((c*rou)^c*rou)*P0/(factorial(c)*(1-rou));
end