function y=QQ(L,r,k,UCL,LCL)
d=(UCL-LCL)/k
m1=0;
for i=-k/2+1:1:k/2
m1=m1+1;
k1=1;
m=k/2;
for j=-m+1:1:m
temp=((j*d)-(1-r)*(i-0.5)*d)/(r);
T1=normcdf(temp,0,1);
temp=((j-1)*d-(1-r)*(i-0.5)*d)/(r);
T2=normcdf(temp,0,1);
R(m1,k1)=T1-T2;
k1=k1+1;
end;
end;
PM=0*ones(1,2*m);
PM(k/2)=1;
SS=PM*inv(eye(2*m)-R)*ones(2*m,1);
y=SS;

我定义了function y=QQ(L,r,k,UCL,LCL)这个函数，但是运行时最终结果和L的值无关，这是用MATLAB解决马尔科夫链转移概率问题，由于对MATLAB了解较少，所以一直没搞清楚，希望有这方面的高手给指点一下，我这个程序问题出在哪了？万分感谢！！！！