lb=1;
ub=10;
teix=0.01;
[X,Y] = meshgrid(lb:teix:ub,lb:teix:ub);
A=abs(X-5)+abs(Y-5);
Z=sin(A)./A;

% % % mesh
mesh(X,Y,Z)
axis([1 10 1 10 -1 2]);

title('unc#nn#sin1');
xlabel('x');
ylabel('y');
zlabel('z');


% % % % %contour
% 
% v=[1 0.99 0.95   0.5  0.12 0.09   0  ];
% [c,h]=contour(X,Y,Z,v);
% clabel(c,h);
% title('unc#nn#sin1');
% xlabel('x');
% ylabel('y');