import numpy as np
from matplotlib import pyplot as plt
import scipy.stats as stats
#import seaborn as sb
a = np.linspace(0, 4, 100)
expo = stats.expon
lambda_ = [0.5, 1]
for l, c in zip(lambda_, color):
plt.plot(a, expo.pdf(a, scale=1./l), lw=3,
color=c, label="$\lambda = %.1f$" % l)
plt.fill_between(a, expo.pdf(a, scale=1./l), color=c, alpha=.33)
plt.legend()
plt.ylabel(u"Функция плотности вероятности в точке $z$")
plt.xlabel("$z$")
plt.ylim(0,1.2)
plt.title(u"Функция плотности вероятности экспоненциальной случайной переменной при различных значениях $\\lambda$ \\");