function g = gradf(x)
% autor: Favian Arenas
g = [-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1));
    200*(x(2)-x(1)^2)];
end