Add beta-equilibrium ye(rho) *gnu script