Problem:

Find the expectation of $3.3 x^2-0.4 x+3.4$ where $x \sim $ \text{PoissonDistribution}[3.9]

from sympy import *
from sympy.stats import *

x = Poisson('x', 3.9)
print(E(3.3*x**2-0.4*x+3.4))
