Loc, Scale = ss.norm.fit_loc_scale(unknown)
unknown_cdf = lambda x: ss.norm.cdf(x, loc=Loc, scale=Scale)
plt.plot(np.linspace(-10, 10), unknown_cdf(np.linspace(-10, 10)), '-')
Есть ли способ выполнить арифметику со случайными переменными SciPy?