function x=F(u,i)
i=max(1,i);
n=numel(u);
i=min(n,i);
x=0.5*u(i).*u(i);