import matplotlib.pyplot as plt
import numpy as np
import scipy.stats as st
import arviz as az
plt.style.use('arviz-darkgrid')
x = np.linspace(0, 4, 500)
for xi in [-1.2, -1.0, -0.8, -0.5, 0.0, 2.0]:
    pdf = st.genpareto.pdf(x, c=xi, loc=0.0, scale=1.0)
    plt.plot(x, pdf, label=rf'$\xi$ = {xi:g}')
plt.xlabel('x', fontsize=12)
plt.ylabel('probability density', fontsize=12)
plt.ylim(0, 2.5)
plt.legend(loc='upper right', title=r'$\mu = 0,\ \sigma = 1$')
plt.show()