function U=TraFouRap(u,n)
if n>0 then
   N=2^n
   w=exp(-2*%i*%pi/N)
   P=zeros(1,N/2)
   for i=1:(N/2)
      P(i)=w^(i-1)
      end
   I_pai=1:2:N
   I_imp=2:2:N
   i_inf=1:(N/2)
   i_sup=(N/2+1):N
   U=ones(1,N)
   U(I_pai)=TraFouRap(u(i_inf)+u(i_sup),n-1)
   U(I_imp)=TraFouRap(P.*(u(i_inf)-u(i_sup)),n-1)
   else U=u
   end

