Problem:

Find the expectation of $-0.8 x^3-1.8 x^2-4. x+1.5$ where $x \sim $ \text{NormalDistribution}[1.,1.5]

from sympy import *
from sympy.stats import *

x = Normal('x', 1., 1.5)
print(E(-0.8*x**3-1.8*x**2-4.*x+1.5))
