

f='2*x.^2-1';  % è la funzione di cui cercare lo zero
df='4*x';
a=-0.5;
b=3;
niter=0; nitermax=100;
toll=0.0001;
x=2.5;
X(1)=x;
diff=toll+1;
while abs(diff)>toll & niter<=nitermax  
  niter=niter+1;
  diff= -eval(f)/eval(df);
  x=x+diff;
  X(niter+1)=x;
end;
if niter<=nitermax  
  zero=x;
  disp(['zero di f è x=' num2str(x)]);
  for i=1:niter-1
    p = 2;
    disp(['p=' num2str(p) '; C=' num2str(abs(X(i+1)-zero)/abs(X(i)-zero)^p) ]);
  end;
end;

% Esercizio: provare con f='3*x.^2-1'  e con f='4*x.^2-1';