contour+mesh+surf+max&min(14)

f([-10 10],[-10 10],100,3);
function []=f(intervalpox,intervalpoy,korak, crtanje)

x=linspace(intervalpox(1),intervalpox(2),korak);
y=linspace(intervalpoy(1),intervalpoy(2),korak);

[X,Y] = meshgrid(x,y);
Z=sin(X.^2+Y.^2)./(X.^2+Y.^2);

[Zmins, idx] = min(Z);
[~, Zjmin] = min(Zmins);
Zimin = idx(Zjmin);

[Zmaxs, idx] = max(Z);
[~, Zjmax] = max(Zmaxs);
Zimax = idx(Zjmax);

if(crtanje==1)
    mesh(X,Y,Z);
elseif(crtanje==2)
    surf(X,Y,Z);
elseif(crtanje==3)
    contour(X,Y,Z,20);
end
    hold on;grid on;
    plot3(X(Zimax,Zjmax),Y(Zimax,Zjmax),Z(Zimax,Zjmax),'ro');
    plot3(X(Zimin,Zjmin),Y(Zimin,Zjmin),Z(Zimin,Zjmin),'ro');
    xlabel('x');
    ylabel('y');
    zlabel('z');
    title('$z=\frac{\sin(x^2+y^2)}{x^2+y^2}$','Interpreter','latex')
end