function y=f_rosenb(x);

y = 100*(x(2)-x(1).^2).^2+(1-x(1)).^2;