Problem:

Find the expectation of $-4.1 x^3+0.6 x^2+3.3 x-2.8$ where $x \sim $ \text{BetaDistribution}[1.3,0.6]

from sympy import *
from sympy.stats import *

x = Beta('x', 1.3, 0.6)
print(E(-4.1*x**3+0.6*x**2+3.3*x-2.8))
