from sympy import *
init_printing(use_latex='mathjax')
x, y, z = symbols('x,y,z')
from sympy.stats import *
X = Normal('X', 0, 1)
X
# Density
# Density of expressions
# Probability
# Sampling
A = Die('A', 6)
B = Die('B', 6)
density(A+B)
from sympy.stats import Gamma
Gamma?
k = Symbol('k', positive=True)
theta = Symbol('theta', positive=True)
A = Gamma('A', k, theta)
A
density(A**2)(x)