Problem:

Find the expectation of $4.6 x^2+4.9 x+2.7$ where $x \sim $ \text{BetaDistribution}[0.4,1.4]

from sympy import *
from sympy.stats import *

x = Beta('x', 0.4, 1.4)
print(E(4.6*x**2+4.9*x+2.7))
