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