#!/usr/bin/python from scipy import signal from scipy import linalg import scipy as sp u = sp.array([4,5,6]) Eu = linalg.norm(u)**2 ruu = signal.correlate(u,u) Eu = ruu[len(ruu)/2]