1 #!/usr/bin/python

     2 

     3 from scipy import signal

     4 from scipy import linalg

     5 import scipy as sp

     6 

     7 u = sp.array([4,5,6])

     8 Eu = linalg.norm(u)**2

     9 ruu = signal.correlate(u,u)

    10 Eu = ruu[len(ruu)/2]