function b_root=steady_helper(coef_a,coef_b1,coef_c)

p = [coef_a coef_b1 coef_c];  % set up a polynomial ax^2+ b1*x + c = 0
rts = roots(p);   % find the quadratic roots
b_root = max(rts);  % pick the bigger one
