Помогите переделать учебный скрипт рисующий поверхность, так чтобы он рисовал поверхность, заданную координатами площадок, а не функцией. Координаты площадок храняться в текстовом файле, например в data.txt:

0 1 0
1 0 0
0 -1 0
-1 0 0 

...

в интернете популярен пример

#!/usr/bin/python
import pylab
from mpl_toolkits.mplot3d import Axes3D
import numpy

def makeData ():
    x = numpy.arange (-10, 10, 0.1)
    y = numpy.arange (-10, 10, 0.1)
    xgrid, ygrid = numpy.meshgrid(x, y)

    zgrid = numpy.sin (xgrid) * numpy.sin (ygrid) / (xgrid * ygrid)
    return xgrid, ygrid, zgrid

x, y, z = makeData()

fig = pylab.figure()
axes = Axes3D(fig)

axes.plot_surface(x, y, z)

pylab.show()

не знаю логики интерполяции поверхности (проволочная поверхность не нужна) и проблемы с синтаксисом чтения из файла

P.S. порекомендуйте наиболее полный учебник по matplotlib и матричным вычислениям в Питоне