function fx=DRosenbrock(x)
fx=[ 400*(x(1)^2-x(2))*x(1)+2*(x(1)-1);
     200*(x(2)-x(1)^2)];
end
