Динамическая легенда в цикле

for b=1:1:6
            for K10=100:10:1000
                K1c=K10.*(-log(1-u1)).^(1/b);
                for i=1:1:100000
                    if K1(i) >= K1c(i)
                    n=n+1;
                    end;
                end;
                nn(1+(K10-100)/10)=log(n/1000);
                n=0;
            end
            plot(xx,nn)
            grid on
            title(['b=0.5:0.25:2 ','Sigm=',num2str(sigm),'S=',num2str(S)])
            xlabel('Ln(K10)')
            ylabel('Ln(n)')
            hold on
end

u=rand(100000,1);
u1=rand(100000,1);
razn=1:1:100000;
nn=0:1:90;
n=0;
x=100:10:1000;
xx=log(x);
sigm=40;
b=0.5;
S=20;
for S=20:20:20
    a=-0.14*sqrt(S).*log(1-u);
    for sigm=30:10:30
        figure
        K1=0.7*sigm*sqrt(pi.*a);
        for b=4:1:4
            for K10=100:10:1000
                K1c=K10.*(-log(1-u1)).^(1/b);
                for i=1:1:100000
                    if K1(i) >= K1c(i)
                    n=n+1;
                    end;
                end;
                nn(1+(K10-100)/10)=log(n/1000);
                n=0;
            end
            plot(xx,nn)
            legend('b=1','b=2','b=3','b=4','b=5','b=6')
            grid on
            title(['Sigm=',num2str(sigm),'S=',num2str(S)])
            xlabel('Ln(K10)')
            ylabel('Ln(n)')
            hold on
        end
    end
end

clear, clc
 
t = 0:0.01:pi; % массив времени
f = 1:3; % массив значений часоты
for i = 1:length(f)
    y(:,i) = sin(f(i)*t); % сохраняем в столбец матрицы
    str{i} = ['f_0 = ',num2str(f(i))]; % сохраняем подпись
end
plot(t,y)
legend(str)

t = 0:0.01:pi; % массив времени
f = 1:3; % массив значений частоты
plot(t,sin(f'*t))
str = {strcat('f_0 = ',num2str(f'))}
legend(str)

clear, clc
 
u=rand(100000,1);
u1=rand(100000,1);
razn=1:100000;
 
x = 100:10:1000;
xx = log(x);
 
for S=60
    a = -0.14*sqrt(S).*log(1-u);
    for sigm=30:10:30
        figure
        K1 = 0.7*sigm*sqrt(pi.*a);
        for b = 1:6
            for K10=100:10:1000
                K1c=K10.*(-log(1-u1)).^(1/b);
                n = sum(K1>=K1c);
                nn(b, 1+(K10-100)/10) = log(n/1000);
            end
            str{b} = ['b = ',num2str(b)];
        end
        plot(xx,nn), grid on
        legend(str)
        title(['Sigm=',num2str(sigm),' S=',num2str(S)])
        xlabel('Ln(K10)')
        ylabel('Ln(n)')
    end
end