function [x1 x2] = quadform2(a,b,c)

d = discriminant(a,b,c);
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);

end

function d = discriminant(a,b,c)
d = sqrt(b^2 - 4*a*c);
end
