function  R  = racines(a,b,c)
%
% La fonction racines
%
R=zeros(size(2,1));   % initialisation
delta = b^2-4*a*c;   % calcul du discriminant
if((delta>=0) &(a~=0))
    R(1) = (-b +sqrt(delta))/(2*a);
    R(2) = (-b -sqrt(delta))/(2*a);
else                  % racines complexes
    R(1) =NaN;
    R(2) =NaN;
end