from pylab import * def gau(x): y = exp(-pi*(x**2)) y = array(y) return y