function B=interpol1(A);

%        function B=interpol1(A);
%
%Si A est une matrice colonne donnée par un vecteur v, prenant la valeur 0
%sur un sous-ensemble E, cette fonction interpole v par le vecteur obtenu
%en interpolant aux points de E les valeurs prises sur le complémentaire de
%cet ensemble. On peut l'utiliser en particulier pour remplacer un vecteur
%réel v par le vecteur obtenu en interpolant les maxima locaux de abs(.
%Si A est une matrice, l'opération est effectuée sur toutes les colonnes de
%la matrice. Cette routine constitue ici un sous-programme dans gaussq. 

mm=size(A);
B=zeros(mm);
for l=1:mm(2)
   v=A(:,l);
   x=find(~v==0);
   lv=length(v);
   lx=length(x);
   t=zeros(lv,1);
   for kx=1:lx-1
      vc=v(x(kx));
      vp=v(x(kx+1));
      delta_x=x(kx+1)-x(kx);
      t(x(kx):x(kx+1))=retta(vc,vp,delta_x+1);
   end;
   B(:,l)=t;
end   