Com a ajuda do método sympy.stats.Multinomial() , podemos criar uma variável aleatória discreta com Distribuição Multinomial.

Uma distribuição multinomial é a distribuição de probabilidade dos resultados de um experimento multinomial.

Sintaxe: sympy.stats.Multinomial (syms, n, p)
Parâmetros: 
syms: o símbolo
  n: é o número de tentativas, um inteiro positivo
  p: probabilidades de evento, p> = 0 e p <= 1
Retorna: uma variável aleatória discreta com distribuição multinomial

Exemplo 1 :

from sympy.stats.joint_rv_types import Multinomial 
from sympy.stats import density 
from sympy import symbols, pprint 
  
x1, x2, x3 = symbols('x1, x2, x3', nonnegative = True, integer = True) 
p1, p2, p3 = symbols('p1, p2, p3', positive = True) 
M = Multinomial('M', 3, p1, p2, p3) 
multiDist = density(M)(x1, x2, x3) 
  
pprint(multiDist)

Resultado :

/ x1 x2 x3
| 6 * p1 * p2 * p3
| ---------------- para x1 + x2 + x3 = 3
<x1! * x2! * x3!
|
| 0 caso contrário
\

Exemplo # 2:

from sympy.stats.joint_rv_types import Multinomial 
from sympy.stats import density 
from sympy import symbols, pprint 
  
x1, x2, x3 = symbols('x1, x2, x3', nonnegative = True, integer = True) 
M = Multinomial('M', 4, 0, 1, 0) 
multiDist = density(M)(x1, x2, x3) 
  
pprint(multiDist)

Resultado :

/ x1 x3
| 24 * 0 * 0
| ----------- para x1 + x2 + x3 = 4
<x1! * x2! * x3!
|
| 0 caso contrário
\