function c=a_to_cepst(a,NBC);
% Syntax: c=a_to_cepst(a)
%      ou c=a_to_cepst(a,NBC)
%
% a est vecteur colonne de la dimension N+1: [1 a(1) a(2) ... a(N)]'
% NBC est le nombre de coefficients cepstraux desires (default N)
% c sera vecteur colonne de la dimension NBC+1: [0 c(1) c(2) ... c(NBC)]'
% les coeffs sont equivalent au cepstre complex !

if nargin==1,
 NBC=length(a)-1;
end
 
c=zeros(NBC+1,1);
ahelp=[a; zeros(NBC-length(a)+1,1)];  % "prolongation" de a

for n=1:NBC,
 som=0;
 for k=1:n-1
  som=som+k/n*c(k+1)*ahelp(n-k+1);
 end						% end de for k
 
 c(n+1)=-ahelp(n+1)-som;
end					% end de for n  
