clf;
m = 1;
range = 10;

w = 0.:0.02:m+range;

plot(w,N(m,m,w)-N(m,m+2,w)+N(m,m+4,w)-N(m,m+6,w)+N(m,m+8,w));

break;

hold on;
plot(w,N(m,m,w));

plot(w,N(m,m+2,w));
plot(w,N(m,m-1,w),'-.');
plot(w,N(m,m+1,w),'--');

