función sympy.stats.Skellam() en Python

Con la ayuda del método sympy.stats.Skellam() , podemos crear una variable aleatoria discreta con una distribución de Skellam. 

El Skellam es la distribución de la diferencia N1 – N2 de dos variables aleatorias estadísticamente independientes N1 y N2, cada una con distribución de Poisson con los respectivos valores esperados mu1 y mu2.

Syntax:  sympy.stats.Skellam(name, mu1, mu2)

Parameters:
 mu1:  A non-negative value
 mu2:  A non-negative value


Returns: discrete random variable with a Skellam distribution.

Ejemplo 1 :

# import sympy, Skellam, density, Symbol
from sympy.stats import Skellam, density
from sympy import Symbol
  
mu1 = Symbol("mu1", positive = True)
mu2 = Symbol("mu2", positive = True)
  
# using sympy.stats.Skellam() method
X = Skellam("x", mu1, mu2)
skeDist = density(X)(z)
  
print(skeDist)

Producción:

(mu1/mu2)**(z/2)*exp(-mu1 - mu2)*besseli(z, 2*sqrt(mu1)*sqrt(mu2))

Ejemplo #2:

# import sympy, Skellam, density
from sympy.stats import Skellam, density
  
# using sympy.stats.Skellam() method
X = Skellam("x", 1, 2)
skeDist = density(X)(3)
  
print(skeDist)

Producción:

sqrt(2)*exp(-3)*besseli(3, 2*sqrt(2))/4