function ed = edgetest(img,w) 

ed = [];
vectorx = [];
vectory = [];
for i=10:size(img,1)-10,
  for j=10:size(img,2)-10,
%    [A, vz] = nae(img,i,j,0.05,w);
    [A, vz] = nae(img,i,j,10,w);
    ed(i,j) = A;

    vectorx(i,j)=vz(1,1);
    vectory(i,j)=vz(2,1);
  end
end

clf;
hold on;
imshow(ed);
%quiver(vectory,vectorx);
hold off;
