function obj=ros(x,y)
  obj=(1-x).^2+100*(y-x.^2).^2;
end