Implementazione in Matlab (radici multiple)

Quello che segue è il codice per la modifica del metodo di Newton in caso di radici multiple

function [x,i,tolf]=newtonmolt(x0,r,f,df,tolx,nmax)
%NEWTONMOLT Esegue il metodo di Newton modificato per la risoluzione di
%   f(x)=0 anche in presenza di radici multiple, nota la loro
%   molteplicità
%
%   [x,i,tolf]=NEWTONMOLT(x0,r,f,df,tolx,nmax)
%
%   I parametri della funzione sono:
%       x0 -> il punto iniziale
%       r -> molteplicità della radice cercata
%       f -> funzione di cui valutare uno zero
%       df -> la derivata della funzione f
%       toll -> tolleranza per la radice
%       nmax -> limite superiore al numero di iterazioni
%
%   I valori di ritorno sono:
%       x -> la soluzione trovata
%       i -> il numero di iterazioni impiegate per ottenere la soluzione
%       tolf -> la tolleranza sulla funzione
%
%   See Also NEWTON, AITKEN
  i=0;
  err=tolx+1;
  x=x0;
  while(i<nmax & err>tolx)
      fx=feval(f,x);
      dfx=feval(df,x);
      tolf=tolx*abs(dfx);
      if abs(fx)<=tolf
         break
      end
      x1=x-r*(fx/dfx);
      err=abs(x1-x);
      x=x1;
      i=i+1;
  end