function [r,err]=ordine(suc,lim)
  %% [r,err]=ordine(suc,lim)
  %% Stima l'ordine di convergenza della successione suc al limite lim.
  err='';
  N=length(suc);
  if nargin==1
    lim=suc(end);
    N=N-1; 
  end
  if N<3
    err='Pochi dati';
    r=[0];
    return;
  end
  for i=1:N
    mm=log(abs(suc(i)-lim));
    if isinf(mm)
      continue
    else
      m(i)=mm;
    end
  end
  if length(m)<=2
    err='Pochi dati';
    r=[0];
    return;
  end
  s=m(2:end)-m(1:end-1);
  r=s(2:end)./s(1:end-1);
