Problem:

Find the expectation of $-3.3 x^3+0.8 x^2+2.1 x+3.2$ where $x \sim $ \text{BetaDistribution}[1.4,1.1]

from sympy import *
from sympy.stats import *

x = Beta('x', 1.4, 1.1)
print(E(-3.3*x**3+0.8*x**2+2.1*x+3.2))
