function [rho]=steady_state_rho(betta)
rho=-log(betta);
end

