import numpy as np ; import pylab as pl
from scipy.special import jn
x = np.linspace(-5, 15, 100)
for i in range(10):
y = jn(i, x)
pl.plot(x, y, label='$j_%i$' % i)
pl.title('Fonctions de Bessel')
pl.legend()
import scipy as sp ; import pylab as pl
l = sp.lena()
pl.imshow(l, cmap=pl.cm.gray)
pl.axis('off')
import numpy as np ; import pylab as pl
x, y, value = np.random.normal(size=(3, 50))
pl.scatter(x, y, np.abs(50*value), c=value)
import numpy as np ; import pylab as pl
x, y = np.mgrid[-5:5, -5:5]
u = -x
v = y
pl.quiver(x, y, u, v)