I have this computation for singly reinforced rectangular beam and found that
rho_max = (3/7) * 0.85 * f'c * beta_1 / fy
That is based on ACI 318 strain of tension steel of 0.004 and concrete of 0.003.
The same formula actually was posted by user JAE (Structural) in this thread...