import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
from scipy.stats import t


if __name__ == "__main__":
    # PDFs for the exponential distribution
    sns.set_palette("deep", desat=.6)
    sns.set_context(rc={"figure.figsize": (8, 4)})

    x = np.linspace(0.0, 5.0, 100)
    lambdas = [0.5, 1.0, 2.0]
    # lambdas = [0.5, 1.0, 1.5, 2.0, 2.5, 4.0]

    for lamb in lambdas:
        y = lamb*np.exp(-lamb*x)
        ax = plt.plot(x, y, label="$\\lambda=%s$" % lamb)
    plt.xlabel("x")
    plt.ylabel("P(x)")
    plt.legend(title="Parameters")
    plt.show()
