import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import chi2

alpha=0.68;

L=(1-alpha)/2
R=L

k=np.arange(0,12);
a=chi2.ppf(L,2*k)/2
b=chi2.ppf(1-R,2*(k+1))/2

print(a)
print(b)

plt.step(k+1, a, label='Límite inferior')
plt.step(k, b, label='Límite superior')
plt.title("Cinturón de confianza")
plt.legend()
plt.grid(True)
plt.xlim(-0.1, 10)
plt.ylim(-0.1, 14)
plt.show()