#!/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]