Loh's method for quadratic formula

altRoots.m

% https://www.poshenloh.com/quadratic/
function rts = altRoots(p2)
    if p2(1) ~= 1
        p2 = p2 / p2(1);
    end
    
    average_value = p2(2) / 2;
    product_value = p2(3);
    
    distance = sqrt(average_value^2 - product_value);
    rts = -average_value + distance * [-1; 1];