1

M = [286588 434446 388251];
K = [0.1452 0.8242 0.7653];
L = 997783;
temp = [];
n = [100 600 1000];
for i = 1:3
    [M,K,L,temp] = myRand(M,K,L,n(i));
    disp(['Среднее значение n=' int2str(n(i)) ': ' num2str(mean(temp))])
    disp(['Дисперсия n=' int2str(n(i)) ': ' num2str(cov(temp))])
    disp(['Стандартное отклонение выборки n=' int2str(n(i)) ': ' num2str(std(temp))])
end
function [M, K, L, R] = myRand(M, K, L, x)
    R = 1:x;
    for i = 1:x
        n = length(K)+1;
        k = (M(1)*K(n-1)+M(2)*K(n-2)+M(3)*K(n-3))-floor((M(1)*K(n-1)+M(2)*K(n-2)+M(3)*K(n-3))/L)*L;
        K = [K k];
        R(i) = k/L;
 
    end
figure(i)
plot(R)
title(['СВ N = ', num2str(i)])
ylabel('Значения СВ')
xlabel('Индекс СВ')

Y = sort(R);
figure(i+3)
plot(Y)
title(['Отсортированные СВ N = ', num2str(i)])
ylabel('Значения СВ')
xlabel('Индекс СВ')
end