% Jacobian of  Ffun function
  function J = Jfun(x)   
  J = zeros(2,2) ;
  J(1,1) = 2*x(1) ; J(1,2) = 2*x(2) ;
  J(2,1) = 2*(x(1)-1) ; J(2,2) = 2*x(2) ;
 
 
