from sympy import stats
X = stats.Normal("x", 66.4, 15) (stats.P(X>73)*0.4+stats.P(X<73)*0.2).evalf()