Получить поток из 100k сообщений и для сообщения каждого типа определить

k = 5;

t1 = 0; t2 = 0; t3 = 0;

a1t1 = 0; a2t1 = 0; a3t1 = 0; a4t1 = 0;
a1t2 = 0; a2t2 = 0; a3t2 = 0; a4t2 = 0;
a1t3 = 0; a2t3 = 0; a3t3 = 0; a4t3 = 0;

for i = 1:1:100*k
  
  r = grand(1, 1, "unf", 0, 1);
  msg = 0;
  
  if r >= 0 & r < 0.49 then
    t1 = t1 + 1;
    msg = 1;
  elseif r >= 0.49 & r < 0.89 then
    t2 = t2 + 1;
    msg = 2;
  else
    t3 = t3 + 1;
    msg = 3;
  end
  
  r = grand(1, 1, "unf", 0, 1);
  
  select msg,
    case 1 then
	  if r >= 0 & r < 0.05 then
	    a1t1 = a1t1 + 1;
	  elseif r >= 0.05 & r < 0.37 then
	    a2t1 = a2t1 + 1;
	  elseif r >= 0.37 & r < 0.76 then
	    a3t1 = a3t1 + 1;
	  else
	    a4t1 = a4t1 + 1;
	  end,
    case 2 then
	  if r >= 0 & r < 0.41 then
	    a1t2 = a1t2 + 1;
	  elseif r >= 0.41 & r < 0.7 then
	    a2t2 = a2t2 + 1;
	  elseif r >= 0.7 & r < 0.87 then
	    a3t2 = a3t2 + 1;
	  else
	    a4t2 = a4t2 + 1;
	  end,
    case 3 then
	  if r >= 0 & r < 0.3 then
	    a1t3 = a1t3 + 1;
	  elseif r >= 0.3 & r < 0.37 then
	    a2t3 = a2t3 + 1;
	  elseif r >= 0.37 & r < 0.61 then
	    a3t3 = a3t3 + 1;
	  else
	    a4t3 = a4t3 + 1;
          end
    end
end

disp(t1, t2, t3);
disp(t1/(100*k), t2/(100*k), t3/(100*k));

disp(a1t1, a2t1, a3t1, a4t1);
disp(a1t1 / t1, a2t1 / t1, a3t1 / t1, a4t1 / t1);

disp(a1t2, a2t2, a3t2, a4t2);
disp(a1t2 / t2, a2t2 / t2, a3t2 / t2, a4t2 / t2);

disp(a1t3, a2t3, a3t3, a4t3);
disp(a1t3 / t3, a2t3 / t3, a3t3 / t3, a4t3 / t3);