clf;
m = 50;
range = 6;

w = m-range:0.02:m+range;

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

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

break;

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

break;

a = w*0;
for j=m-10:2:m+10
if j~=m 
 if rand(1)>0.5
  a = a+N(m,j,w);
 else 
  a = a-N(m,j,w);
 end;
end;
end;
plot (w,a,'.');

