正确运行程序并获 q-T 图像

 ###### 问题遇到的现象和发生背景

clear ;
clc;

syms T0 toa k T_e n v u b M_T m_T mu  h_ij r z_j f;
f=zeros(10,10);
h_ij=zeros(1,1);
h2=zeros(10,10);
T0=1;
toa=0.2;
k=0;
T_e=10;
n=10;
v=0;
b=0.79;

q=0;
p=0;
h=0;
r=0.2;
z_j=0.1;
e=0;
a=0;
T(1, 1) =0;
M_T(1,1)= 0;
m_T(1,1)= 0;
i=0;

while T<T_e
    T(a + 2, 1) = T(a + 1, 1) + toa;
    a=a+1 ;
    k=0;
    while k<n
        z=1+(2-1).*rand([a,k]);
        sort(z);
        k=k+1;
    end
    
    while 1.0+r<=2
        for i=2:1:9
            e=e+1;
            if z(a,(i-1))<1.2+r && z(a,(i-1))>1.0+r
           h=h+1;
           h2(a,(i-1))=h/e*0.1;
           f(a,(i-1))=h2(a,(i-1))./(z(a,(i-1))-z(a,i));
            end
            r=r+z_j;
        end
    end
    v=v+1;
    k=0;
    z_max=max(max(z));
    z_min=min(min(z));
    u=z_max-z_min;
    if u<=b
        mu=1;
        m_T(i + 2, 1)=m_T(i + 1, 1)+ mu;
        M_T(i + 2, 1)=M_T(i + 1, 1) + 1;
       
    end
    if v<M_T
        
    else
       q=(1-(m_T/M_T)) ;
       p=((m_T)/(M_T));
      
      plot(T,q,'bo') ;
      
        hold on;
        grid on;
    end
end

 ###### 运行结果及报错内容
错误使用 sub2ind (第 43 行)
下标超出范围。

出错 sym/subsref (第 904 行)
R_tilde = sub2ind(size(L), Idx.subs{:});

出错 polika (第 56 行)
m_T(i + 2, 1)=m_T(i + 1, 1)+ mu;

 ###### 我想要达到的结果

程序正常运行后得到图像