matlab 结果>=0,求助！Matlab中为什么结果输出的都是0啊！

function[]=msququeright(m,n,r,u)

outfile=fopen('MMmn.txt','at+');

fprintf(outfile, '排队系统\n');

r=10;u=8;

fprintf(outfile, '\n平均到达时间间隔%12.3f minutes',r);

fprintf(outfile, '\n平均服务时间%12.3f minutes',u);

global s1;

global s2;

global P0;

global Pk;

global Lq;

global Ls;

global Wq;

global Ws;

for m=1:3

for n=m:10

fprintf(outfile, '\n装载机数量%12d', m);

fprintf(outfile, '\n自卸汽车数量%14d',n);

for i=0:m

s1=sum(1/factorial(i)*1/factorial(n-i)*(u/r)^i);

end

for j=m+1:n

s2=m^m/factorial(m)*sum(1/factorial(n-j)*(u/m*r)^j);

end

P0=1/factorial(n)*(s1+s2)^-1;

fprintf(outfile, '\n装载机空闲概率%12.3f',P0);

for  k=0:n

if k<=m

f=factorial(n)/factorial(n-k);

g=f/factorial(k);

Pk=g*(u/r)^k*P0;

else

Pk=factorial(n)/factorial(n-k)/factorial(m)/m^(k-m)*(u/r)^k*P0;

Lq=sum((k-m)*Pk);

Ls=sum(k*Pk);

end

fprintf(outfile, '\n系统中有k辆车的概率%12.3f ',Pk);

end

fprintf(outfile, '\n平均排队长%12.3f ',Lq);

fprintf(outfile, '\n平均队长%12.3f ',Ls);

fprintf(outfile, '\n平均等待时间%12.3f minutes',Wq);

fprintf(outfile, '\n平均逗留时间%12.3f minutes\n',Ws);

end

end

fclose(outfile);

%%%%%%%%%%%%%

再就是输出系统中有k辆车的概率的时候这个k能改变吗？

[本帖最后由 rankfirst 于 2009-6-14 21:52 编辑]