In [76]:

    

import numpy as np
import scipy.special as sp
import pylab as plt
%matplotlib inline


    

In [77]:

    

N=10
n=np.arange(N)
n1=np.linspace(0,N,N*10)


    

In [78]:

    

p=0.5
q=1-p


    

In [79]:

    

def binominal(n,p,q,N):
    return sp.binom(N,n)*p**(n)*q**(N-n)
print binominal(n,p,q,N)


    

    




[ 0.00097656  0.00976562  0.04394531  0.1171875   0.20507812  0.24609375
  0.20507812  0.1171875   0.04394531  0.00976562]

In [80]:

    

def gauss(n,p,q,N):
    return 1/np.sqrt(2*np.pi*(N*p*q))*np.exp(-0.5*(n-N*p)**2/(N*p*q))
    #return 1/np.sqrt(2*np.pi*(N*p))*np.exp(-0.5*(n-N*p)**2/(N*p))


    

In [81]:

    

plt.plot(n,binominal(n,p,q,N))
plt.plot(n1,gauss(n1,p,q,N))


    

    
Out[81]:





[<matplotlib.lines.Line2D at 0x3aa7510>]






    

In [81]: