function a=cepst_to_a(c,N);
% Syntax: a=cepst_to_a(c,N)
%      ou a=cepst_to_a(c,N)
%
% c est vecteur colonne de la dimension NBC+1: [0 c(1) c(2) ... c(NBC)]'
% a sera vecteur colonne de la dimension N+1: [1 a(1) a(2) ... a(N)]'
% N est le nombre de coefficients LPC desires (default NBC)
% N ne peut pas etre > NBC

if nargin==1,
 N=length(c)-1;
end
 
a=zeros(N+1,1);
a(1)=1;

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