Problem:

Find the expectation of $1.1 x^3-2.3 x^2+1.5 x-1.6$ where $x \sim $ \text{PoissonDistribution}[2.]

from sympy import *
from sympy.stats import *

x = Poisson('x', 2.)
print(E(1.1*x**3-2.3*x**2+1.5*x-1.6))
