NumPy可以将多种数据处理任务表述为简洁的数组表达式，而这种用数组表达式代替循环的做法，通常称为矢量化。

创建二维矩阵

points = np.arange(-5, 5, 0.01) # 1000 equally spaced points
xs, ys = np.meshgrid(points, points)
ys

计算平方根，并绘制图表

from matplotlib.pyplot import imshow, title

import matplotlib.pyplot as plt
z = np.sqrt(xs ** 2 + ys ** 2)
z
plt.imshow(z, cmap=plt.cm.gray); plt.colorbar()
plt.title("Image plot of $\sqrt{x^2 + y^2}$ for a grid of values")
plt.draw()