def gauss2D(x, mean, cov): # x, mean, cov??numpy.array?? z = -np.dot(np.dot((x-mean).T,inv(cov)),(x-mean))/2.0 temp = pow(sqrt(2.0*pi),len(x))*sqrt(det(cov)) return (1.0/temp)*exp(z)