Problem:

Find the expectation of $3.2 x^2-1.8 x-0.4$ where $x \sim $ \text{BetaDistribution}[0.2,1.8]

from sympy import *
from sympy.stats import *

x = Beta('x', 0.2, 1.8)
print(E(3.2*x**2-1.8*x-0.4))
