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

alpha=0.80;
n=10;
L=(1-alpha)/2
R=L

k=np.arange(0,n+1);
a=beta.ppf(L,k,n-k+1)
b=beta.ppf(1-R,k+1,n-k)

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, n)
plt.ylim(-0.01, 1.01)
plt.show()