Central Limit.py

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from scipy.stats import binom

n = [1,2,10,15,25,50]
p = 0.5
j = 0
fig,ax = plt.subplots(nrows=1, ncols=len(n),figsize=(25,2))

for i in n:
    values = range(0,i+1)
    bi = binom.pmf(values,i,p)
    ax[j].bar(values, bi)
    j += 1