function m=retta(y0,y1,np);

% m=retta(y0,y1,np);
%
%   function m=retta(y0,y1,np);
%Réalise l'interpolation linéaire entre 
%y0 and y1 avec pas np. Cette routine constitue ici un sous-programme 
%dans interpol1, donc dans gaussq. 

c_ang=(y1-y0)/(np-1);
v=[];
for x=0:np-1
   v=[v y0+c_ang*x];
end;
m=v;