equ2bis(a,b,c)

function [x1, x2, nb] = equ2bis(a,b,c)
%  EQU2BIS rsolution d'une quation du second degr version relle scalaire
%  function [x1, x2, nb] = equ2(a,b,c)
%                  2
%  rsolution de ax + bx + c
%  x vecteur des racines
%  nb nombre de racines
if nargin ~= 3
   error('equ2 demande trois arguments');
else
   if a ==0
      [x1,nb] = equ1(b,c);
      x2 = [];
   else
      delta = b.*b - 4.*a.*c;
      if delta > 0
         delta = sqrt(delta);
         if b > 0
            x1 = (-b - delta)./(2 .*a);
         else
            x1 = (-b + delta)./(2 .*a);
         end;
         if x1 == 0.0
            x2 = 0.0
         else
            x2 = (c./a)./x1;
         end;
         nb = 2;
      else
         nb = 0;
         x1 = [];
         x2 = [];
      end;
   end;
end;