
k=-(n/2):((n/2)-1); 
k=k';  dx = 2*pi/n;

ltk = k.^2;
lk  = 2*(1-cos(k*dx))/(dx*dx);

plot(k,ltk,'ro',k,lk,'bo')

