function f=radonatt_aux(A); 

%Routine auxiliaire pour radonatt et invradonatt

[m,m]=size(A);
f=zeros(m);
for i=1:m
    f(i,:)=sum(A(i:m,:),1);
end