Исправить код

% Ввод исходных данных
lambda=14.5; mu=6.4;
m=5; n=8;
N=m+n+1; % число возможных состояний системы
% Определение вспомогательных переменных
Num=input('Число изменений состояния системы за время модели-рования');
s=1; % индекс текущего состояния системы
T=zeros(1,Num);% вектор моментов изменений состояний
Qt=zeros(1,Num);% состояния процесса Q(t) в моменты
% изменений состояний системы
Time=zeros(1,N);% вектор, в котором накапливаются
% суммарные времена пребывания системы
% в состояниях 0, 1, ..., m+n
for i=2:Num
    switch Qt (i-1)
        case 0
            v=Inf; u=-log(rand)/lambda;
        case m+n
            u=Inf; v=-log(rand)/(m*mu);
        otherwise
            u=-log(rand)/lambda;
            qs=min(m,Qt(i-1));% число работающих приборов
            v=-log(rand)/(qs*mu);
    end;
    delta_time=min(u,v);
    T(i)=T(i-1)+delta_time;
    Time(s)=Time(s)+delta_time;
    if u<v
        Qt(i)=Qt(i-1)+1;
    else
        Qt(i)=Qt(i-1)-1;
    end;
    s=Qt(i)+1;% индекс текущего состояния
end;
z=zeros(1,Num);
mvect=ones(1,Num)*m;
qt=max(z,Qt-mvect);% очередь за время моделирования
% Построение на одном экране графика состояний системы и
% числа клиентов в очереди
subplot(211);plot(T,Qt);
subplot(212);plot(T,qt);
disp('Оценки характеристик системы MIM1518,');
disp('про лученные на основе моделирования');
Psost=Time/T(Num);
disp('Вероятность простоя'); disp(Psost(1));
disp('Вероятности состояний 1, ..., m+n');
disp(Psost(2:N));
disp('Вероятность ожидания'); Pw=1-sum(Psost(1:m));
disp(Pw);
disp('Вероятность отказа'); PL=Psost(N); disp(PL);
disp('Среднее число заявок в системе');
k=0:(N-1); Qsr=k*Psost'; disp(Qsr);
disp('Среднее число заявок в очереди');
z=zeros(1,N); qsr=max(z,k-m)*Psost'; disp(qsr);