#http://materias.df.uba.ar/estadisticaa2017c1/files/2012/07/Poisson_limite_exacto.py
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()