function AlgoNewton

function[x]=AlgoNewton(f,x)
h=0.0001;
d2f=hessf(f,x);
  if (det(d2f)==0)
      disp('la méthode de newton n''est pas applicable pour les matrices non inversibles');
  else
      df=gradf(f,x);
      while ((norm(df))^2>h)
        df=gradf(f,x);  
        d2f=hessf(f,x);
        d=-(inv(d2f))*df;
        x=x+d;
      end
      disp('la solution optimale du problème est: ');
      x;
      disp('la valeur minimale que peut donc atteindre f est: ');
      f(x)
  end