function g = rosgrd(x, data)

% y = 100*(x(1)^2-x(2))^2+(1-x(1))^2;
g = zeros(2, 1);
g(1) =  400*x(1)*(x(1)^2-x(2))-2*(1-x(1));
g(2) = -200*(x(1)^2-x(2));