#
# pylabimage.py
# (c) Neil Gershenfeld 2/10/08
# plot sin(r)/r with pylab
#

from pylab import *
from numpy import *

x = arange(-20.5,20.7,.2)
y = arange(-20.5,20.7,.2)
(X,Y) = meshgrid(x,y)
R = sqrt(X*X + Y*Y)
Z = sin(R)/R
gray()
imshow(Z)
show()