function [en]=normu(x,k)

n = length(x);
ns = 0;

for i=1:k:n
    for j=i:i+k-1
        ns = ns + x(j)*x(j);
    end
end

en = ns^(1/2);