Solve for positive root of quadratic ax^2 + bx + c = 0 when it exists
quad_pos_solve(a, b, c)
a
coefficient of x^2
b
coefficient of x
c
scalar coefficient
real positive root if it exists