quad_pos_solve function

Solve for positive root of quadratic ax^2 + bx + c = 0 when it exists