Problem:

Find the expectation of $1.7 x^2-1.3 x-3.8$ where $x \sim $ \text{BetaDistribution}[0.4,0.4]

from sympy import *
from sympy.stats import *

x = Beta('x', 0.4, 0.4)
print(E(1.7*x**2-1.3*x-3.8))
