function f=R4etoile(x);

%                  f=R4etoile(x);
%
% Action de l'opérateur adjoint de R4 sur un vecteur ligne x. 

N=length(x);
H=zeros(1,8);
H(1)=.230378;
H(2)=.714846;
H(3)=.630881;
H(4)=-.027984;
H(5)=-.187035;
H(6)=.030841;
H(7)=.032883;
H(8)=-.010597;
ff=zeros(1,2*N+2);
for k=9:2*N,
	if (mod(k,2) == 1)
		ff(k)=H(1)*x((k-1)/2)+H(3)*x((k-3)/2)+H(5)*x((k-5)/2)+H(7)*x((k-7)/2);
	else
		ff(k)=H(2)*x((k-2)/2)+H(4)*x((k-4)/2)+H(6)*x((k-6)/2)+H(8)*x((k-8)/2);
	end;
end;
tt=1:2*N-7;
f(tt)=ff(tt+9);