function x = quadraticRoots(a,b,c)
% quadraticRoots  Compute the two roots of the quadratic equation
%
% Input:  a,b,c = (scalar) coefficients of the quadratic equation
%                  a*x^2 + b*x + c = 0
%
% Output:  x = vector of the two roots
%
d = sqrt( b^2 - 4*a*c);
x1 = (-b + d)/(2*a);
x2 = (-b - d)/(2*a);
x = [x1, x2];
end