In [1]:
    
%pylab inline
import numpy as np
n=100
plt.plot(np.random.random(n),np.random.random(n),'ro')
    
    
    Out[1]:
    
In [11]:
    
n=200
plt.plot(np.random.random(n),np.random.random(n),'ro')
    
    Out[11]:
    
In [12]:
    
n=10000
plt.plot(np.random.random(n),np.random.random(n),'ro')
    
    Out[12]:
    
In [8]:
    
def uniform(x,xi=0.5,alpha=1):
    if np.abs(x-xi)<alpha/2.:
        return 1/alpha
    else:
        return 0
    
In [16]:
    
xx=np.linspace(-1,2,200)
plt.plot(xx,map(uniform,xx))
plt.ylim(-0.1,1.1)
plt.grid()
plt.xlabel('$x$',size=20)
plt.ylabel(r'$f(x;\xi,\alpha)$',size=20)
    
    Out[16]:
    
In [5]:
    
x[np.abs(x-0.5)<1./2]
    
    Out[5]:
In [2]:
    
def cube(x): return x*x*x
map(cube, range(1, 11))
    
    Out[2]:
In [ ]: