function [f,g]= rosenb(x);
n=10;
t1=sum( (x(2:n)-x(1:n-1).^2).^2 );
t2=sum( (ones(n-1,1)-x(1:n-1)).^2 );
f=1+100*t1+t2;
if nargout>1;
    g=[-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1));
       200*(x(2:n-1)-x(1:n-2).^2)-400*x(2:n-1).*(x(3:n)-x(2:n-1).^2)-  ..
          2*(ones(n-2,1)-x(2:n-1));200*(x(n)-x(n-1)^2)];
end
