df = @(x) [6*x(1) 2*x(2)];
x0 = [1 1];
eps=0.0 1;

function x_koz = GradiensModszer(df, x0, eps)
 % x_k+1 = x_k - df(x_k)*alfa
  x_koz = x0;
  alfa = 0.1;
  
  norma = 2*eps;
 
  while (norma > eps)
    x_uj = x_koz - df(x_koz)*alfa
    norma = norm(x_uj - x_koz);
    x_koz = x_uj;
  end
 end

 GradiensModszer(df, x0, eps);