Problem:
Find the probability that $1.25 < 1.7 x^2+2.5 x+0.2 < 2.08$ where $x \sim $ \text{ExponentialDistribution}[1.6].
from sympy import *
from sympy.stats import *

x = Exponential('x', 1.6)
print(P((1.25 < 1.7*x**2+2.5*x+0.2) & (1.7*x**2+2.5*x+0.2 < 2.08)))
