Какая вероятность, что из вытащенных шариков: 1 красный, 1 зеленый, 1 синий, 1 желтый, 1 фиолетовый 1 циановый

Table[Module[{sel, pairs, r, g, b, tr, y, v, c},
   sel = RandomChoice[{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
       0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
      1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2}, 30];
   r = 25 - Count[sel, 0];
   g = 15 - Count[sel, 1];
   b = 10 - Count[sel, 2];
   pairs = Partition[sel, 2];
   tr = Map[
      If[#[[1]] != #[[2]], 
        Switch[#, {0, 1}, {3, 3}, {1, 0}, {3, 3}, {0, 2}, {4, 4}, {2, 0}, {4, 4}, 
   {1, 2}, {5, 5}, {2, 1}, {5, 5}], #] &, pairs] //Flatten;
   r = r + Count[tr, 0];
   g = g + Count[tr, 1];
   b = b + Count[tr, 2];
   y = Count[tr, 3];
   v = Count[tr, 4];
   c = Count[tr, 5];
   r*g*b*y*v*c/15890700
   ], 5000000];

Table[Module[{sel, pairs, r, g, b, tr, y, v, c}, 
   sel = RandomChoice[{0, 0, 0, 1, 1, 1, 2, 2, 2}, 6];
   r = 3 - Count[sel, 0];
   g = 3 - Count[sel, 1];
   b = 3 - Count[sel, 2];
   pairs = Partition[sel, 2];
   tr = Map[
      If[#[[1]] != #[[2]], 
        Switch[#, {0, 1}, {3, 3}, {1, 0}, {3, 3}, {0, 2}, {4, 4}, {2, 0}, {4, 4}, 
      {1, 2}, {5, 5}, {2, 1}, {5, 5}], #] &, pairs] // Flatten;
   r = r + Count[tr, 0];
   g = g + Count[tr, 1];
   b = b + Count[tr, 2];
   y = Count[tr, 3];
   v = Count[tr, 4];
   c = Count[tr, 5];
   r*g*b*y*v*c/84], 5000000];

restart;
with(combinat);
d := [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3];
d := [1, 1, 1, 2, 2, 2, 3, 3, 3];
R := nops(d);
tr := 0;
tg := 0;
tb := 0;
for i to R do
    if d[i] = 1 then tr := tr + 1; end if;
    if d[i] = 2 then
        tg := tg + 1;
    end if;
    if d[i] = 3 then
        tb := tb + 1;
    end if;
end do;
print*("R G B=", tr, tg, tb);
W := 0;
V := 0;
P := 3;
N := 5000000;
for cnt to N do
    s := randperm(R);
    r := tr;
    g := tg;
    b := tb;
    y := 0;
    c := 0;
    f := 0;
    for i to P do
        k := 2*i;
        next; if d[s[k]] = d[s[k - 1]];
        if d[s[k]] = 1 then
            r := r - 1;
            if d[s[k - 1]] = 2 then
                g := g - 1;
                y := y + 2;
            else
                b := b - 1;
                f := f + 2;
            end if;
        end if;
        if d[s[k]] = 2 then
            g := g - 1;
            if d[s[k - 1]] = 1 then
                r := r - 1;
                y := y + 2;
            else
                b := b - 1;
                c := c + 2;
            end if;
        end if;
        if d[s[k]] = 3 then
            b := b - 1;
            if d[s[k - 1]] = 1 then
                r := r - 1;
                f := f + 2;
            else
                g := g - 1;
                c := c + 2;
            end if;
        end if;
    end do;
    v := r*g*b*y*c*f;
    W := W + v;
    if 0 < v then
        V := V + 1;
    end if;
end do;
print(1.*V/N);
print(6.!*W/(N*R*(R - 1)*(R - 2)*(R - 3)*(R - 4)*(R - 5)));