function z = method(x, y, e, n)

        n=n-1;
      z =  func2min(x,y);

        [xn, yn] = evk_norm(x, y, e);


        if func2min(xn,yn) < func2min(x,y)
           % fprintf('меньше\n')
            x = xn;
            y = yn;
           % z = func2min(x,y);
            z=method(x, y, e, 10);
            return
        end

   while n>0
       % fprintf('n=%d\n',n)
       z= method(x, y, e, n);

       return
   end

   fprintf('z=%.5f\nx=%.5f\ny=%.5f\n',z,x,y)
   return
end

function [x1, y1] = evk_norm(x, y, e)
    x1 = x + (rand(1)-0.5)*2*e;
    y1 = y + (rand(1)-0.5)*2*e;
   % fprintf('x1=%.5f y1=%.5f\n',x1,y1 )
    if sqrt((x1-x)^2 + (y1-y)^2) > e
        evk_norm(x,y,e);
    end
end