// residu du'un AR de coefficients phi
function eps = residu_ar(X, phi)
  T = length(X);
  ordre = length(phi);
  eps = X(ordre+1:T);
  for i=1:ordre
    eps = eps - phi(i) * X(ordre-i+1:T-i);
  end
endfunction