function x = roots(A, B, C)
% roots of the quadratic
x(1) = (-B + sqrt(discrim(A, B, C))) / (2 .* A);
x(2) = (-B - sqrt(discrim(A, B, C))) / (2 .* A);