% hold off
step = 0.1;
[x,y] = meshgrid(-5:step:5, -5:step:5);
% x = [-5:step:5];
% y = [-5:step:5];
z = myfun(x,y)

[gx, gy] = gradient(z)


stX = -1
stY = -1

hold on
contour(x,y,z, 128);
quiver(x,y,gx,gy);
plot(1,1,'X');


X = stX, Y=stY
plot(X,Y,'*');

% for i=1:5
%     print 'Help'
%     [gx, gy] = gradient(z, X,Y)
%
%     X = X - gx
%     Y = Y - gy
%
%     plot(X,Y,'*');
%
% end


% hold on
hold off