# coding=utf-8
import numpy as np
from scipy import stats
import matplotlib.pyplot as plt
n = 10
p = 0.3
k = np.arange(0,21)
binomial = stats.binom.pmf(k, n, p)
plt.plot(k, binomial, 'o-')
plt.title('Binomial: n=%i, p=%.2f' % (n, p), fontsize=15)
plt.xlabel('Number of successes')
plt.ylabel('Probability of successes', fontsize=15)
plt.show()