Problem:

Find the expectation of $0.4 x^3-2.5 x^2+1.8 x-1.5$ where $x \sim $ \text{ExponentialDistribution}[0.4]

from sympy import *
from sympy.stats import *

x = Exponential('x', 0.4)
print(E(0.4*x**3-2.5*x**2+1.8*x-1.5))
